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EXECUTIVE SUMMARY 



Remote sensing techniques applied to amiospheric monitoring problems offer a spatial and temporal database 
coverage which cannot be duplicated by conventional ground based measuring programs. However, in spite of these 
obvious advantages, the methods of interpretation and inversion of remotely sensed optical data are often more 
qualitative than quantitative. The sensitivity of existing techniques for the inversion of satellite imagery data is not 
well understood in terms of discriminating constituent types and threshold levels of detection. This lack of technical 
support has retarded the defmition of a comprehensive strategy for satellite sensor development geared towards the 
remote sensing of tropospheric aerosols. 

The objective of the present study was to simulate the atmosphoicaUy scattered signal received by a satellite sensor 
and hence to evaluate the sensitivity of this signal to atmospheric aerosols in general and to atmospheric pollutants 
in particular. The emphasis of the investigations was on broacfttand scattering and continuum absorption effects (Le. 
the level of information which can be extracted from typical remote sensing satellites) rather than the high spectral 
resolution properties characteristic of structured molecular absorption phenomena. In the long term this type of 
research is geared towards the development of an operational satellite monitoring system within the context of the 
long range transport of atmospheric pollutants (LRTAP) programs. 

The simulations were performed as a function of environmental, geometrical and operational parameter constraints. 
The environmental modelling included the effects of urban and rural type aerosols where the absolute and relative 
concentration of water soluble (sulfate) particles was allowed to vary and where particle growth and refractive index 
effects of relative humidity were incorporated. Combined aerosol modes which included fme particle (sulfate based), 
coarse particle (quartz based), and absorptive aerosols (carbon based) were incorporated as external (independent) 
aerosol mixtures. Some tests were performed in terms of internal (colayercd) aerosol mixtures in order to evaluate the 
effect of the mixture assumption on the results of the sensitivity study. 

The variability of the simulated signals indicate that information on the vertically integrated contribution of 
submicron scattering and absorbing aerosols can be useAiIly extracted. Inasmuch as the dominant optical variations 
occur in the lower troposphere this extracted infcmnation can be correlated with ground level concentrations. 
Furthermore, investigations into the dependence of the satellite signal on vertical inhomogenity (altitude dependent 
mixing ratios) indicated a degree of insensitivity which significantly simpliBes the radiative transfer calculations and 
the inversion methodology. 

The optical effects of ground level sulfates (predominandy in the form of ammonium salts) are strongly influenced 
by the ambient relative humidity. Such induced variations in sulfate optical activity result in a ground level mass 
scattering efficiencies which may vary by a factor of 4 or S across the naturally occurring range of relative humidity. 
Furthermore the scattering contribution by non ailfate fine panicle constituents yields scattering efficiencies (per 
unit mass of sulfotes) which, although remarkably consistent, may vary significandy from one region to the nexL 
Because large angle aerosol scattering effects are dominated in the visible spectral region by Rne particle aerosols and 
because fine particle scattering is predominantly a volume or mass effect die variation of the mass scattering 
efficiencies due to size changes in the dry dimensions of the fine particles is minimaL 

The implicadons of the former variations in tenns of satellite remote sensing is that one must, within the context of 
achieving a sulfate measuring precision of ± S ]ig/m}, eithor directly measure the ground level mass scattering 
efficiencies or develop a methodology which permits their computation from ground based meteorological 
measurements combined with parameters txtnuctai from the satellite data. 

Fine particle mode (submicron) particles can be distinguished from coarse particle (supermicron) particles by means 
of two or more bands placed in the spectral region between 0.4 and 1.5 |im. The specificity of sulfates within the 
fme pardcle mode is more a function of the dominance of fine particle (water soluble) aerosols in terms of optical 
effects and in turn the predominance of the sulfate ftaction in the water soluble component. At ground level the 
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optical contribution due to sulfates is typically the largest particularily in conditions of high relative humidity and 
total extinction. Accordingly, even in the absence of information relating to the contribution of other components 
(but with information on relative humidity) one can extract an estimate of dry sulfate mass to within a factor of 2 or 
better. 

The study also showed that inversion algorithms should be restricted to reflectances well below and above the 
Kaufman critical surface reflectance. The critical surface reflectance condition, which is characterized by a 
counterbalance between atmospheric backscattcr and attenuation of surface reflected signals, induces inversion 
singularities since the satellite signal is effectively transparent to the mass loading of scattering particles. It was also 
found that the satellite signal was fairly insensitive to the presence of carbonaceous soot over small surface 
reflectances (scattering atmosphere) but conversely was strongly dependent over high surface reflectances. This 
suggests that an inversion methodology would be best served by formulating the sulfate inversion around low 
reflectance pixels in the image and the soot invasion around high reflectance pixels. In the former case one would 
consider water pixels and vegetation pixels (in the chlorophyll absorption band near 0.65 [m) while in the latter case 
one could employ for example snow or bare soil surfaces. 

The extraction of mass loading parametere from die total signal received by the satellite sensor implies some apriori 
knowledge of the ground surface reflectance in the satellite image. To a certain degree this requirement can be 
circumvented by restricting measurements to spectral regions and scene surfaces whose reflectance is small. However 
to fully exploit the potential of the satellite data requires a calibration program involving the establishment of 
reflectance calibration sites and the multi-temporal analysis of repetitive satellite images. 

TTie operational parameter constraints were keyed to the specifications of selected remote sensing satellites which are 
eithor operational or will be by at least 1995. These included the NOAA and GOES meteorological sateUite scries, 
the Landsat and SPOT resource monitoring satellites and the oceanographic satellite SEAWIFS. The meteorological 
satellites offer a relatively high rate of image repetition over the same area but arc consffaincd to one or two coarse 
spectral bands and a large spatial resoluUon. The 10 bit NOAA sensors, although less than adequate with respect to 
the number and resolution of spectral bands, offer the best dynamic range and radiometric resolution characteristics 
of existing spacebome sensors. The current SPOT and Landsat satellites acquire high spatial resolution images (order 
of tens of meters) at comparatively high spectral resolution but with a coarse repetition cycle (order of weeks) and a 
radiometric resolution which is marginal for detecting differences induced by significant changes in sulfate 
concentration. 

The Landsat SEAWIFS sensor, due to be launched in 1991, offers a specttal resolution, repetition rate (2 days), an 
unagc overlap capability and a radiomeuic resolution which more closely a^Jroximatc die requirements for passive 
atmospheric remote sensing. Within the near future the HIRIS and MODIS sensors possibly coupled with vertical 
struciune data from the laser sounder (LAS A) represent the realization of a sensor package for which one of the major 
design constraints was the remote sensing of atraoq)hcric constituents. The latter sensors will be part of the EOS 
space aation package to be launched in 1995. 

A comparison of radiative transfer models indicated that die analytical 5S model developed by D. Tami and co- 
workers represented Uie most optimal compromise of speed and accuracy in terms of applying a pixel by pixel 
mversion aJgoritiim to large satellite images. Comparison with a multiple scattering, multiple layer model 
demonstrated however that, for non nadir conditions, modifications to the model were required in order to achieve an 
accuracy commensurate with more rigorous solutions. These modifications are relatively easy to impleanent, and add 
littie to the computation time. 

The simulation analysis was sufficienUy promising to warrant a program of simultaneous ground based, airbcmc and 
satellite measurements as a subsequent phase to the present work. This program should include as well the 
acquisition and analysis of high spectral resolution airborne data to address the specificity jroblem. 



vn 



Table 2.1.1.1 
Table 2.1.3.1 
Table 2.1.4.1 

Table 2.1.4.2 
Table 2.1.4.3 
Table 2.2.1.1 

Table 2.2.1.2 

Table 2.2.2.1 

Table 3.1 
Table 3J2 
Table 3.4.5.1 

Table 3.4.6.1 

Table 4.2.1.1 

Table 42.3.1 
Table 42.3.2 

Table 42.4.1 
Table 4.3.1 



LIST OF TABLES 

Log Normal Parameters 

Aggregate Aoosol Models 

Ground Based Measurements of SulMe Concentrations and Aeroscd Scattering 
CoefiHcients 

Maximum Surface Number Densities of Aggregate Aerosol Models 

Number Density Values (N) Used in the Radiative Transfer Calculations 

Cbmputed Optical Parameters (Particle Size Distribution (PSD) indqiendent of 
altioide) 

Refractive Indices of water and dry component aerosols which were employed in all 
optical calculations 

Sample refractive index calculations as a function of relative humidity (wavelength 
of O.SS pjn) 

Radiative Transfo' input/ouQHit parameters 

Radiative Transfer Calculations (input parameter variations) 

Apparent reflectance variation corresponding to an SO4- variation from to SO |Xg-m'^ 
normalized to the total satellite sensor noise oror tabulated in Table 42.4.1 

Observa's zenith angle and scattering angle in the solar plane for a solar zenith angle of 
4S degrees 

Table of measured ground reflectances and their uncertainties for a number of surface types 
derived from Landsat data (from Royer et al. 1988) 

Single scattering albedo as a function of wavelength, relative humidity and aerosol model 

Error rate (5p*/5o)o)f,p for nadir geometry, 45 degrees solar zei^ angle and 0.633 (un 
wavelength (units of percentage change pa- .01 incremoit in (Oo) 

Noise Equivalent Reflectances (N£R) and Quantitization Error Reflectances (QER) 

Input parameter precision required to achieve a sulfate concentration precision of 1.6 
^g/m^ per parameter 



Table 6.1 



Corrected Number Densities 



1. INTRODUCTION 



The principal objective of this study was the simulation of the atmospherically reflected signal received by a satellite 
imaging sensor and an evaluation of the sensitivity of this signal to atmospheric aerosols. Within this context we 
sought to evaluate the sensitivity of the signal to pollutant aerosols with an emphasis on fine particle sulfates. The 
motivation for an analysis of this nature is to evaluate the potential for employing satellite imagery as a means 
towards monitoring spatial and temporal variations of anthropogenic aerosols. Such an evaluation is essential 
towards understanding how remote sensing technology can be applied to enhance or complement operational 
measuring techniques employed in long range pollutant (LRTAP) problems. 

The sensitivity study was then employed to define an error budget fw an inversion procedure which sought to extract 
altitude integrated concentrations of sulfates given realistic values and uncertainties in both operational and 
environmental parameters. The operational parameters included radiometric specifications of a variety of existing and 
planned satellite sensors while the environmental parameters encompassed the solar, atmospheric and surface 
characteristics likely to have a significant effect on the total measured satellite signal. The incentive for this 
sensitivity analysis was to permit a more judicious selection of operational parameters and an understanding of 
o-adeoff considerations given different levels of uncertainty in the environmental parameters. 

1.1 Background ^ 

The exploitation of ground based optical measurements for the quantification of atmospheric turbidity and (at least 
implicidy) aerosol concentration variations has been carried out by astronomers since the early 1900s (Roosen et al., 
1973). In more recent times, meteorologists have recognized the utility of optical measurements as a means of 
quantifying the aerosol presence in the atmosphere (Middleton, 1952). 

One of the most apparent optical effects of aerosols is the attenuation of visible light reflected from structures at 
distances of Idlometers from an observer. This phenomenon which is predominantly due to the scattering of light out 
of the line of sight of the observer is to a great extent influenced by the fine particle (accumulation) mode in the 
radius size range between .1 and 1 ixm. Numerous studies have shown that ground level scattering coefficients and 
measures of visibility are strongly correlated with the sulfate component over urban and rural sites (see for example 
the references detailed in Table 2.1.4.1.), Other researchers have observed that scattering by carbonaceous products 
and ammonium nitrate can be important at least in urban areas (Sloane, 1983) while absorption by carbonaceous 
soot may significantly effect the light scattering budget in rural areas (Japar et al., 1986; Ackerman and Toon, 1981) 
and in particular over urban sites (Rosen et al., 1978; Ackerman and Toon, 1981; Sloane, 1983; Richards and 
Bergstrom, 1986). 

These types of optical correlations have inspired a few reseachers in the remote sensing field to investigate the 
applicability of satellite optical imaging as a means of monitoring at le^t the most optically important atmospheric 
pollutants. In fact the use of remote sensing for the detection of both anthropogenic and natural aerosols from space 
has, in the past ten years, begun to be recognized as a viable monitoring option (Lyons, 1980; Quenzel,1982; 
Chung, 1986). This approach represents a natural evolution of ideas since the remote sensing community has for 
some time been concerned with the estimation of aerosol optical effects as a means of predicting and eliminating the 
atmospheric contribution to surface reflected signals (Turner and Spencer. 1972; Ahem et al., 1977; O'Neill et al. 
1978: Gordon. 1978). 

The principal application of satellite detection methodologies to date has involved the detection of massive aerosol 
dust clouds where the optical effects are actually dominated by supemnicron particles (Kaestner et al., 1982; Legrand 
et al.. 1985; Takayama and Takashima, 1986; Deuze et al., 1988; Egan, 1986; 1988; Legrand et al., 1988; Tanr6 et 
al., 1988). Since no satellite sensor, designed specifically for downward looking aerosol detection exists, the analysis 
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of remote sensing techniques has been performed in terms of meteorological and earth resources satellites These 
investigations include work with NOAA data (Kindman et al., 1985; Freemantle et al.. 1986; Takayama and 
Takashima, 1986; Carson and Needling, 1977; Ahmad et al., 1988; Rao et al., 1988). GOES data (Griggs 1979- 
Fraser et al.. 1986; Tsonis and Leaitch, 1986; and Tsonis, 1987). METEOSAT (Kaestner et al., 1982; Deuze et al ' 
1988; Legrand et al. 1985,1988), and the Landsat earth resources satellite (Meckler et al., 1977; Fraser 1986- 
Otterman et al., 1982; Royer et al., 1988; Tam€ et al., 1988). In the Soviet Union results have been reported for the 
Meteor satellite (Ivanchik et al., 1985). 

In terms of the fine particle mode a number of researchers have attempted to extract indicators related to sulfate 
loading from sateUite images (Lyons, 1980; Fraser and Kaufman, 1984; Tsonis and Leaitch, 1986; Tsonis. 1987) 
These investigations were primarily empirical in nature and sought to compare satellite derived estimates of sulfate 
content with a comparatively few ground based measurements. The simulation approach employed in this study 
permits a more flexible and fundamental understanding of the important variables in the error budget for extracting 
sulfate content. Provided the models employed are physically realistic and are based on realistic input data this 
technique will assist in the interpretation of empirical results while enabling a comprehensive evaluation of input 
parameter optimization. Even in those cases where the modeling is inexact in an absolute sense the relative sulfate 
sensitivity as a function of operational and environmental input parameters wUl provide informauon useful for the 
development of an inversion methodololgy. 

i.2 General Approach 

The simulation approach adopted in this work entails the development of models which are sufficiently realistic to be 
applicable on a general scale without being overly encumbered in unnecessary detail. Whether a particulate model 
atmbute can be considered a detail depends ultimately on its contribution to the satellite measured signal relative to 
the noise level of that signal. 

Initially the modelling efforts included the definiuon of aerosol size distributions for aerosol constituents known to 
be opucally active m the tropospheric boundary layer. In this regard we were aided by existing climatological models 
denved from extensive optical and in situ measurements. As well airborne particle size distribution measurements 
made dunng different seasons over rural and urban regions in southeastern Canada were employed in an effort to 
obtam regionally representative p^ticle characteristics. 

Aggregate aerosol models consisting of three externally mixed (independent) particle size modes (carbonaceous soot 
water soluble and dust like) were used as standard inputs throughout the simulations. The hygroscopic component of 
these aggregate aerosol models was assumed to be primarily due to the narrow fine particle water soluble mode which 
itse f was assumed to consist primarily of sulfate particles. The effects of relative humidity were then incorporated by 
applying a panicle growth model to the sulfate component of the fine particle mode and computing the modified size 
distnbution. "^^ 

Standard values of refractive index were employed in onler to define the optical characteristics of each particle size 
mode. In die case of the sulfate component the refractive index was computed as a volume weighted mean of the dry 
sulfate refracuve mdex and that of pure water. Assuming spherical panicles, a Mie scattering code was utilized to 
compute the scattenng and absorbing coefficients characterizing a single aggregate aerosol particle (or the additive 
effects of a number of such identical particles) and the relative angular scattering distribution (phase function) as a 
function of relative humidity. ■^ j "^ a 

These point volume optical parameters were then effectively integrated or averaged over altitude with the assumption 
that the aggregate aerosol particle size distribution was itself independent of alUtude. The columnar optical parmeters 
so obtained were then entered into a radiauvc transfer equation (5S) which had been preselected on the basis of speed 
and accuracy consideraUons. The 58 model was utilized to simulate the signal received by a satelUie sensor owa 
vanety of operauonal and environmental parameters including columnar sulfate concentration, sensor passband 
surface reflectance, aggregate aerosol model, carbonaceous soot absorption and relative humidity 



I 

I 



Given the results of the radiadve transfer simulations, the partial derivatives {error rates) relating changes in sulfate 
concentration to changes in the opo^onal and environmental parameters were computed. The uncertainties in these 
latter parameters were then evaluated in order to obtain an error budget for sulfate concentrations computed from 
satellite optical data. This error budget was analyzed in order to characterize those operational and environmental 
parameter ranges which ven optimal for the development of satellite based inversion methodobgy for the extraction 
of sulfate concentration. 
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CHAPTER 2 

AEROSOL MODEL DEFINITIONS AND 
OPTICAL PARAMETER COMPUTATIONS 



2.1 Aggregate Aerosol Models 

In this section we characterize the component particle size modes and the manner in which they are combined in onler 
to define the two basic aggregate aerosol models (standard and urban). The aggregate aerosol model specifically refers 
to an average aerosol panicle defmed by the fractional number density conoibution of each component (see also 
Appendix A.l). 

2.1.1 Component Aerosols Defined 

Table 2.1.1.1 is a compilation of aerosol component model parameters which characterize the gross size distribution 
mode features of dry urban and continental aerosols. All the models employ the log normal size distribution as a 
conveniant and standardized means of specifying distinct modal contributions (c.f. Appendix A.l). Nearly all the 
models include three distinct modes (i) a water soluble mode whose optica] properties are. at least within the scope of 
this study, indistinguishable from sulfate properties, (ii) a soot mode characteristic of carbon based particles which 
are prcdominanUy responsible for aerosol absorption effects (Rosen et al.. 1978)) and (iii) a large particle mode "dust 
like" component which depending on its origin may or may not be hydroscopic (see Section 2.1.2 below). 

The entries labeUed SRA represent standard models defmed for the World Meteorological Organization by working 
groups concerned with aerosol modelling for cUmatic impact studies (WCP-55). The tabulations labelled SRA++ 
were an auempt by the working group to improve the somewhat unrealistic characterization of the water soluble and 
dust like components defmed in the initial (SRA+ in Table 2.1.1.1) tabulations. The Leaitch model parameters 
(Lcaitch and Issac, 1988) represent the results of a series of aircraft particle size measurements made at two different 
alutudes dunng different seasons over sites in Ontario. Quebec, New York state, Nova Scotia and the Northwest 
tcmtoncs. Finally the Shetde and Fcnn parameters were extracted from a report which dealt extensively with the 
optical effects of relative humidity (ShctUe and Fenn, 1979). 

For the simulations presented in this report we settled on a hybrid model which combined the water soluble and soot 
results of Leaitch and Issac (1988) with the dust like SRA definitions. The Leaitch and Issac models were employed 
principally because they represented size distribution data acquired over an area relevant to Canadian LRTAP interests 
(ajutheastem Canada). The universality of Leaitchs data and the differences between his model paramctm and those 
of the SRA model are points worth considering. However for our purposes we merely note that a remarkable 
consistancy was obtained in the fine particle mode of Leaiich's data over a diversity of sites. Similar results in terms 
of the value of rN were obtained by Kim ct al. (1988) over a pan of central USA but for less monodispcrsive 
aerosols (a in the order of 1.8). 

The fme particle accumulation mode, which in our three component model described above includes the water soluble 
mode and the soot mode, is (at least in the visible spectral region) the primary aerosol component responsible for 
the large angle scauenng and attenuation effects which influence the satellite signal. Sloane (1983 1984 1986) has 
noted that the fine particle mode consists primarily of sulfates and nitrates in the form of ar^moni'um salts ( 
(NH4)2S04 and NH4NO3 ). carbonaceous soot, organic carbon and residue which includes soU derived particles. 

Within the context of optical scattering properties, one can justifiably argue that the fine particle mode is composed 
primarily of sulfates. Although this concept is not essenUal to the development of a remote sensing inversion 
methodology it certainly simplifies matters when and if it applies. 



The ammonium salts are hydroscopic and can be regarded as the principal components of the water soluble mode. Of 
the remaining non hydroscopic components, only carbonaceous soot is distinctly different in optical properties from 
the ammonium salts. However, above the deliquescence point of the salts (sec Appendix A.2) hydroscopic particle 
growth effects increase their optical mass scattering efficiency relative to the non hydroscopic partilc^. Within the 
water soluble mode itself numerous authors have argued that the sulfate salts are more efficient scatterers per unit 
mass than the nitrate salts (Tang et al., 1981; Sloane, 1982). 

The implication of these assertions is that the optical scattering effects of sulfates are dominant whenever the mass 
of sulfates is comparable to the other fine particle mode components. In regards to mass distributions. Fraser and 
Kaufman (1984) have observed that, for the eastern United States, the ratio of sulfate mass (S04°) to the other fine 
particle mode constituents falls within the range 0.4 - 1.5. Furthermore, numerous optical scattering studies over a 
variety of meteorological conditions have demonstrated a predominant correlation between generic aerosol scattering 
and sulfate mass (Husar and Patterson, 1980; Pierson et aL, 1980; Macias et al., 1980; Leaderer and Stolwijk. 1980; 
Fraser and Kaufman, 1984; Vossler and Macias. 1986). 

Accordingly within the scope of the work presented below the fme particle mode and thus the water soluble mode tat 
identified with sulfates. This argument, it should be noted, applies to information extracted from scattering 
measurements in Uie visible spectral region and to externally mixed aerosol (see Section 2.1.3 below). 

2.1.2 Particle Size Dependence on Ambient Relative Humidity 

An important influence on the optical activity of aerosol particles is the ambient water vapour concentration or 
relative humidity which acts to cause deliquescence of hydroscopic aerosols. This deliquescence produces changes in 
intrinsic optical properties associated with a change in physical dimensions of the hydroscopic nuclei. The changes 
in intrinsic optical properties include a water volume induced decrease in refractive index along with an increase in 
scattering power due to relative humidity induced particle growth. 

In Appendix A.2 we describe a simple model, based on the standard log normal size distribution, which accounts for 
the optical effects of ambient relative humidity. Figure (2.1.2.1) shows some results for the particle growUi model in 
the form of size disuibution plots for the Lcaitch water soluble model at different relative humidities. Since tiie total 
number of particles remains constant the area of the curve dN / dlogr is fixed while the volume curve dV / dlogr 
grows in magnitude with increasing relative humidity. 

The effects of relative humidity on the other aerosol components were neglected. This is a reasonable approach given 
that the soot particles as well as the soil derived component of supermicron particles are non hydroscopic (Sloane. 
1983). Hydroscopic supermicron particles such as NaCl marine salts (Hindman et al.. 1985) were assumed to 
contribute negligibly to the optical effects in a continental or urban atmosphs^. 



Table 2.1.1.1. 
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Standard deviation 
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Soot 5 12 
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(internally mixed) 



Leaitch et al. (1988) * 

.138 
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1.24 

1.24 
1.24 



Shettle and Fenn (1979) •* 



Rural water soluble .027 

Rural dust like 4$ 



.19 
5.5 



2.24 
2J1 



Urban water soluble .025 

Urban dust like M 



.16 
5.1 



♦ + 



2.24 
2.51 



standard model WCP-55. (Table 22 in tack Appendix A) 

WCP-55; (p. 31 of their Table 4.5 which replaces their Table 2.2 of Section 2) 

Average of his 5 sites (.85 S r S 1 2 inn), 0% r.h. 

0%rJi. 

Personal communication (Leaitch. 1988); soot like water soluble centered around 3m diameter 
See Appendix A.1 for discussion of log normal parameters 



Note: 1. In rv = In IN + 3 In 2 c for log normal distribution 2. loga = stand, dev. of log(r/rN) ( or log(rAY) ) 



2.1.3 Dermition of Models Employed 

In order to understand the optical effects of sulfate concentration variations over a realistic and representational range 
of atmospheric aerosol models one must make a judicious choice of standard models. These "aggregate aerosol 
models" should, according to our viewpoint, permit investigations into (a) the effects of relative humidity on the 
combined contributions of hydroscopic and non hydroscopic aerosols (b) the optical effects of internal and external 
mixing and (c) the effects of varying degrees of correlation between the sulfate component and the other aerosol 
components. 

Two fundamental aggregate aoxjsol models were defined to represent (A) typical urban and (B) continental conditions 
(Table 2.1.3.1). The volume fraction parameters for the dry vCTsion of these models (0% relative humidity) were 
extracted directly from the modified SRA standard models described in the appendix of the WCP-55 report (see 
Table2. 1 .3. 1 and the WCP reference below). The log normal parameters employed for each of the three components 
of the aggregate models (Table 2.1.1.1) were actually a combination of the the SRA parameters (for the dust like 
component) and data based on the findings of Leaitch et al. (1988) for the water soluble and soot components. 

For each fundamental aggregate aerosol model the effects of relative humidity variations w«e considered in terms of 
6 nominal values (0, 50. 70. 80. 90. and 95 %). Computations performed for the 0% case were intended to play the 
role of standard calculations, independent of the relative humidity particle growth model, which could be more 
directly linked with other work based on the dry particle parameters. The other relative humidity values represent a 
realistic range of values in the continental and urban regions of eastern North America 

Questions relating to the manner in which the component aerosols physically coexist is the basis for considering the 
externally and internally mixed models listed in Table 2.1.3.1 Qabelled I and n respectively). For both the 
continental and urban aggregate models we investigated their optical variations in terms of component aerosols 
which could be treated as (I) independent constituents (externally mixed ) or (II) as components which could coexist 
as the same particle (internally mixed). In the latter case (which is treated in Appendix C) the aggregate aerosol was 
imagined to consist of a soot core surrounded by a water soluble mantle (and further surrounded by a water mantle in 
the case of non zero relative humidity) and a distinct dust like component However with the exception of Appendix 
C we employed the exiemallv mixed interpretation for all opti cal and radiative transfer calculations performed in this 
work. 

The parameters of Table 2.1.3.1 are fractional coefficients which fix the values of unitless optical parameters 
(refractive index, phase function, single scattering albedo etc.). Cumulative effects which depend on particle number 
density (number of particles per cm3 as described in Appendix A.l) ultimately affect the radiative transfer 
calculations and hence the magnitude of the signal computed at satellite altitudes. If one maintains the number 
density or volume fractions defined in Table 2.1.3.1 then variations in the number of each type of particle are 
significantly constrained in that component number density concentrations are perfectly correlated. Peiturbatioiis 
from this rather rigid model were considered by pemitting the water soluble component to vary independently of the 
other components. These departures from the standard urban and continental models defined of Table 2.1.3.1 (termed 
urban variant and continental variant) thus allowed a much more comprehensive analysis of optical effects due to 
variations in the water soluble component The details of the urban and continental variant mo(kls are given below 
in Section 2.1.4. 
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Tahlf 2.1.:;.1. Aggregate Afrnsol ModeLs 

I: Externally Mixed Case 

A) 



Contineatal Model 



Component / r.h. 


7i 


m 


50% 


70% 


Ci 


80% 




90% 


95% 


Water soluble model 
Dust like (dry) 
Soot 


.96566E00 

.10433E-2 
.33299E-1 


.29 

.70 
.01 


.3072 
.6830 
.0098 


.3533 
.6376 
.0091 




.4495 

.5427 
.0078 




.5715 
.4775 
.0060 


.6466 

.3484 
.0050 


B) 








Urban/industrial Model 






Component / r.h. 


7i 


0% 


50% 


70% 


ci 


80% 




90% 


95% 


Water soluble model 
Dust like 

Soot 


.73487E00 
.91668E-4 
.26504E00 


.61 
.17 

.22 


.6861 
.1368 
.1771 


.7760 
.0976 
.1264 




.8155 
.0804 
.1041 




.8821 
.0514 
.0665 


.9320 
.0296 
.0384 


II: Internally Mixed Case 


















A) 








Continental 


Model 








Component / r.h. 


0% 


50% 




ci 
70% 


80% 




90% 




95% 



Water soluble /soot .29/. 1 ♦ 
Dust like (dry) .70 



.3072/.0098 .3533/.0091 .4495/.0078 .5715/.0060 .6466/.0050 
.6830 .6376 .5427 .4225 .3484 



B) 



Urban/industrial Model 



ci 



Water soluble / soot .61/.22 .6861/.1771 

Dust like (dry) .17 .1368 



.7760/.1264 .8155/.1041 
.0976 .0807 



.8821/.0665 .9320/.0384 
.0514 .0296 



Note: 1 . Water soluble component is that of Leaitch and Issac (Table 2. 1 . 1 . 1) 

2. * "Insoluble fraction" for e.g. 3 .01/.3 (volume fractions same as for externally mixed model) 

3. The volume fraction coefficients (ci) are computed firom a knowledge of the number density fractions 
ifi) (c.f. Appendix A.I) and the log normal parameters derived from the particle growth models defined 
in Appendix A.2. 

4 . Dust like refers to SRA-h- model parameters (Table 2.1,1.1) 
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2.1.4 Surface Number Density Considerations 

(a) Extreme sulfate Variations 

Table 2.1.4.1 is a compilation of sulfate concentration and optical scattering coefficient measurements made at or 
near ground level over a variety of temporal and spatial parameters. The table lists (when available) extreme 
conditions along with averages over sampling periods which vary from days to months and over urban and rural 
geographical regions. Some of the listings include, as well, readings of relative humidity made in conjunction with 
the ground based measmem^ts. 

The extrema for the data compiled vary from 1 to 50 tigm-3 of (dry) S04= or from (1 to 50) x [132/96.1] dry 
CNH4)2S04. This ammonium salt form of sulfate is, according to many researchers (Sloane , 1983, 1984, 1986 for 
example), the dominant sulfate compound. Employing the notation of Appendix A.l we have for dry CNH4)2 S04: 

HN < 69 X 10 -12 (2-1) 

where p. is the average mass of an ammonium sulfate particle and N is the number density (cm-3). Accordingly: 

N < 69 X 10 -12 / [D (4/3) n rjS] particles per cm3 (2-2) 

where D is the density of ammonium sulfate particles ( approximately 1.8 gm/cm3; Vossier and Macias, 1986), and 
(4/3) n 13^ is the average volume of a log-normal distribution particle (c.f. Appendix A.l) 

For the log normal (Leaitch, water soluble) parametas defined in Table 2.1.1.1 we obtain 13 = .129 ^m and hence: 

Nmax = 4300 cm-3 



Equating the sulfate component to the water soluble component we proportion the particle size distribution according 
to the Yi values of Table 2.1.3.1 and hence obtain the following maximum number densities at ground level 
(extonally mixed model): 

Table 2.1.4.2 Maximum Surface Number Densities of Aggregate Aerosol Models 

Continental Urban / industrial 

# doisity [cm-3] # density [cm-3] 

4300 
0.53638 
1550.85 
5851.38 



Water soluble 


4300 


Dust like 


4.6457 


Soot 


148.28 


Total 


4452.93 



Tobli 2.1.4.3 
Numbir daniitu valuta (N) uatd In tha radlBtlva iranafar calculattona 
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(A) 



Corresponding voTgf ne scattertng coefflclentsig. f1Q*'«-4 m**-]] 

surface • densllu (N) 

[cm»»-3l [cm»»-3l NdnUlal)" delta N 'slgma at InrHrated relwHyn humtrtHi.c ff ) 

dBllBN NdnltlBl) Sigma 50 70 60 90 95 



CONTINENTAL VARIANT 



welersol. 
dust like 
toot 
total 

scat effec [m""2/g! 

URBAN VARIANT 



430 



430 



0,215732 0,257445 0,346665 0.410242 0,658605 1,184306 

2,3229 0,145383 ' 

74,14 0.038732 

76,46 0,184116 0.215732 0,257445 0.346665 0,418242 0,658605 1,184306 

4,257477 5.060693 6,841437 8,254007 12,99757 23,37228 



water sol. 
dust like 
soot 
total 



"^O 0.215732 0.257445 0,346665 0.418242 0,658605 1,184306 

0,26819 0.016765 

775,42 0.405097 

430 775,688 0,421882 0.215732 0,257445 0,346665 0,418242 0,658605 1,164306 



scat effec [m»«2/8l 
CONTINENTAL STANDARD 



4,257477 5,080683 5,841437 8.254007 12,99757 23,37228 



water sot. 


430 








dust like 


0,46457 








soot 


14,828 








total 


445,2926 









•cat effec [m"»2/fll 
URBAN STANDARD 



0,215732 0,257445 0,346665 0.418242 0.658605 1.164306 

0.029076 0.029076 0,029076 0,029076 0,029076 0,029076 

0,007746 0,007746 0,007746 0,007746 0,007746 0,007746 

0,252554 0,294267 0,383487 0,455064 0,695427 1.22H29 

4,984170 5,807377 7.568130 8,980701 13.72427 24,09897 



water sol. 
dust Ilka 
soot 
total 



430 





0,053638 





155.085 





585.1386 






•cat afftc [m""2/gl 



0.215732 0,257445 0,346665 0,418242 0,658605 1,184306 

0,003357 0,003357 0,003357 0,003357 0,003357 0,003357 

0,081019 0,081019 0,061019 0,081019 0,061019 0,081019 

0,300109 0.341622 0,431042 0,502619 0,742982 1. 268683 

5,922657 6,745663 6,506617 9,919188 14.66275 25.03746 



Notaa: 

•cat affac . changa In tha total voluma acat coaf par unit change In tha mass of dru aulfate 
croas •actions obtatnad bg moltlplulnB Mia affaclenclaa (Q) bu p)T2»»2 whara r2 U the affactlvi 
mean aurfaca radtua of a log normal dlatrlbutlon 



• Sample croas sections (stgma) [units of I0"< 
(wavelength of 0.555 micrometers) 



•10 cm»"21 



Model 


Oust llks. 


Soot 


Sigma tscat) 


625,071 


5,22423 


Sigma Cobs) 


233.51 1 


6.7545 



Water soluble (sulfate) 

01 • SOX 70X • BOX ■ 90X • gsx 

5,01703 5,9871 B.06I98 9,72656 15,3164 27,542 

0,18244 0,15912 0,12833 0,113436 0,085325 0,060371 
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b) Surface Number Density Variations 

In an effort to characterize the extremes of surface number density variation for sulfate like particles as well as the 
other constituents of the aggregate aerosol models we defined two different cases of surface number density variation 
for the continental and urban aggregate aerosol models discussed in Section 2.1.3 (Table 2.1.3.1). The detailed 
surface number density descriptions of these 4 different cases {Standard continental , urban continental, continental 
variant and urban variant ) along with sample optical parameter computations are given in Table 2.1.4.3. 

In the. first (variant) case we considered depanures from the standard aggregate okdsoI models by permitting the water 
soluble (sulfate like) constituent to vary independently. Accordingly the sulfate like number density was varied 
between and 4300 cm-3 (10 equal increments) while maintaimng the other constituents at half of the extreme 
values tabulated above. As a consequence the standard number density fractions defined in Table 2.1.3.1 were 
reobtained at a sulfate number density of 43(X)/2 = 2150 cm-3. This approach thus affords a distinctly different 
viewpoint wherein the optical effects of aerosols are dominated by sulfate variations which are totally uncorrelated 
with the other constituents. 

In the second (standard) case the number densities of the continental and urban industrial aggregate aerosol models 
were varied uniformly while maintaining the number density propcations Ti defined in Table 2.1.3.1. This approach 
thus assumes correlated changes between the number densities of all constituents at all relative humidities. In 
amplication the number densities were allowed to vary from to the maximum values defmed above in increments of 
one tenth of the maximum value. 

2.2 Optical Parameter Variations , 

2.2.1 Optical Parameter Computations 

Having defined the aggregate aerosol models one is in a position to compute the point optical (Mie) parameters at a 
given radius and wavelength and subsequently to integrate these parameters over the size distributions defined in 
Section 2.1 (Royer et al. 1985). Table 2.2.1.1 lists the required input parameters to the optical scattering 
computations along with the derived output parameters. The upper half of the table deals with point opdcal 
parameters which in principle represent a 3 dimensional volume function in the atmosphere but which in fwactice are 
permitted to vary only in the vertical direction or as assumed in this work, are independent of position. The second 
half of the table lists columnar parameters which represent vertically integrated analogs to the point qptical 
parameters. 

For a given radius and wavelength the Mie algorithms are employed to compute the scattering and attenuation cross 
sections and the phase function (c.f. Appendix A.l) for each aerosol component characterized by refractive index mi. 
The integration of the phase function over solid angle weighted by the cosine of the scattering angle yields the 
assymetry factor (g). This coefficient which is indicative of the anisotropy of the phase function is a required input 
to the SS radiative model (c.f. Section 3.1 below). Given the size distribution functions discussed in Section 2.1 the 
particle size distribution averaged optical parmeters are computed by performing a numerical quadrature across the 
optkMy active radius region for each aerosol component The aggregate aerosol models are then optically combined 
in terms of a fractional number density weighted avoage to obtain total average point optical parameters. 

Knowledge of the number density distribution and the radius dependent particle sized distributions as a function of 
altitude enables one to calculate the cumulative properties of each optical parameter for a complete atmospheric 
column. A significant simplification is obtained if it is assumed that the particle size distributions are independent of 
altitude and that the altitude profile (hi(z) in Table 2.2.1.1) is the same for all aggregate aerosol model components. 
The first assumption permits a direct link between cumulative parameters such as optical depth and point optical 
parameters (or equivalently between columnar number densities and ground concentrations) while the latter 
assumption guarantees a constant optical mixing ratio or hence a scattering phase function which is independent of 
altitude. Neither assumption is particularity valid in an atmosphere grossly divided into highly forward scattering 
Iropospheric aerosols and nearly isotropic molecular scattering above the troposphere. Fortunately however the 
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cumulative implications of these simplifications at the atmospheric boundary layers are not strong. In Section 3.3 
we discuss in more detail the level of enxjrs encountered when altitude dependent atmospheric parameters are replied 
by representative me^is. 





altitude) 








Type 


Parameter 


Symbol 


Derived optical parameters 


Symbol 


Point 


reCtactive index 


mi= mil -mi2j 


attenuation cross section 


o (mi, a) 




wavelength 


X 

A 


scattering cross section 
conservative phase function poi 


crsca (mi. a) 
(mi, ct, cos X ) 




unit particle size 
distribution 


ni(r) 


PSD averaged unit attenuation 
cross section 

PSD averaged unit scattering 
cross section 

PSD aveij^ed conservative 
phase function 


ai(X) 

Oi!^(X) 

poi (cos T ) 



Columnar 



number 






density fracdon 


Yi aggregate aerosol model parameters a(K) 


fraction by volume 


Ci 


Po(cosz) 
g* 


altitude profile 


hi(z) = Ni(z) / Ni(0) altitude optical profile 


Ki(z,X)/Ki(O.X) 


equivalent column 






length 


Hi qptical depths 


Ti = Ni(0)Hi<yi(X) 


surface # density 


Ni(0) 


tiSca„Ni(0)HiCTi««(X) 



t=N(o)XTi(o)a,(J) 
t^=N(o)X'^)<f(^Hj 



* asymmetry factor 



IS 



Table 12A.1 is a compilatton of standard refractive indices which we employed throughout this work to characterize 
the intrinsic dry optical properties of the water soluble, soot, and dust like components. The effect of particle growth 
as a function of relative humidity on the refractive index (water soluble case) was accomodated by computing the 
volume weighted mean of dry particle and pure water indices as-detailed in Appendix A^. 

Table 2.2.1.2 Refractive indices of water and dry component aerosols which were employed 
in all optical calculations 

wavelength 

(microns') water Water-Soluble Dust-Like ^SiSA 

0.400 1.339 -il.86E-9 1.530 -iO.O05 1.530 -i0.008 1.750-10.460 

0.488 1.335 -i9.69E-10 1.530 -iO.005 1.530-10.008 1.750-10.450 

0.515 1.334 -il.l8E-9 1.530-10.005 1.530-10.008 1.750-10.450 

0.550 1.333 -il.96E-9 1.530 -i0.006 1.530-10.008 1.750 -i0.440 

0.633 1.332 -il.46E-8 1.530-10.006 1.530-10.008 1.750 -i0.430 

0.695 1.331 -i3.05E-9 1.530 -i0.007 1 .530 -i0.008 1 .750 -i0.430 

0.860 1.329 -i3.29E-7 1.520-10.012 1.520-10.008 1.750-10.430 

1.536 1.318 -i9.97E-5 1.510 -i0.023 1.400-10.008 1.770-10.460 

2.250 1.292 -13. 90E-4 1.420 -iO.OlO 1.220-10.009 1.810-10.500 

3.750 1.369 -i3.50E-3 1.452-10.004 1.270 -10.011 1.900-10.570 

Source: Colwell (1983) 

2.2.2 Influence of Relative Humidity 

Figure 2.2.2.1 illustrates some typical variations of the real (scattering) and imaginary (absorption) parts of the 
refractive index as a function of the relative humidity for the Leaitch water soluble distribution. These examples 
calculated for the wavelengths ranging from 0.4 \im to 3.75 )J,m show the typical decrease in the real part of the 
refractive index from the dry particle value (Table 2.2.1.2) to the pure water value of less than 1.4. The somewhat 
intricate variation of the complex part of the refi^ctive index results firom the non monotonic spectral behavior of the 
dry panicle at 0% relative humidity and the subsequent transition with increasing relative humidity to the strongly 
monotonic (but comparatively small magnitude) water index. Table 2.2.2.1 shows explicitly some calculations of 
refiractivc index carried out for a wavelength of 0.55 (im. 

Table 2.2.2.1 Sample refractive index calculations as a function of relative 

humidity (wavelength of 0.55 ^m \ 

(Leaitch model, water soluble, Dry density = 1.8 g/cm^ (NH4)2S04 ) used in growth model) 



Relative 


Re&active index 


humidity 


rsLoan 


complex part 


0% 


1.530 


-i.006 


50% 


1.515 


-i .0055 


70% 


1.481 


-i.0045 


80% 


1.433 


-i .0030 


90% 


1.395 


-i .0019 


95% 


1.378 


-i .0014 
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Figiffc 2.2.2.2 shows the effect as a function of particle size parameter ( a = 2ra- / X ) of a change in relative 
humidity on the componenu of the scattering coefficient for the Lcaitch water soluble model (wavelength of 0.55 
\m). Tlie three graphs iUustrate weU that the increase in volume scattering coefficient with relative humidity is more 
due to the physical increase in the size of the particles than to a change of intrinsic optical properties. This follows 
from the expression for the differential contribution of a given particle size to the total scattering coefficient (ci 
Appendix A.I) viz: 

dKscat/da = Q(a)dS/da (2-4) 

where CKa). the scattering efficiency is plotted in Figure 2 J.2.2(a). dS / da the differential surface contribution due 
to pamclcs of size parameter a is plotted in Figure 2.2.2(b) and die product is plotted in Figure 2.222(c). Although 
some of the increase in moving from 0% to 95% relative humidity is due to a movement up the principal (first) peak 
of die scattering efficiency curve it is clear diat die large disparity in Figure 2J2.2.2(c) is more attributable to the 
surface magnitude disparity evident in Figure 2.2.2.2(b). 

In terms of die maximum number densities defined in Table 2.1.42 the volume scattering coefficient for dry sulfate 
particles at 50% relative humidity is given by: 

Kscai = Qrtr22N(0) ^5j 

= 2.6 X 10-* m-1 

where r2 = .126 jim is the mean log normal surface radius (Appendix A.1) calculated using tht Leaitch parameters of 
Table 2.1.1.1 and the scattering efficiency Q (= 1) was calculated by integrating die differential in equation (,2A) 
(Figure 2.2.2.2c) over size parameter. For a typical aerosol scale height of 1 km diis value of i£scat(0) then yields a 
maximum dry sulfate scattering optical depth of .22. 

If on die other hand we consider the wet sulfate at 95% relative humidity Uw integrated volume scattering coefficient 
IS enhanced considerably. Assuming die same maximum number density we obtain: 

K»«« = 11.8 X 10-4 

where r2 has now increased to J207 lun and die scattering efficiency is approximately 1.8. For a scale height of 1 
km the wet sulfate volume scattering coefficient corresponds to a maximum wet sulfate optical depdi of 98 This 
value IS consistant widi die optical depth magnitude reported in Tsonis & Lcaitch (1986). ITie volume scattering 
coefficient maxunums can be compared with the entries of Table 2. 1.4.1 (after multiplying die latter by 10). 

In Figure 22.2.3 we have plotted die variation of die scattering phase function ( po(cosx) defined in Appendix A.1) 
computed for die Leaitch water soluble model as a function of scattering angle and relative humidity and at two 
different wavelengdis. The fairly large increase of die forward (0°) scattering peak widi increasing relative humidity 
and die correspondmg decrease of backscatter at angles larger dian about 40 degrees is generally symptotic of bodi an 
mcrease m average particle size and a decrease in refractive index associated wiUi die increasing water content of die 
scattcrmg particles. It should be noted diat die implications of die phase function induced decrease in backscatter on 
die signal scattered toward die satellite sensor are generally outweighed by die corresponding increase in die volume 
scatienng coefficient. An increase in wavelengdi (0.55 nm to 1.536 m in Figure 2.2.2.3) effectively defines 
particles which are optically smaUer and more isotropic in dieir scattering properties. 

Figure 2 2.2.4 is a plot of single scattering albedo computed for die Leaitch water soluble model as a function of 
relative humidity for wavelengdis ranging from 0.4 ^lm to 3.750 \im. The consistentiy smaller values of single 
scaitenng albedo at die larger wavelengdis result because die scattering particles in diis wavelengdi regime are 
opucally small (average size parameter a = 2nr / X < \ ). For optically small particles Rayleigh (or dipole) 
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absorption and scattering phenomena dominate. Since Rayleigh absorption varies as a while Rayleigh scattering 
varies as a^ (Wickramasinghc. 1973; Kattawar and Plass. 1967) it is clear that an increase in wavelaigth will, 
depending on the actual value of the refractive index, induce an increasing dominance of absorption effects over 
scattering or hence a decreasing value of the single scattering albedo. At increasing values of relative humidity the 
single scattering albedo is larger because the imaginary part of the refractive index is decreasing (Figure 2.22.1) 
and/or because the oc^ variation becomes increasingly more important with increasing particle size. 

Figure 2.2.2.5 indicates the variation of the assymetry factor (g; defmed in Section 2.2.1) for the sulfate type 
particles as a function of relative humidity and at different relative humidities. This parameter is effectively an 
indicator of directional anisotropy in the scattering phase function (a value of zero indicates an isotropic phase 
function while a lai^e positve value is induced by a strong forward scattering peak). It illustrates well the tendeiKy 
towards a more isotropic phase furiction as (a) the waveloigth increases and (b) the relative humidity dexnases. Both 
effects which result at least in part &om an increase in t^tical particle size can be inferred qualitatively from the 
behavior of the phase function in Figure 2.2.2.3. 

Figure 2.2.2.6 shows the variation of the spherical albedo as a function of relative humidity and at incremental 
values of sulfate concentration (0.55 \xm wavelength). This parameter which is strictly not an intrinsic optical 
parameter indicates the effective (angularity integrated) reflectance of the atmosphere for an observer looking up. It is 
included because it demonstrates well the near linearity of atmospheric scattering with the mass loading of sulfates 
and because the classical non linear enhancement of scattering effects with increasing relative humidity is clearly 
evident 

2.2.3 Correlation of scattering and absorption coefficients with component 

concentrations 

(a) Point volume coefficients 

An important parameter in die analysis of the optical behavior of sulfates is tiie dry sulfate mass scattering efficiency 
given by: 

e(X) = AKScit(X) / AS04" [in2/g] (2-6) 

where ic*cat(X) is the total aerosol volume scattering coefficient [m-1] and SO4- is the mass density of dry sulfate 
particles [g/m^]. In terms of number density Qf) the dry sulfate mass doisity is given by: 

S04= = n D (4/3)jr r33 N (2-7) 

where T\ is the ratio of gram molecular weights of 504' and ammonium sulfiate ( (NH4)2S04 ) and represents the 
conversion factor to S04~ given that the sulfate particles are assumed to be in the form of aronwnium sulfate. The 
factor D is the density of a single ammonium sulfate particle (taken as 1.8 g/cm^) while the spherical volume term 
(4/3)r rs^ N is the total volume occupied by all the sulfate particles, assumed to be described by a log ncsmal 
particle size distribution with an average volume radius r3 {cS. Appendix A.1). 
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The entries of Table 2. 1 :4.3 are Mie computations of the component and total aerosol volume scattering coefficients 
as well as the sulfate mass scattering efficiency for different aggregate aerosol models (standard and variant surface 
number density cases) as a function of relative humidity. Employing the above expressicms with a value of r3 = .129 
\im (derived from the log nomial parameters of the Leaitch rtiodel in Table 2.1.1.1) yields: 

£(X) = AKScat(>.) / ( (96.1/132) x 1.6186E-14 [g] AN [cm-3] x 10<5 [m-3/cm-3] ) 

= 8.486E7 AKScat(X) / AN [m2/g] (2-8) 

where AN is taken as 430 cm-3 in Table 2.1.43. This relation yields the scattering efficiencies of Table 2.1.4.3. 

The range of sulfate mass scattering efficiencies listed in Table 2.1.4.3 are commensurate with ground based 
measurements such as those detailed in Table 2.1.4.1 . At a dry relative humidity of 50% the expression above 
produces values of t{X) ranging from 5.1 to 6.7 nfi/g. In comparison the in situ measurements of Table 2.1.4.1 
range bom 2.0 to a high of 30.0 with an apparent higher frequency in the neighbourhood of 10 - 15 m2/g. 

It should be noted that the values in Table 2.1.4.1 are open to interpretation in terms of the effects of relative 
humidity on the measured volume scattering coefficient As Kaufman and Fraser (1983) point out the measuring 
device (nephelometer) heats the scattering aerosols so that the ambient humidity is reduced with an accompanying 
decrease in the total aerosol scattering coefficient. They go on to say that the nephelometer based measurement can 
only be reliably used for relative humidities less than 70%. In only one of the cases for which we have 
documentation did the authors take pains to eliminate the heating effect and obtain scattering coefficients at the 
ambient relative humidity (Vossier & Macias, 1986). Except for these data then, one should generally regard the 
entries of Table 2.1,4.1 as estimates of the sulfate mass scattering efficiency at relative humidity values somewhere 
below 70%. 

The occurance of higher values of mass scattering efficiencies t(K) in the in situ measurements is in our viewpoint 
not due to the existance of larger fine panicle modes such as those reported for some urban atmospheres (Vossier & 
Macias; Hering & Friedlander, 1982) or less monodispersive fine particle modes such as those observed by Tang et 
al. (1981) These deviations from the assumed aggregate aerosol models, can add a significant contribution to the 
optical scattering signal but do not add significantly to the sulfate mass scattering efficiency (c.f. Appendix A.3). 
This reasoning along with the fact that these types of panicle modes were not measured in Leaitch's (1988) airborne 
program led us to exclude their incorporation in the aggregate aerosol models discussed in Section 2.1. 

The more likely reason why the in situ measurements appear to produce larger scattering efficiencies relates to the 
fact that our fine particle mode is assumed to consist entirely of sulfate plus a comparatively small amount of soot. 
In reality other fine particle aerosol components may contribute in the order of 30 to 90 % of the volume scattering 
coefficient (Table 2.1.4.1). The implication of a significant non sulfate contribution is clearly that the scattering 
coefficients calculated for a pure sulfate fine particle mode will be proportionately smaller than the nephelometer 
measured coefficients. This disparity between Mie calculations based on the sulfate component only and regression 
slopes derived from nephelometer measurements was pointed out as well by Vossier & Macias (1986; p. 1239). The 
excellent correlation coefficents which are typically obtained between the aerosol scattering coefficient and sulfate 
mass are in pan due to non sulfate particles which covary with the sulfates. 

Given a representative value for t(K) and a measurement of the aex)sol scattering volume coefficient k*c«K^) one can 
in principle calculate the sulfate mass density ^ig/m^] firom: 

S04= ^ 106 X ( t(X) )-l KScat(X) 0^ 

For the urban variant model of Table 2.1.4.3 a relative humidity of 0% yields .17 g/m2 for ( £(0.55 M-m) )-l. In 
comparison Fraser and Kaufman (1984) present a regression based expression which at 0% relative humidity and 0.55 
|im yields 0.10 g/m2 for e-l. As discussed above, the difference between the two values is probably due to the 
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covariance of other fine particle constituenis with sulfate (in addition to a non covarying quantity which they 
modelled as a covarying contribution with associated errors). 

(b) Vertically integrated coefficients 

If both sides of the point volume ei^ession above are integrated over altitude one obtains: 

IS04= = lO^x ( E(X) )-l x»c«(>-) (2-l(^ 

where ZS04= represents a mass density per unit cross sectional area of an atmospheric column [ng/m2] whUe 
Tsc«KX) is the total measured aerosol scattering optical depth. If the integrated sulfate expression is divided by 1000 
one obtains an equivalent sulfate surface concentration [M-g/m^] in terms of a forced aerosol scale height of 1000 
meters (the scale height is defined in Appendix A.1). The equivalent surface concentration represents the surface 
conceniradon required to yield the actual columnar concentration ES04= if the scale height of the aerosol mixing 
layer is fixed at 1 km. 

The relevance of this parameter in terms of the actual surface concentration is clearly related to the natural variation 
of the scale height and indeed whether aerosols can be realistically modelled in tarns of a mixing layer characterized 
by a exponentially decreasing number density. Elierman (1970) developed a vertical scale height model referenced to 
the horizontal visual range, which implied a very strong correlation between columnar quantities (optical depth) and 
surface quantities (visual range and volume scattaing coefficient) and which served as common departure point for a 
decade of radiative transfer woiic. His investigations, performed for horizontal visual ranges from 2 to 13 km, showed 
that scale heights ranging from .84 to 1.23 km could adequately represent the variation of aerosol volume attenuation 
coefficients between ground level and an effective upper boundary to the mixing layer at about 5 km. This spread of 
visual ranges, meant to encompass the so called haze regime between a clear atmosphere and fog, is clearly of 
particular imsest in sulfate sisveillance problems. 

The work of Isaac et al. (1988). however, indicated some rather significant variations in die altitude profile of the 
mixing layer with important implications on Uie computed scale height Kauftnan and Fraser (1983) presented results 
for visual ranges from 5 to 16 km (roughly corresponding to Elterman's haze regime) which also implied that 
integrated and surface measurements were often not well correlated. They point out that variations in scale height of 
the order of a kilometer produce large errors in vertically integrated aerosol mass estimates given that one has ground 
level measurements of visibility but no information on the scale height The analogue to this observation in terms 
of satellite based estimates of columnar mass loading is that large errors in the estimate of ground based 
concentrations can be expected if no information about the scale height is available. Given a large uncertainty in the 
mixing layer scale height, the relevant parameter extracted from satellite observations is thus the total integrated 
sulfate concentration (IS04=). However we will continue to use the concept of equivalent sulfate surface 
concentration in order to betto' relate our results to standard ground based pollution measurements. 
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2.2.4 Effects of Changes in Aggregate Aerosol Models and Equivalent Surface Number 

Concentration 

(a) External Aggregate Aerosol Model 

Figure 2.2.4.1 is a graph of total aerosol scattering optical depth as a function of equivalent sulfate surface 
concentration for all the aggregate aerosol models and surface number density cases listed in Table 2.1.4.3 
(wavelength of 0.55 jim). The computations, which were performed using the volume scattering coef^ctents listed in 
Table 2.1.4.3 and an effective scale height of 1 km, explicitly show the uncertainty in the total aerosol scattering 
optical depth ((^ alternatively the equivalent sulfate surface concentration) associated with the manner in which 
constituent number densities are combined to yield total number densities. By way of illastration the graph shows an 
uncertainty of approximately 10, 8 and 2 (ig/m^ sulfate for relative humidities of 0% 80% and 95% respectively if 
one has no information on which of the four models in Table 2J2.4.3 was employed to derive the volume scattoing 
coefficients. This error is in effect associated with with an uncertainty in e(X) (see also the discussion in Section 
4.2.2 (c)) 

In general the standard models yield larger aoosol scattering optical depths than the variant cases simply because all 
aerosol constituents are permitted to vary (as opposed to only sulfates) and hence the total number densities are larger 
at every sulfate increment The urban aerosol models yield larger aerosol optical depths than the continental models 
principally because of a larger soot component which contributes a small but significant portion to the total 
scattering. The dust like or large panicle contribution is not a significant factor in determining the gross visible 
wavelength scattering attributes of a particular model. 

(b) Internal Aggregate Aerosol Model 

The results ^ffesented in Appendix C indicate less of a dependence on external versus internal aerosol model than on 
the manner in which particle growth effects were modelled. Indeed the differences at the larger relative humidities 
were traced to the functional dependence of the wet to dry particle radii and in turn to the ratio of water mass to dry 
particle mass in terms of water activity. The mass ratio model which was employed for the calculations throughout 
this text (c.f. Appendix A.2) were found to give lowo' results than those of Appendix C as well as those-obtained by 
Tang et al. (1981). Presumably this difference relates to the fact that our model was based on an average aerosol 
while those of Tang et al. and Appendix C were stricdy derived for pure sulfates. In light of the rather dispanue 
results presented by Tang et al. (1981) the lack of difference between die external and internal model calcukticms 
bears ftirther investigation. 
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CHAPTER 3 

RADIATIVE TRANSFER SIMULATIONS 

Given the optical parameter variations typified by the results presented in Chapter 2 one is in a position to apply 
these numbers as input parameters to a radiative transfo model which effectively integrates the contribution from 
each point in the atmosphere to obtain the signal measured by the satellite. This problem is non trivial since solar 
photons incident on aerosols in the lower atmosphere may undergo several scattering events before being directed 
into the field of view of the satellite sensor. In addition the integrated effect (tf multiply scattered photons redirected 
into the sensor field of view is complicated by the distinctly different optical [HXipeTties of the upper atm<^here in 
comparison to the lower boundary layer whose properties we wish to investigate. One is then confit>nted with a 
classical radiative transfer dilemma of having to compromise between arbitrarily accurate but time consuming 
solutions and approximate solutions which rationalize the intrinsic optical parameters and hence reduce execution 
time. This question of computation time is further aggravated by considerations of the large size of remote sensing 
image databases and the specific requirement that an inversion must be perfonned on each point in the image. 

The principal objective of Section 3 is to demonstrate how the output parameter of the radiative transfer model 
(exoatmospheric reflectance) is influenced by the input environmental and operational parameters. Once this is 
understood we can proceed to the inverse problem of evaluating sulfate sensitivity in terms of the satellite signal and 
how uncertainties in the other input parameters will affect inversions which are performed. 

We briefly describe the input and output parameters required for the radiative transfer calculations as well as the 
choice of radiative transfer model employed The relevance of the optical parameters described in Section 2.2 is made 
apparent by indicating how these same parameters act as critical inputs to the radiative transfer computations. 

3.1 The 5S Model and Required Input/Output parameters 

The model employed for the radiadve transfer calculadons was the SS model (Appendix B). This model is a multiple 
scattering analytical solution to the radiative transfer problem posed in a plane parallel vertically homogeneous 
atmsophere. It's principal feanue is an impressively rapid execution time suited to the data intensive requirements of 
remote sensing. The limitations of the model which stem from the simplifying assumptions employed in its 
derivadon are discussed in Section 3.2. 

Table 3.1 indicates the operational and optical (environmental) parameters required as input to the SS model, (detailed 
descripuons of the optical input parameters are given in Appendix A.l). For each input vector of parameters one 
outputs an apparent reflectance: 

p* - tmol(e) { <p> T(e,t, g, ©o) TOs.t, g, (Do)/ [1 - <p> S(X, g, OOo) ] 

+ pt( ♦, e, 0», t, a)opo(cosx) ) ) (3-1) 

which is simply the upwelling radiance at the top of the atmosphere noimalized to the incident solar radiation on a 
plane parallel to the earth's surface (see Appendix B for example). The first term accounts for the attenuation of 
radiation reflected of the terrestrial surface while the second term is due to atmospheric backscattering only. The 
multiplicative transmission term (tmol(6)) models the effects of the most important molecular absorbers. All non 
geometrical parameters are in general a function of wavelength. 

One can perfomi computations of p* at a single wavelength or over a satellite passband by interpolating between 10 
rcfffcnce wavelengths (Table 3.2) for which fairiy precise radiadve transfer calculadons were carried out by the 
aulhcvs of 5S. The geomedical parameters include the observers zenith angle (9) measured relative to nadir and the 
observers azimuth angle measured relative to the solar azimuth ((}>). An illustrative diagram showing these angles as 
well as the solar zenith angle for a typical satellite geometry is included in Appendix B (Figure 2.1). 
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The effective target reflectance <p> is an area weighted average of the target reflectance (pO md the average 
background reflectance (pb). The parameter F in Table 3.1 is the weighting function which defines the relative 
influence of the target and the average back^und reflectance on the photons which have been reflected and 
subscqucnUy scattered. R defines the average size of the homogeneous target surrounding the pixel being viewed. 

The point optical volume attenuation parameters in generaUy enter into the radiative transfer equation as explicit 
variables. However in the case of 5S where one effectively assumes an constant optical mixing ratio (c.f. Section 
3.3) these parameters ento- indiracdy in tarns of their integrated analogue of optical dq>th. The scattering phase 
function ( po{cosx) ) which is a weighted mean of the molecular and aerosol phase functions is computed diiecdy in 
the molecular case and entered in tabular form at fixed scattering angles for each of the aerosol component modes. 
The assymetry coefficient g is computed from the phase function (ci. Appendix A.1 and Section 2.2.1). 

The single scattering albedo (o>o) represents the fraction of scattering to total energy in an elemental volume (K^cat / 
K). If one assumes constant number density mixing ratios with altitude its value is simply given by the ratio of 
scattering to total optical depths. The total optical depth t represents the sum of the aerosol scattering, molecular 
scattering, aerosol absorption and molecular absorption contributions. The molecular contributions (ozone, water 
vapour etc.) are computed for standard atmospheric models while the aerosol concributions vary with the type of 
aggregate aerosol model, the relative humidity and the surface number density of each aerosol component It should 
be noted as weU that absorption by molecular water vapour was tied to the relative humidity variations at a fixed 
average temperature. 

The spherical albedo (s) and total transmission term (T) represent cumulative atmospheric radiative ffansfer effects 
and are derivable in approximate analytical form from the radiative transfer equation. 



Table 3.1 Radiative Transfir innii t /oMtniif naram^fpr«^ 
Type Parameter 

Spectral wavelength 



Geometrical 



Atmospheric boimdary 



Point optical 



observer zenith angle 

observers azimuth relative sum 

solar zenith angle 

background reflectsKe 
target reflec&nce 
relative spatial albedo fimction 
effective target reflectance 

aggregate optical parameters 



Symbol 
X 



pt 

F(R) 
<P> 

kO.) 
KSca(X) 

Po (cos X ) 



Columnar parameters 

Molecular (absorptive) transmissioii 



T 
tfnol 



Output 



apparent reflectance 
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Table 3.2. Radiativt Tmiifgr ralcnlattniis Hiinut aaramrter varintinnsl 



Parameter 
Bands* 



Variations 

SEAWIFS 

SPOT-1 

TM 

MSSLandsat 
METEOSAT 
GOES 

AVHRR(N0AA9) 
(NOAA 10) 

0, (30). 180 



numbo' 

S 

I 

4 
4 
1 
1 
t 
2 



Comment 



ev 


(10) 70 


es 


30, 45, 60. 70, 75 


Pb 


Vegetation, water, sand 




(0. .2, .5) 


pt"pb=p 




F(R) 


5S mean 


a (X). a «a (X) 


Continental standard 


Po (cos X ). «»o (X) 


Cwitincntal Variant 




Urban / industrial standard 




Urban / industrial Variant 


t. I»c« 


Standard models 




Variant models 



i 

5 
3 

3 



standard 5S r^fectance model 
(low, medium and high lef.) 



tmol»c». tmol**'s 



6 aggregate aax)sol models® 
6x11 aggregate aerosol models 
6 aggregate aerosol models@ 
6x11 aggregate aerosol models 

1 1 increments of water soluble (sulMe) 
number density (c.f. Table Zl.4.3) 

1 1 incremoits of total aerosol number 
density (c.f. Table 2.1.4.3) 

SS standard values 



@ One model per relative humidity (0%, 50%, 70%, 80%, 90%, 95%) 

* The 5S computations are interpolated between the 10 reference wavelengths for which optical (Mie) parameters are 
required as input in order to obtain an effective continuum m wavelength. A quadrature sum is then performed over 
the satellite passband of interest The 10 reference wavelengths are .400, .488, .515, .550. ,633, .695, .860, 1.536, 
2250. and 3.750 (im. 
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Table 3.2 is a tabulation of the sets of input parameters employed in this study to investigate the variation of the 
exoatmospheric signal (apparent reflectance) as calculated by the 5S model. Apparent reflectances in different 
satellite bands are obtained by convoluting the 5S monochromatic (apparent) reflectances with spectral response 
functions stored in the 5S program. The geometrical angular parameters were chosen to obtain a representative range 
of values at a resolution which previous experience had shown was sufficiently dense to ensure reasonably linear 
segments between calculation points. For simplicity the target reflectance was chosen to be the same as the 
background reflectance. The background reflectances selected for investigation were fixed values between and .5 in 
addition to standard tabulations for low turbidity water, vegetation and sand. A host of aggregate aerosol models 
(defined in Section 2. 1 .4) were in addition employed to represent the range of agpegatc aerosol models which might 
be encountered in the lower atmosphere. 

3JL Model Dependency 

In a study complementary to this work the effects on the computed satellite signal due to differences in the radiative 
transfer models employed was evaluated (Royer ei al., 1988; Appendix B). Investigations into the optimal model 
which best combined the requirements of speed and accuracy indicated that the 5S model with certain qualifications 
was the best compromise for purposes of satellite data inversion. The limitations of this model are fairly well 
imderstood and it is sufflciendy fast that a simple iterative inversion procedure can be used to extract atmospheric 
optical depths and subsequendy sulfate densities. 

In brief the investigations showed that the simplifying assumptions incorporated in 5S induce a variety of small 
errors in the calculated apparent reflectance which become progressively larger with increasing solar or observers 
zenith angle, wiUi decreasing wavelength and as a function of the proximity of the actual ground reflectance to the 
critical ground reflectance (Section 3.4.3). These errors were thought to be primarily due to the manner in which 
aerosol and molecular scattering are decoupled in the model and as such were linked with the related factors of 
interparticle scatter between molecules and aerosols and the vertical inhomogeneity of the atmosphere. In later work 
we noted that the separation of these two phenomena was not obvious. It was observed that these errors could be 
considerably reduced by effectively transforming the 5S model into a two layer rather than a single layer 
representation of aerosol and molecular scattering. However in the interests of simplicity this correction was not 
incoT3X}rated in the present analysis. 

3<3 Effect of Vertical InbomogeDeities 

The kemal of the radiative transfer equation is the scattering phase function whose form changes with altitude 
depending on the relative contributions of the different aerosol and molecular phase functions (see Appendix A.1 for a 
description of the scattaing phase function and Appendix B for a general discussion of radiative transfer concepts). 
These relative contributions which are deteamined by optical mixing ratios of volume scattering coefficjents to total 
volume extinction are dominated by aerosols in the troposphere and by molecules above the troposphere. Bo:ause 
aerosol scattering is so strongly anisotropic in comparison to molecular scattering, the scattering phase function 
undagoes significant changes in form as a function of altitude. 

Since such vertical inhomogeneities have important implications in terms of the degree of complexity and 
computation time of a radiative transfer code it is imperative that one understands tiieir influence on the radiative 
transfer mechanisms in the atmosphere. In particular, is it reasonable to hypothesize that the radiation leaving the 
boundary layers of the atmosphere can be effectively evaluated using an equivalent radiative transfer model whose 
kemal phase function is a vertically averaged mean. Such an approach which clearly can only apply to lines of sight 
which view the entire atmosphere considerably simplifies the radiative transfer computations and hence the inversion 
methodology. 

The results of the study described in Appendix B corroborate the concept of equivalent mean atmospheres. Except for 
very large solar zenith angles or observer's zenith angle a single layer radiative transfer model computes apparent 
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reflectances to values better than .005 in moderately loaded atmospheres (aerosol optical depths less than 0.5). In 
cases where this may be insufficient (see the noise discussion in Section 4.2.4) simple bias relations referenced u> 
the single layer model can be developed. This latter approach is currently being investigated. 

With respect to accuracy the SS model performs less well than an arbitrarily accurate single layer model. This 
degraded performance is a non trivial result of the mmna in which the single layer characterisitics of this model are 
coupled with the analytical approximations employed in its derivation. 

3.4 Selected Results 

3.4.1 Variation of aggregate aerosol model 

Figures 3.4.1.1 and 3.4.1.2 are plots of apparent reflectance versus equivalent sulfate surface concentration for the 
four number density models (urban standard and variant and continental standard and variant as discussed in Section 
2.1.3) . surface reflectances of 0.0 and O.S and for two different relative humidities. These apparent reflectance 
variations illustrate the different effects of changing the number density model in the context of a predominantly 
scattering atmosphere (zero surface reflectance) and a predominantly attenuating atmosphere (surface reflectance of 
0.5). 

In the case of zero surface reflectance a maximum apparent reflectance difference of .006 occurs between the 
continental variant and urban standard models, at a relative humidity of 90% and at the maximum equivalent sulfate 
concentration of 50 iig/m^. At smaller relative humidities and surface concentrations this value decreases 
considerably (see also the discussion in Section 4.3 below). The lack of sODng relative variation between number 
density models is not surprising given the similar type of behavior observed in the aerosol scattering optical depth 
plots of Figure 2.2.4.1. 

However a second more subtle reason derives from the fact that the main difference between the models is the 
variation in the optical effects due to the soot component The addition of soot particles (as occurs for example in 
changing from the continental standard to the urban standard atmosphere) causes a roughly equivalent inaease in both 
the scattering and absorption volume coefflcients of this component The enhanced atmospheric backscattering effect 
of the increase in scattering coefficient is accordingly offset by the increase in absorptive attenuation. The 
consequence of this is that the change in aggregate aerosol model is nearly transparent to a satellite measurement of 
apparent reflectance. 

Such is not the case for a large surface reflectance where the attenuating properties of the atmosphere dominate the 
apparent reflectance measurement. Here the absonive and scattering increases due to the soot particles are actually 
combined since scattering in this context plays an attenuating role (outscattering between the ground and the observer 
in addition to outscattering <rf surface reflected sunlight which has be^ scattered one or more times into the receiver 
field of view). Counter balancing this attenuation effect, but of lesser importance, is the scattering contribution such 
as that observed in the zero albedo case plus the scattering of sur&ce r^lected radiation. The net effect is a decrease in 
apparent reflectance which is however less severe than the purely attenuating decrease given by Ap* /p* -dx (where 
Ax is the change in optical depth due to the change in aggregate aerosol model). 

One can observe also that the standard aggregate aerosol curves are considerably more linear than the variant model 
curves. This behavior would appear to stem from the fact that the variant model goes through a relative increase in 
scattering properties (odo increases) with increasing sulfate concentration. Note that since the same is true for the 
zero reflectance surface the effect must be more prevelant in the sur&ce reflected/scattered radiation. 
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3.4.2 Ground albedo variations 

The contribution to the total satellite signal of photons reflected at the ground and directly transmitted or rescattered 
into the sensor field of view is one of the single most important influences on the behavior of the apparent 
reflectance as a function of aerosol or sulfate concentration. Figure 3.4.2.1 illustrates this behavior for a variety of 
surface reflectance values and relative humidities, a nadir viewing geometry and a wavelength of 0.55 |im. 

The most notable effect of variations in surface reflectance is clearly the change in slope from postive to negative as 
the surface reflectance is increased. This phenomenon which has serious implications in terms of inversion 
algorithms for sulfate concentration is due to the competing influences of atmospheric attenuation and absorption 
(c.f. Appendix B). The reflectance for which the slope is actually zero (Kaufman. 1982) is indeed a singularity at 
which the two influences are perfectly balanced independent of atmospheric optical thickness. For surface reflectances 
below this critical surface reflectance the satellite signal is dominated by scattering effects and accordingly increases 
with an increase in the number of scattering particles (sulfate particles). I f the actual surface reflectance is greater 
than the critical surface reflectance than attenuation effects of surface reflected photons dominate and the satellite 
signal actually decreases with an toarease in the number of scattering particles. 

The critical reflectance singularity is best iUuscrated by a Kaufman type plot of apparent reflectance (less surface 
reflectance) as a function of surface reflectance. In Figure 3.4.2.2 the effectively linear segments which correspond 
each to a particular value of optical depth clearly intersect at a well defined critical reflectance. Since the position of 
the intersection point is a strong function of the average atmospheric single scattering albedo the singular behavior 
of the apparent reflectance in the neighbourhood of the point can be turned to advantage by utilizing it to quantify 
atmospheric absorption (see Section 4.2.3 below). 

However the key point to note is that the satellite signal is only weakly sensitive to atmospheric aerosol 
composition over those areas for which the surface reflectance is close to the critical reflectance. This point is further 
developed in terms of inversion errors for sulfate concentration in Section 4 below. 

3.4.3 Relative humidity effects 

Referring again to Figure 3.4.2.1 one can observe the significant effects of relative humidity on the computed value 
of apparent reflectance.over a range of relative humidities. The three ground reflectances (0, .2 and 3 in (a), (b), and 
(c) respectively) correspond to values wcU below, at and well above the critical ground reflectance and as such the 
variations in slope of the three curves with different ground reflectances illustrate well the tradeoff effects of 
atmospheric attenuation and scattering discussed in Section 3.4.2 above. 

At a surface reflectance of zero ( Figure 3.4.2.1(a) ) the increase in apparent reflectance with increasing sulfate is 
most significant when the relative humidity is high and the scattering cross sections of a given concentration of 
sulfate particles are maximum. At large reflectances above the critical reflectance ( Figure 3.4.2.1 (c) ) the effects of 
attenuation dominate and the larger cross sections associated with higher relative humidities effect a decrease in 
ai^arent reflectance as the relative humidity increases. 

At values of reflectance near the critical reflectance one expects a flat or nearly flat response curve for apparent 
reflectance versus sulfate concentration. Figure 3.4.2.1(b) does show this lack of sensitivity along with a weak 
minima for the larger relative humidities (note that the apparent reflectance scale is strongly magnified). The weak 
minima results from the fact that the Kaufman singularity is not well defined at the higher relative humidities. 
Figure 3.4.3.2 (relative humidity of 95%) illustrates this point and in comparison to Figure 3.4.2.2 clearly 
demonsuates how a value of reflectance chosen in the neighbourhood of the cluster of intersection points induces a 
non linear relationship between apparent reflectance and optical depth (sulfate concentration). 
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3.4.4 Absorption Effects 

Figure 3.4.4.1 shows the variation of apparent reflectance as a function of equivalent soot concentration. The 
insensitivity of the apparent reflectance to increasing soot absorption for zero surface reflectance is the same 
scattering/absorption offsetting effect previously discussed in Section 3.4.1. Sinrularily the comparatively large slope 
observed for high surface reflectance is indicative of the dominance of attenuation effects (Section 3.4.1). The sulfate 
and soot inversion implications of these variations is discussed in Sections 4.2.3, 4.2.5 and 4.3. Conversion of the 
soot concentrations to single scattering albedo values is discussed in Appendix A.3 (Figure AJ.3). 

Note that the soot concentration variation is arranged so that the urban standard and continental standard fractional 
number densities are obtained at 13.8 M.g-m-3 of sooL The soot concentrations reduce by a factor of .436 if the 
remarks of Chapter 6 are taken into account. 

3.4.5 Dependence on satellite bands 

Figure 3.4.5.1 illustrates the variation of the apparent reflectance (nadir geometry) as a function of equivalent sulfate 
concentration for a low reflectance target (5S lake standard) over a range of satellite bands located in the visible and 
near infrared spectral region. The magnitude differences between tte various curves is primarily due to the X-" type of 
falloff expected in the aerosol scattering optical depth (see Section 4.2.5b for example). Figure 3.4.5.2 and 3.4.5.3 
are plots of the same parameters but for vegetation and sand targets. The vegetation results pass from a typical 
scattering atmosphere response for visible bands to an attenuating atmosphere response in the near IR. The 
insensitivity of the sand target results because the sand reflectance is near Uie critical reflectance of the atmosphere. 
The reader is referred back to Section 3.4.2 for a more complete discussion of reflectance effects. 

A more relevant characterization of the low reflectance lake data (Figure 3.4.5.) is presented in Table 3.4.5.1 where 
the extreme variations of apparent reflectance, due to a change of equivalent sulfate concentration from to 50 
Mg/m3, are normalized to the noise equivalent reflectance and radiometric resolution reflectance (extracted from Table 
4.2.4.1) for each of the satellite sensors. A discussion of the implications of these signal to noise ratios is left to 
chapter 4. It suffices here to note ^t, for small surface reflectance, only those sensors with 10 bit resolution (or for 
which the dynamic range straddles the low reflectance range characteristic of atmospheric phenomena) can yield 
sufflcient on the signal variations due to sulfates. 
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TABLE 3.4.5.1 
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3.4.6 Iliumination Geometry 

(a) viewing geometry 

Figure 3.4.6.1 shows the variation of the apparent reflectance as a function of sulfate concentration for different 
values of observers zenith angle. The calculations were made fw a relative humidity of 95%, a surface reflectance of 
zero and the standard 45© value for the solar zenith angle. The left hand side of each graph corresponds to look 
directions on the anti solar side (observers back to the sun) while the right hand side of each ^ph represents the 
solar side. The distinction is best explained in terms of Table 3.4.6.1 below: 

Table 3.4.6.1 Observers zenith angle and scattering angle in the solar plane for a solar 

zenith angle of 45 degrees 

solar side anti solar gj je 

observers zenith 50 40 30 20 10 10 20 30 40 50 

angle 

scattering angle 95 105 115 125 135 145 155 165 175 185 195 

Note: the scattering angle is measured from the direction of the sun's rays to the direction of flic ray scattered into 
the sensor field of view. In the solar plane the calculation of the scattering angle can be performed by 
geometrical means. Off the solar plane the calculation is best performed using spherical geometry (see for e g 
Turner & Spencer( 1972) ). 

For zero surface albedo the apparent reflectance dependence on viewer geometry is predominandy a function of zenith 
angle and can be taken as being - 1/cose where 6 is the observers zenith angle (Deschamps et al. 1983). Second order 
depaitiffes from this first order variation are induced by the angular form of the scattering phase function. The 
scattering phase function, which for low order scattering considcnttions varies as the scattering angle from the sun, 
induces variations as a functicm of die observers zenith angle and the dbservci's azimuth angle. 

With respect to the shorter wavelength (0.55pm) variation in Figure 3.4.6.1 we note that the extreme end of the 
solar side of the soot phase function dominates (by about a factor of thr^) die extreme end of the anti solar side and 
that the phase function decreases monotonically between the two extremes. Since, for the urban variant aggregate 
aerosol model, the soot particles are optically dominant at zero sulfate concentration (see for example Table 2.1.4.3) 
the variation with l/cos9 on the solar side is much stronger than on the anti solar side wh«e the decreasing phase 
function acttially competes with increasing 1/cose. These differences are less ^parent at large sulfate concentration 
since die sulfate scattering phase function is much more symmetric about an observer's anith angle of degrees. At 
larger wavelengths (Figure 3.4.6.1 (b)) the aggregate aerosol scattering phase function is much more Rayleigh like 
(more isotropic) and die disparity between solar and anti solar directions is less obvious. 

Al large surface rcflt^tancc the apparent reflectance response curves are somewhat complicated by the competing 
effects of scattering and attenuation. At zero sulfate concenttBtion the l/cos9 and phase function variation mentiwied 
above exceed attenuation effects on die solar side while die opposite is true on the anti solar side. The l/cosG effects 
of scattering and/or die changes in single scattering albedo widi increasing sulfate concentration connibute to 
exaggerated non linearities of die apparent reflectance curves at die larger observer's zenith angles. 

Figure 3.4.6.2 shows in Uie case of large surface reflectance a more promising variation to consider is die derivative 
of die apparent reflectance widi respect to die cosine of die observer's zenidi angle. Using diis a^Jroach we can 
observe a comparatively well behaved monotonic increase as a function of sulfate concentration. It should be noted 
however diat die apparent reflectance variation is not large and accordingly implies a fine signal sensitivity. 
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(l») solar geometry 

Given that the apparent reflectance for zero surface reflectance varies approximately as lA;ose at fixed solar zenith 
angle (8,) one can appeal co the principle of reciprocity (Chandrasekhar. 1960) to argue that the apparent reflectance 
varies approximately as ^< . l/cosO, at fixed observer's zenith angle. This type of behavior is readily apparent for 
the zero reflectance case of Figure 3.4.6.3. The same argument of reciprocity can be applied to relate the solar zenith 
angle dependency to the observer's zenith angle dependency in the case of large surface reflectance (surface albedo « 
0.5 in Figure 3.4.6.3). 
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CHAPTER 4 
INVERSION ANALYSIS 

The inversion methodology for extracting optical parameters and eventually sulfate concentrations requires that the 
satellite measured total radiance be ccmvened to an apparent reflectance. This apparent reflectance along with all the 
operational and auxiallary environmental parameters is then entered into an iterative algorithm whose kemal is the 
5S model. In rcfq^nce to equation (2-10) and (3-1) we can describe the iterative operation as: 

t = f(p*(x).<t), e. eo,,X,<p>.tOoPo(cosx)) (4.1a) 

t,er*<=at = x(X, y, N(0), Haer, f, COo) - Tacr*''* - Xmol(X) (4-lb) 

504= = IS04=/Haer 

= 106 X ( t(K Y. f , (Oo) )-l Taer«=« / H«er (4-lc) 

where S04= is the equivalent sulfate surface density, Tacr^*» is the aerosol absorption optical depth and Xmol is the 
absorption and scattering optical depth for molecules. In equation (4-lb) the total optical depth (t) is shown 
explicitly to be a function of a vector of fractional number densities (y), the total surface number density (N(0)), the 
aerosol scale height Haer. the relative humidity (f) and the single scattering albedo ((Oo). The dependency of the 
aerosol absorption optical depth (which is the same as tile scattering optical deptii) has been suppressed. In order to 
investigate die sensitivity of sulfate extraction from the inversion procedure the differential dS04= wiU be 

investigated in terms of those parametCTS whose uncertainty can cause significant errors in the sulfate concentration 
calculation. Symbolically one can write: 

.SSO^. ,gSO:\ /8S0'\ 



6p 'P'-E ^ ge 



(4-2) 



'8©^ CP 



where die key parameters of surface albedo, apparent reflectance and sulfate mass scattering efficiency are show 
expliciUy while those variables considered as noise or of secondary importance are expressed in terms of the 
conEribution to the uncertainty in apparent reflectance (p*). As is the custom, the subscripts of each partial derivatii 
(hereafter termed error rates) indicate which parameters are held constant. This list of parametos is, for the sake < 
simpliciUy shortened to include only the three key parameters mentioned above. Whether one expresses the em 
directiy in toms of a given parameter of interest or as part of a contribution to Ap* is clearly subjective since tl 
latter cases can be equivalenUy expressed as: 



The three key error rates are evaluated in Section 4. 1 below while the separate error contributions (Ap*. Ap and Ae) 
are analyzed in Section 4.2. 

4.1 Error Rate Calculations ( SSO4-/SX ) 

Figure 4.1.1 is a typical calculation of sulfate equivalent surface concentration uncertainty per unit percent 
uncertainty in the apparent reflectance (5S04 / 5p*). These set of curves were computed for the urban variant model 
from the plots of apparent reflectance versus equivalent sulfate concentration shown in Figiue 3.4.2.1. The variation 
in Figure 4.1.1(a) shows clearly why the sulfate error is expected to be smaller at large equivalent sulfate 
concentrations and small ^ound surface reflectances (scattering atmosphere). Alternatively Figure 4.1.1(c) illustrates 
that the error per unit uncertainty in apparent reflectance decreases with decreasing surface concentration when the 
surface reflectance is large (attenuating atmosphere). This behavior clearly implies one technique for inversion qrp r 
minimization whic|i involves tlje choosing of tigght targets for tow sulfate loading and d aA targets for high sulfate 
liMiig. 

In the case of surface reflectances near the critical reflectance (Figure 4.1.1(b) ), one not unexpectedly, finds a large 
error potential about the sulfate values which display a minima in Figure 3.4.2.1(b). The error rate curves rapidly 
deteriorate to large negative positive values below and above the assymptope defmed by the minimum in Figure 

3.4.2.1(b). 

Figures 4.1.2 (a), (b) and (c) which were computed from the database used to create Figure 3.4.2.1 show the sulfate 
equivalent surface concentration deviation per unit increment in sulfate mass scattering efficiency ( (SSO4'' / 5e)p* 

modified by the factor • e). when the uncertainty in mass scattering efficiency is due to an uncertainty in relative 
humidity. These error rate values were effectively exacted by measuring the distance (ASO4) between relative 
humidity curves at constant apparent reflectance and dividing by the corresponding change in mass scattering 
efficiency. Their sign (in the absence of the multiplicative factor (- e) ) is always negative since an increase in 
relative humidity (which perforce implies an increase in scattering) must necessarily be accompanied by a decrease in 
sulfate concentration if the apparent reflectance is to remain constant 

The employment of the multiplicative factor ( •£) for the error rate was inspired by the results presented in Appendix 
A.3 where it was shown that the mass scattering efficiency error rate is to the first order givoi by - (S04~ / e). In 
the absence of other influences one expects the modified error rate to be equivalent to the sulfate concentration. Vix 
small surface reflectances (scattering dominated atmosphere) the acuial error rates lies below the y = x line due to 
second order effects induced by variations in the scattering phase function. If e increases by virtue of an increase in 
relative humidity the scattering phase function decreases (in the neighbourhood of large scattering angles typical of 
satellite remote sensing problems) because of the attendant decrease in re£ractive index and the increase in particle size 
(see Figioe 22,23 for exmnple). Because of this decease the required negative ofifset in sulfate concendaticm required 
to maintain p* ccMistant is reduced. This phenomenom becomes less apparent with increasing (qKical depth since 
larger order multiple scattering ejects tend to obscure the phase function influence. 

If the surface reflectance is lai^e (attenuating atmosphere) the phase function effects are less important at snail 
sulfate concentration and the curves all lie extremely close to the y = x line. At larger sulfate concenirations the 
scattering effects induce a greater scatter in the curves along with a deoease in the absolute magnitude of the slope. 
At surface reflectances close to the critical reflectance the transition betwe^ scattering and attenuating atmospheres 
cause rapid slope changes. 
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The variation between the y = x standard and the other error rate curves may be considered as a range of error rate 
values whose relevance depends on the inversion methodology and type of uncertainty influencing e. For example if 
the optical input parameters to the radiative transfer equation are known then in principle the inversion for aerosol 
optical depth (equation 4-lb) produces negligible error and the uncertainty in the sulfate concentration can be 
associated with the y s x line. Departures from this somewhat idealistic case include the relative humidity influence 
discussed above. 

The error rates (5S04= / 5p)p*^ computed for the same reflectance parameters as Figure 3.4.2.1 are illustrated in 
Figure 4.1.3 for two different relaUvc humidities. These values which represent the uncertainty in S04= per unit 
percent increment in surface albedo were computed by calculating the S04= increment (AS04=) between iso surface 
reflectance curves fcs- constant values of apparent reflectance. 

When the atmosphere is dominated by scattering effects (well below the criUcal reflectance) an increase in ground 
reflectance neccessitates an offsetting decrease in sulfate concentration (offsetting decrease in atmosphoic scattering) 
if the apparent reflectance is to remain constant The error rate in this case is accordingly less than zero and as weU 
decreases in absolute magnitude with increasing sulfate concentration as the atmosphere better masks the effects of 
the surface reflectance increase. 

Alternatively if the atmosphere is dominated by attenuation effects (well above the critical reflectance) an increase in 
surface reflectance necessarily requires an increase in sulfate concentration in order to offset (attenuate) the surface 
reflectance incease. In this case. then, the error rate is positive and increases with increasing sulfate concentration. 
The latter effect results because an increase in sulfate represents a gradual change from an attenuation dominated 
atmosphere to a scattering dominated atmosphere (influence of surface reflectance become increasingly less evident). 
Accordmgly, if one is at a point where attenuation still dominates but where scattering is significant (i.e. large 
sulfate concentration) the amount of sulfate which must be added to offset a reflectance increase is comparatively 
larger than that required at low sulfate concentration. Here again a mini mi/atinn of nnr^ntial error can in nrinrip IP h^ 
CfrcCtea t?Y an a ppropnate SflgCUOn of dark surfaces for large sulfate loadi ng and hriyht surfaces for low ^.ilfafP. 

loasliiiiL 

The passage from attenuation to scattering dominated atmosphere is marked by a rapid increase in positive error rate 
with increasing sulfate concentration if attenuation is dominate and a progressively more negative error rate with 
decreasmg sulfate if scattering is dominate. The two regions can in fact lie on eiihor side of an assymptotic sulfate 
value depending on the actual surface reflectance and the relative humidity (as is evident in some of the curves of 
Figure 4.1.3). 

44 Error contributions (AX) 

4.2.1 Effects of surface reflectance 

A critical requirement for any inversion scheme which attempts to extract an atmospheric parameter from the total 
measured satellite signal is some apriori knowledge of the contribution due to ground reflectance. In the past this 
cntena has been surmounted by restricting atmospheric inversion algorithms (or their complement, atmospheric 
correcuons) to clear water regions in the image where the known low reflectance of water contributes relatively Uttle 
to the total satellite signal (Ahem et a!.. 1977, Royer et al.. 1988b). This oftentimes severe constraint has led 
researchers to mvesugate the use of other terresmal surfaces for which the surface reflectance is small in certain 
wavelength regions (Kaufman and Sendra, 1988) or for which ground calibrations and/or sateUite temporal snidies 
can be used to extract a value of ground reflectance. Another alieni:,i ive is m vjpw the same nixe.l in iw^ 
QverlaPDiPg images acquired m two riiffprt^nt p a^^^. This approach assumes that the aerosol mass loading has not 
changed substantially between the acquisition time of the two images. The NOAA satellite sensors and in the future 
SEAWIFs are ideally suited to this purpose since the frequency of consecutive image acquisition is only of the ordo- 
of an hour. 
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Wiih no apriori information on the water qusdity parameters which characterize the optical turbidity of water one can 
assume a nominal value of reflectance associated with an uncwtainty characteristic of the natural variation of 
chlorophyll. Employing simulations such as those of O'Neill and Miller (1988) the magnitude of water leaving 
reflectance (ignoring the Fresnel reflection of skylight and assuming a scene target well removed &om the specular 
sun gutter direction) is of the order of .008 ± .005 in the blue (400 nm), .009 ± .003 in the neighbourhood of 500 
nm and .002 ± .001 near 675 nm. These ranges assume that the water can be categorized as being low in inorganic 
content (case 1 of Morel & Prieiff, 1977). 

Table 4.2.1.1 indicates the typical means and standard deviations which can be obtained over a variety of surfaces and 
over an extended time period if one anploys atmospherically and radiometrically corrected satellite data to compute 
surface reflectances. The procedure for extracting such data is distinct from the inversion procedure for obtaining 
optical depth and sulfate concentrations even though, in principle die same data could be used to do both. The entries 
of Table 4.2. 1 . 1 thus represent apriori information used to specify the values of the surface reflectance parameter to 
the radiative transfer calculations employed in the sulfate inversion procedure. The magnitude of the standard 
deviations reflect to some degree the effects of radiometric uncertainty but are probably more indicative of the 
refl«:tance precision which one can hope to achieve given the temporal and spatial variability of natural surfaces 
across a satellite scene. 

In tenns of amiOSOheric remote sensing the lower rcnftrtance.s of Tabl e 4.2.1 .1 such as those of water and forest 
vegetation in the absorbing spectral region of chlnmn h vll /MSS hand2> are tvpicallv ontimal fnr minimizing «tiilfat». 
JnYersion CrratL. in regions where the surface reflectance is large (Such as the reflectance peak of vegetation found in 
MSS bands 3 and 4) the standard deviation is correspondingly large. The inversion error effects of large uncertainti^ 
in the surface reflectance depend on the actual surface reflectance of the pixel being evaluated and the average 
reflectance of surrounding pixels which act as a secondary soiffce of reflected sunlight for aimospheric scattering. 

4.2.2 Sulfate Mass Scattering EITiciency 

The sulfate mass scattering cfflcicncy is influence! by a number of factors including relative humidity (0. effective 
particle radius at zero relative humidity (reff(0)) and the relative contribution of sulfates to total aerosol scattering. 
Symbolically we may write the total error as: 

where y represents the fractional number density due lo sulfates. The use of this latter parameter is somewhat of a 
simplification since the variation of e due to variations in the contributions of constituent aerosols is clearly a non 
trivial phenomenom. 

Depending on the apriori information available the above equation can represents a complex variety of sulfate 
inversion schemes. If one computes the scattering efficiency from first principles as was done in our simulations 
then the uncertainties in each component represent measurement uncertainties. On the other hand if ground based 
measurements of e(f jeff(0).Y) were available in a given scene then the uncmainties are jBsociated with the natural 
variation of each parameter across the satellite image on which the sulfate inversion is being performed. 



Table 4.2.1.1 Table of measured ground reflectances and their uncertainties (in 

parenthesis) for a number of surface types derived from Landsat MSS data 
(from Royer et al. 1988) 
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LAND 

COVER 

TYPES 


MSSl 


MSS2 


MSS3 


MSS4 

MSS 

ALBEDO 


MEAN 


Agriculture 


0.067 
(0.016) 


0.065 
(0.021) 


0.231 
(0.036) 


0.289 
(0.071) 


0.175 
(0.030) 


Forest 


0.042 
(0.010) 


0.031 
(0.007) 


0.197 
(0.027) 


0.306 
(0.054) 


0.163 
(0.017) 


Urban 

(high density) 


0.092 
(0.021) 


0.092 
(0.019) 


0.112 
(0.018) 


0.117 

(0.019) 


0.105 
(0.018) 


Gravel Pits 


0.127 
(0.026) 


0.124 
(0.022) 


0.136 
(0.020) 


0.137 
(0.001) 


0.130 
(0.018) 


Water 


0.010 
(0.003) 


0.005 
(0.002) 


0.002 
(0.001) 


0.001 
(0.001) 


0.005 
(0.001) 
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(a) Variations in Relative Humidity 

The results of both chapters 2 and 3 indicate an evident requirement for priori infonnaiion on relative humidity if 
the inversion methodology depends on a relative humidity model for the calculation of the sulfate mass scattering 
efficiency. Ideally this information should be extracted from strategically located spectral bands on the same sensor or 
sensor platfonn which attempts lo measure the sulfate induced signal. 

The most promising satellite based method developed to date has been the split window technique which employs a 
differencing algorithm between two or more infrared bands to extract an integrated water content measure (Deschamps 
et al., 1980). Dalu (1986) reported retrieval accuracies of ± 3 and ± 5 kg/m^ for the IRIS (Nimbus 7) spectrometer 
and AVHRR respectively (.3 and .5 cm equivalent water column). 

The form of the water vapour profile with altitude is sufficiently stable that a soong correlation exists between 
integrated water vapour content and relative humidity at the surface (Hay, 1971). One can acconlingly transform 
satellite derived water content measures into surface relative humidity which in turn can be directly related to the 
explicitly required parameters employed in the inversion algorithm. For the great lakes region Hay's relation 
becomes: 

WS.15 + .0013 fes (4-5) 

where w is the integrated water vapour content in equivalent cm. f is the relative humidity in percent and Cs is the 

saturation vapour pressure. For a surface temperature of IS^C the eiror in the computed value of the surface relative 
humidity (Af) corresponding to the satellite retrieval error (Aw) is simply: 

Af= lOOAw/2.39 (4^ 

or the order of 10% and 25% for the IRIS and AVHRR sensors respectively. 

(b) variations in dry sulfate log normal parameters 

In Appendix A.3 it is demcmstraied that the aerosol scattering efficiency per unit mass of dry sulfate does not vary 
greatly over a substantial range of the dry log iK»rmal geometric mean radius and the log normal geometric standard 
deviation (iN(0) and a defined in Appendix A.2 and Table 2.1.1.1). 

This phenomenon which results because aerosol scattering is predominantly volume or mass dependent in the fine 
particle (sulfate) range implies that the oror analysis arguments developed above for fixed iN(0) and o and referenced 
to dry sulfate mass are applicable to log normal parameter values larger than the nominal values given in Table 
2.1.1.1 (.12 and 124 for rN(0) and a respectively). Accordinylv ffor identical a mbient conditions! a unit vnliime nf 
la rge Slllfale panicles produces aPDroximatelv equivalent ontical/ra d iarive tran-;fPT effecLs as a unit volume nf ■small 
but more numerous sulfates if their mass is identical. This assumes that changes in directional scattering properties 
(phase function) are second order widi respect to their effect on the exoatmspheric reflectance (refer to the discussion 
in Appendix A.3). The range of particle size over which the approximate constancy of e(X) holds is a function of 
relative humidity and wavelength. However in coarse terms one may vary a between 1.2 and 1.8 for rN(0) = .12 and 
between 1.2 and 1.4 for rN(0) = J24 without affa:ting a change of greater than 20 % between any of the computed 
values (relative humidity < 90%, wavelength = 0.55 \im). 

Outside of this region or for more precise calculations one may employ all the error rate versus sulfate cwicentration 
curves reported above by utilizing an approximate normalizing factor. The method of computing this normalizing 
factor is demonstrated in Appendix A.3. Generally a decrease in mass scattering efficiency implies a proportionate 
increase in the magnitude of the computed error rate. 
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(c) Contributions due to otlier fine particle mode constituents 

The discussion presented in Section 2.2.3 indicated that the most likely reason for the rather large discrepency of c 
values derived from ground truth measurements was due to the covariance of other fine particle constituents with 
sulfate. Comparisons of measurements made by different authors show a strong discrepency which would a^jpear to 
be due to the local nature of the fine particle aerosol mode. For the most part however the measurements made by 
any given author almost invariably indicate a strong aerosol scattering coefficient to sulfate mass correlation and 
hence a correspondingly small uncertainty associated with the value of e (typically ~ 1 m^/g). 

4.2.3 Influence of. Optical Absorption (Single Scattering Albedo) 

Kaufman (1987) developed a method for measuring the average aerosol atmospheric single scattering albedo ((Oo) 
based on a multi image determination of the critical surface reflectance (c.f. Sections 3.4.2 and 3.4.3). His technique 
exploits the presence of optical depth singularities such as those of Figure 3.4.2.2 by employing the corresponding 
critical reflectance as an indicator of the single scattering albedo. Since the critical surface reflectance is independent 
of optical depth it's variation is strongly correlated with atmospheric absorption and accordingly facilitates the 
development of an inversion methodology for the extracung of (Oo- The weakness of the technique resides in its 
dependence on an invariant values of coo and surface albedo over two or more images acquired at different dates. 
Given two images for which (Oo is approximately constant Kaufman reports a computed error of less than 10% for 
reasonably typical environmental and operational parameters (solar zenith angle = 45 degrees, refractive index = 1.43 
-.00351, log linear size distribution, wavelength ^ 610 nm) and associated errors. 

By the same token localized ground measurements of carbonaceous soot absorption on the filter of a high volume 
sampler (such as those used by the Ontario Minisicy of the Environment for their pollution index) may be employed 
to compute a ground level value of absorption coefficient (Rosen and Novakov, 1983) from which a value of oJo 
may be extracted if scattering or extinction coefficient measurement are available. Alternatively, if estimates of the 
relative humidity and the number density (or mass) ratio of soot particles to sulfates are available a value of the 
aerosol single scattering albedo can be extracted (c.f. A^wndix A.3 (b) ). The tacit assumption of this approach, that 
externally mixed (independent) carbonaceous soot panicles play a dominant role in atmospheric absorption, is 
strongly supported in the literature (Bergstrom, 1973; Rosen et al., 1978; Japar et al., 1986) . 

The accuracy of the absorbing and scattering coefficient measurements which one can infer from the literature to be 
of the order of 10-20 % (Japar et al., 1986; Ruby and Waggoner, 1981) implies a similar uncertainty in ©o. A 
vertical profile of tOo must then be assumed or measured in order to establish an average atmospheric value (the 
simplist approach is to assume that the mixing ratio of aerosol scattering to absorbing particles is effectively 
independent of aldtude). 

If only measures of soot absorption are available one can stiU inf« an average atmospheric value of coo given some 
estimate of the scale height of the soot panicle along with a ground based measure of total optical depth. In the 
absence of any absorption information whatsoever one can set coarse limits of .8 ± .2 in an urban atmosphere and .9 
± .1 in a rural atmosphere for visible wavelengths of radiation (Ackerman and Toon, 1981, Table 4.2.3.1). 

The computations of Kaufman (1987) indicate apparent reflectance uncertaindes which range from .002- .005 per .01 
in single scattering albedo increment at large values of surface reflectance (.4) down to values less than .001 per .01 
increment in a>o at zero surface reflectance. These calculations were performed at the solar zenith angle, refractive 
index and wavelength parameters given above over a aerosol scattering optical depth variadcm between and .6. This 
variation in aerosol scattering optical depth roughly corresponds to our standard - 50 |ig/m3 equivalent sulfate 
surface concentration variation at 90 % relative humidity at an aerosol scale height of 1 km. In comparison the 
difference between the continental standard and urban standard results for which the single scattering albedo goes 
through a significant change in Wo (approximately .164 at 90% relative humidity and 0.633 |im wavelength) 
indicated uncertainty magnitudes of less than .0005 and .004 respectively for and .4 surface reflectance across an 
coo increment of .01. 
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Table 4.2.3.1 

Single scattering albedo as a function of wavelength, relative humidltu and aerosol model 



Relative humiditu 



OS 



0,436 
0,728 
0,965 

0,725 
0.898 



50 % 



0,436 
0,726 
0,974 

0.752 
0,914 



70 % 



0,436 
0,728 
0,984 

0,794 
0,936 



60 % 



0.436 
0,728 
0.968 

0,819 
0.946 



90 X 



0,436 
0,728 
0,994 

0,87 1 
0,966 



95 « 



0,436 
0.728 
0.998 

0.921 
0,981 



0.55 micrometers 

Component narnsnls 

soot 

dust like 

water sol. 

Optical averages 

urtoen/lndustrial standard 

continental standard 

3.75 micrometers 

Component aerosols 
soot 

dust like 
I water sol. 
Cptlcal averages 
urtan/lndustnal standard 
continental standard 



scat effec = change In the total volume scat coet per unit change In the mass of dry sulfate 
cross sections obtained by multiplying Mie effeclencles (Q) by p1*r2*«2 where r2 Is the effective 

mean surface radius of a log normal distribution (see below) 
Optical average computed as the weighted mean of scattering cross sections normalized to the weighted mean of 
total cross sections (Appendix B). The weights are computed as fractional number densities (Table 4) 



Dust like Soot 



0,015 


0,015 


0,015 


0,0)5 


0,015 


0,015 


0,880 


0,880 


0,880 


0,880 


0,880 


0,880 


0,258 


0,289 


0,346 


0,363 


0,475 


0,594 


0,193 


0,194 


0,199 


0,204 


0,224 


0,283 


0,834 


0,833 


0,829 


0,826 


0,819 


0,809 



Model 

r2 (microm.) 

Cross sections (sIgma) [units of I0«*-10 cm*»2l 



Water soluble (sulfate) 
■ OX • 50f 70f BOX 90X 95X 

1,096704 0,125683 0,125683 0,134909 0,150531 0,160488 0,185800 0,222392 



0.55 micrometers 

Sigma (scat) 625,871 5,22423 5,0)703 5.9671 8,06198 9,72656 15,3164 27,542 

Sigma (abs) 233.511 6.7545 0.(6244 0.15912 0.12633 0,H3436 0,085325 0,060371 

3.75 mlcrometsrs 

Sigma (scat) 622,9660 0.014128 0,003305 0,004752 0,008482 0,012009 0,027173 0,076048 

Sigma (abs) 85,34760 0.944820 0,009523 0,011668 0,016047 0,019380 0,030049 0,052040 



m 



Table 4.2.3.2 below summarizes calculations of the error rate (6p* / 5o)o)rJi..p performed for nadir geometry and a 

wavelength of .55 \un. The calculations were performed using the data points plotted in Figure 3.4.1.3 (along with a 
sa»nd set for 90% relative humidity). The significantly greater values of the error rate at the larger value of surface 
reflectance (0.4) illustrate weU the inherent diffo^ice between the effects of soot particles at small and targe surface 
reflectances (c.f. Sections 3.4.1 and 3.4.4). 



Table 4.2.3.2 Error rate (5p* / 5fflo)r.h.,p for nadir geometry, 45 degrees solar 

zenith angle and 0.55 ^.m wavelength (units of percentage change in p* per .01 
increment in (Oo) 

Surf. ref. 0.5 



Rel. hum. 70 90 W W 

Sulfate concentration 

0. 



OJ 


m 


0,0 


m 


-.011 


-.028 


0.20 


0.34 


-.tm 


-.039 


0.38 


0.63 



4.2.4 Noise and Radiometric resolution 

Table 4.2.4.1 is a compilation of radiometric reflectance errors as a function of wavelength band for the principal 
satellite sensors being investigated. The noise equivalent reflectances (NER) represent actual noise speciflcations for 
the given sensor which are eithor evaluated in pre launch calibrations or during the operational lifetime through the 
utilization of standard sources and/or homogeneous ground targets. The quantitization etnx reflectance (QER) is 
simply the quantitization resolution of the sensors translated into an apparent reflectance error. 

Since the satellite sensor specifications are typically given in terms of radiance or digital levds one employs the 
apparent reflectance relation: 

p* = jtL/nsEt (4-7) 

to write: 

Ap* s nAL/M.tEf (4-8) 

where AL is an error specification in terms of radiance, \it is the cosine of the solar zenith angle, and Eg is the 
exoatmospheric solar irradiance. The i)arameter \is was wherever possible taken to be the nominal value of 0.5 (soiar 
zenith angle of 60 degrees) in order to yield a maximum noise ovo' a realistic range of solar zenith angles. The solar 
irradiance was eithor extracted from the specifications for the sensor of interest or when this was not available 
computed from siandju-d tables (Wolfe & Zissis, 1978) assuming a square spectral response function. 

Often the noise figure for a particular sensor is given in terms of signal to noise ratio (SNR) at a particular 
reflectance (p*). In this case the NER can be obtained by ratioing the above two equations: 



TflbU 4.2.4.1 Nolst Equlvslant Riflictanct* (NER) and 
QuantUlzallon Errcr RafUclanctt (QER) 

SaUIIHt Sansor Band Bendwtdth 

[micrometers! 



NER 



bus/ 

pixel 



QER 
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Rer. Total error Notes 



Land»Bt -5 



Tn 



1 

2 
3 

4 
5 

7 



.45 - .52 
.52 - .60 
.63 - .69 
.76 - .90 
1.55- 1.75 
2.08 - 2.35 



0,0065 
0,0057 
0,0057 
0,0033 
0,0168 
0,02 



0,0011 
0,0023 
0,0019 
0,0020 
0,0018 
0,0028 



1,8 



0,0066 ( 

0,0062 

0,0058 

0,0043 

0,0169 

0,0202 



SI»flT 



nss 



mv: 



l 

2 
3 
4 

XS1 
XS2 
XS3 



.5 -.6 
.6 -.7 
.7- .8 
.8- 1.1 

.50 - .59 

.61 - .69 
.79 - .90 



0,0057 
0,0057 
0,0065 
0,007 

0,0036 
0,0056 
0,0031 



0,0036 
0,0029 
0,0033 
0,0041 

0,0021 
0,0025 
0,0035 



1,0 



2,9 



0,0068 
0.0064 
0,0073 
0,0081 

0,0044 • •• 

0.0061 

0,0047 



NOAA AVHRR 1 .58 - .68 0,0017 10 0,00058 3 

(9,10,11) 2 .715 -.982 0,0017 10 0,00058 



0.0018 •,! 
0,0018 



heteosat 

QOES VAS/VtSSER 



.4- 1.1 
.54- 75 



0,0013 

NA 



0,0023 4,7 
0,0025 4,6 



0,0027 • 
0,0018 •,•• 



LandS9i-6 



SEAWIFS 



I 

2 
3 

4 
5 
6 



.443 - 453 
.490 - .510 
.555 - .575 
.655 - .675 
.745 - .785 
.843 - ,887 



0,00049 
0.00035 
0,00031 
0,00023 
0,00019 
0,00016 



10 
10 
10 
10 

to 

10 



0,00023 
0,00017 
0,00016 
0,000092 
0,000081 
0,000069 



0,00054 
0,00039 
0,00035 
0,00025 
0,00039 
0,00017 
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Notes 

♦ calculated for a maximum digital level of 64 (radiance or nflcctance varies 

as the square of the digital level 
" NER r quadrature sum of 60 deg. sol. zen. angle along column and cross 
diode noise In ref. 2 

♦ NER s reflectance / SNR (reflectance either quoted or computed from nominal 

radiance and E(0) where 1(0) = the exoatmospherlc solar irradlance 
$ nominal values, slight differences exist between each NOAA tensor 
QER s reflectance error given an error of one quontttlzatlon Itvel s pl'Iradlance p 

digital level! / lcos(SZA) • E(0) • SQRT(3) 1 
C Total error * SQRT(N£R-2 ♦ QER"2) 
** solar zenith angle z 60 degrees 

♦ Launch scheduled for 1990, values shown ere expected specifications 



m 



Ap* = p*/SNR (4-9) 

The radiometric resolution uncertainty (QER) represents a band of values encompassed by ± I about any digital 
level. Accordingly, if the sensor is characterized by "n" bits per pixel, and radiances Lmin and Lmw cwresponding to 
zero and the maximum number of digital levels respectively, the quantitization error radiance is given by: 

AL2 = [ (Lm« - Ltnin) (l/2n)]2 / 3 . (4-10) 

where the divisor 3 results from a computation of the root mean square value of a pillbox probability encompassed 
by the digital level limits of -1 and 1. After conversion to reflectance the two contributions to the degradation of 
radiometric resolution are combined according to: 

Ap* = [ Ap*(NER)2 + Ap*(QER)2 ] 1/2 (4-11) 

where one assumes that the two sources of noise are independent. 

The NER and QER values of Table 4.2.4.1 computed for SEAWIFS were deduced from the nominal values detailed 
in Table A-2 of Baker et al. (1987). Lmax in this case was assum«i to be the saturation radiances while Lmin was 
taken to be approximately zero. It is worth noting that the SEAWIFS target radiometric parameters are at least a 
factor of two better than any of die other reflectance error values given in Table 4.2.4.1. In general the current 
generation of 7 and 8 bit sensors are characterized by a NER of the order of half a percent while die 10 bit NOAA 
sensors are down in the .2 % range. In all cases the QER values are commensurate widi NER tabulations. 

4.2.5 Aggregate aerosol model (discrimination of component aerosol contributions) 

One of the most critical elements in die development of an inversiOT methodology specific to sulfates is clearly how 
to distinguish the radiadve transfer effects of sulfate from those of other aerosol contributors. In principle this can be 
achieved given a sufficient number of spectral bands and a disdnct difference between the spectral behavior of the 
major opucal contributors. 

The most easily accomplished discrimination is effected between the dust like and sulfate components. Figure 4.2.5 
shows some calculations of apparent reflectance as a function of wavelength for two different values of total aerosol 
content (defined in terms of the total aerosol optical depth at O.SS M.m), a relative humidity of S0% and at zero 
surface reflectance. In order to better isolate die reladve influence of dust and sulfate panicles the separate curves have 
been calculated for pure samples of each type of panicle (in the nomenclature of Table 2.1.3. 1 two externally mixed 
aggregate aerosol models were defined widi (i) Yi(S04") = 1. Tf2(dust) « and (ii) Yi » 1, Tf2 = 0). As well die pure 
Rayleigh contribudon (zero aerosol opdcal depth) has been subtracted in order to enhance the contrast between the 
two curves. 

As die wavelength increases from the visible spectral region to the near IR (.4 - 3.75 pm) one can characterize a 
transition of sulfate particles from optically small (1 < 2jcreff / A. < 10) to optically very small (2jcreff / X < 1) and 
dust like particles from optically large (2nreff / X > 10) to opdcally ^nall. The optical effects associated with these 
transitiOTS are such diat the sulfate volume scattering coefficient become increasingly less important widi increasing 
wavelengdi while that of die dust like particles remains relatively unaffected In the range 0.4 to 0.86 jam the sulfate 
extinction which is approximately die same as the scattering optical depdi decays roughly as X-2-4 whDe at longer 
wavelengths the extinction decay is even more rapid with Rayleigh absorption playing a progressively more 
important role (see Section 2.2.2 for a discussion of Rayleigh absorption). 
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The net effect of these spectrally different optical properties is to produce a dominance of sulfate sclinering in the 
visible and a dominance of dust scattering in the near IR. The crossover point docs not occur at 0.55 nm (i.e. where 
optical depths are equal) but rather at a larger wavelength due to the fact that the optically large dust like particles 
suffer from weak backscattering properties (comparatively small phase function values outside of the anisotropic 
cones at and 180 degrees scattering angle). A sudden drop in the dust like backscatter phase function for the 225 
\xm and 3.75 jim wavelengths produces the comparatively low reflectance values seen in Figure 4.2.5. The phase 
function decrease results from a substantial decrease in the real part of the refractive index (see Table 22.1.2 for 
example). 

Spectral data can accordinylv be employed to distinguish dust like for even cloud> particles from sulfate naTtirlft<! 
without any apriori information on the content of the aeros ol cloud in the field of view. To avoid cumbersome 
multiparameter inversion schemes the problem of properly interpreting spectral data should be approached from a 
semi-empirical correlation search between the spectral curve slopes observed in Figure 4.2.5 and the contribution of 
the major aerosol optical components. If this procedure is found to be reasonably independent of optical depth 
(aerosol mass loading) one can develop an inversion methodology for extracting particle size information and sul£ate 
concentration extraction in two independent steps. 

The concentration of carbonaceous soot relative to the concentration of sulfates can be estimated if an independent 
measure of the aerosol single scattering albedo can be computed from the satellite data (c.f. Appendix A3). This 
could be computed using Kaufman's technique although the requirements of this approach appear to be rather 
stringent (c.f. Section 4.2.3). An alternative £q>proach would be to exploit the comparatively high sensitivity of the 
apparent reflectance to soot concentration over high reflectance surfaces (c.f. Section 3.4.1). In this way hig hl y 
reflecting pixels in a given scene would be the main source o f soot concentration information while low reflectance 
pixels would be used as the primary source of information for sulfate c oncentration. The employment of highly 
reflecting surfaces necessarily implies a calibration program to deduce the intrinsic reflectance of the targeted surface 
or the utilization of multi look angles afforded by overlapping images. 

It is difficult to imagine a satellite based inversion algorithm, even a multispectral inversion algorithm capable of 
detecting realistic concentrations of secondary scattering particles (such as ammonium nitrate) based on broadband 
responses to atmospheric scattering. This lack of sensitivity results because of the lack of difference between the 
intrinsic optical properties of the secondary scattoers and because scattering is largely a volum^iic effect independent 
of the chemical nature of the scattering particle (see Appendix A.3). It is. on tl)e other hand, entirely feasibly ^q 
envisage a high spectral resolution imajging svstem fsuch as the EOS MQDIS a nd HTRIS sensors^ which deterts 
anthropQgenicallv generated gases bv differential absofptio n techniqy ies , In this case the strategic placement of high 
resolution spectral bands about one or more absorption features unique to a given molecular species provides a 
highly discriminating means of trace gas detection. Ground based instruments of this nature have been developed for 
the detection of ozone, nitrogen dioxide and sulfur dioxide (Kerr et al., 1980; Syed &. Harrison, 1980; Kerr & Evans, 
1980). 

4.2.6 Radiative transfer model 

(a) analytical approximations 

The results presented in Appendix B (see also Section 3.2) indicate that at near nadir directions the error expected in 
employing the 5S analytic^ approximate is of the order of .005 in apparent reflectance for solar zenith angles less 
than 60 degrees, wavelengths between 0.4 pm and 1 .6 |im, stirface albedos between and 0.4 and an aerosol optical 
depth in the neighbourhood of 0.5. Assuming a reasonably linear correlation between apparent reflectance and optical 
depth this corresponds to a maximum error of the order of .01 at an aerosol optical depth of unity or hence in 
conditions roughly equivalent to a wet atmosphere of 95% relative humidity with an equivalent sulfate surface 
density of 50 ng/m3 (c.f. Table 2.1.4.3). 

In off nadir directions the apparent reflectance error derived in Appendix B increases rapidly with the cosine of the 
observers zenith angle to values as large as .05 at 60 degrees (solar zenith angle less than 45 degrees, wavelength of 
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.55 (xm) or hence roughly .1 for the extreme case discussed immediately above. This error is more a function of the 
analytical approximations employed in the 5S model than of the problems associated with vertical inhomogeneity in 
the atmosphere (c.f. Section 3.2 above and (b) below). 

(b) vertical inhomogeniety 

Comparisons of the satellite level signal computed for a multilayer inhomogeneous atmosphere and a single layer 
model (c.f. Section 3.3 and Appendix B) indicated that the apparent reflectance error associated with replacing the 
vertically detailed atmosphere by an equivalent average atmosphere were small. At near nadir directions the error was 
of the same order as 5S while for off nadir directions the error was significantly smailo- than that of 5S (less than 
.005 for solar zenith angles up to 60 degrees and a wavelength of 0.55 jim). 

4.2.7 Other Sources of Error 

Sections 4.2.1 to 4.2.6 deal with what we consider to be the major sources of error in a sulfate inversion 
methodology for which different levels of apriori information can be exploited. Other sources of error which are 
typically second order in magnitude exist and can contribute to the error budget. These include polarization 
contributions from the atmosphere and the sensor optics, uncertainty in Rayleigh (molecular) scattering 
parameterization and variations associated with changes in surface temperature and pressure (Penndorf, 1957; Young, 
1981), variations in ozone optical depth (O'Neill, 1982). uncertainties associated with the parameterization of water 
vapour absorption in the visible continuum (Tomasi, 1979) and nitrogen dioxide in the blue spectral region (O'Neill. 
1982) and non Lambertian surface reflectance effects. 

4.3 Cumulative sulfate Errors (AS04=) 

Table 4.3.1 lists typical input parameter error values required in order to achieve a precision of + 1.6 |ig/m3 per 
input parameter in the extraction of sulfate concentration. The choice of 1.6 Hg/m3 is simply a figure of merit which 
is indicative of the average precision which should be attained if a total precision of ± 5 Hg/m3 is targeted in the 
presence of 10 independant and equal errors ( [10*1.62](l/2) _ 5 |ig/m3 ). 

The three sets of entries defined under the heading scattering efficiencies refer to the precision required in the relative 
humidity, the effective radius and the sulfate fractional number density (as per Section 4.2.2) for a sulfate precision 
per component of 1.6 \ig/m^. The manner in which one interprets this required precision is again a function of the 
apriori information available for determining the mass scattering efiiciency. As stated in Section 4.2.2 two pertinent 
examples include (a) the scenario where the mass scattering efficiency is calculated from first principles and the 
dependent parameter uncertainties (in f, reff(0) and y) are measurement errors or (b) measurements of e(X) are 
available in the image and the dependent param^cr uncertainties represent neural variations between the ground 
points. In any case, as shown in Appendix A.3, the mass scattering efficiency error rate fit>m which the prticision 
values were determined is always of the order of - S04= / e(X). 

At larger relative humidities (e.g. 90 % in Table 4.3.1) the required precision in relative humidity ( ~ 1%) is 
unattainable by current methods applied to existing satellite sensors (see Section 4.2.2 above). One must eithor relax 
the precision requirements on this variable or appeal to other measuring techniques. Ground based measurements 
certainly achieve this order of accuracy and could be employed eithor to calibrate existing satellite estimates of water 
content or directly &x>m a surface mapping tied to the measured values. The a;curacy periformance. in particular as 
related to the spatial surface gradient of relative humidity was not investigated. 

In die near future the space station sensor MODIS (Esaias. 1986) will incorporate spectral bands specifically for 
measurements of water content by means of differential absorption. In addition the instrument package will include 
an active laser system (LASA) whose capabilities include the measurement of water vapour profiles by means of 
differential absorption. 
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concentration precision of 1.6 ug/nv**3 per parameter 
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Atmospheric 
/geometric 
input parms. 

solar zen. ang. (deg.) 
obser. zen. ang. (deg.) 
obser. azi. ang. (deg.) 
surface ref. 
rel. humidity 
wavelength (um) 
sulfate cone (ug/m»*3) 
scale height (km) 



Urban variant model 
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45 


45 


45 


45 


45 


45 


45 






























































.5 


.5 


.5 


.5 


70 


90 


70 


90 


70 


90 


70 


90 


.55 


.55 


.55 


.55 


.55 


.55 


.55 


.55 


25 


25 


50 


50 


25 


25 


50 


50 


1 


1 


1 


1 


1 


1 


1 


1 



required Input 
parameter precision 
scattering efficiency 

(a) rel. humidity 

(b) effective red. (um) 
(c)S04fr8c. • density 
surface ref. 

sing. scat, albedo 
noise reflectance 
S04 froc. * density 
scale height (km) 



3.6 


0.5 


1.9 


.22 


4.6 


.40 


1.8 


.17 


.06 


.04 


.03 


.02 


.05 


.05 


.03 


.03 


0.2 


0.2 


.08 


0.1 


0.1 


0.3 


.08 


0.2 


.001 


.002 


.001 


.003 


.001 


.001 


.001 


.0005 


.06 


.05 


.04 


.04 


.005 


.003 


.002 


.0007 


0007 


.001 


.0008 


.002 


.001 


.001 


.0006 


.0004 


.07 


.06 


.06 


.05 


n/a 


n/a 


n/a 


n/a 


.06 


.06 


.03 


.03 


.06 


.06 


.03 


.03 



Notes on parameter 
precision 
scattering efficiency 

(a) rel. humidity 

(b) effective rad. 

(c) S04 frac. • density 



surface ref. 
sing. scat, albedo 
noise reflectance 
sulfate frac. • density 



scale height (km) 



computed from Figure 4.1.2 * d(eps11on)/df computed fi 
computed from Figure 4.1.2 and Figure A.3.1 
computed from Figure 4.1.2 * d(eps11on)/d(gamma) com 
Table 2.1.4.3 (urban standard minus cont. standard) anc 
epsllon values (needed to denormallze Figure 4.1.2) froi 
computed from Figure 4.1.3 

computed from Figures 3.4.1.1,3.4.1.2, 4.1.1, and Tal 
computed from Figure 4.1.1 

computed using Figure 4.2.5 (difference between pure si 
dust like converted to an approximate value for 
d(app. ref. )/d(su1 fate frac.) and Figure 4.1.1) 
using equation (4-lc) 



r- "^"1 
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The precision values listed for the dry effective radius were obtained by computing a value of 6e/Sreff(0) from the 
weak slopes observable in Figiu-e A.3.1 (approximately 1.1 and 2.2 per .1 ^m increment in effective radius for 70 
and 90 % relative humidity respectively). The computed values are within the naturally occurring variation about the 
mean reported by Young et al. (1988) and Lcaitch and Isaac (1988) (note that otc must COTvert their geometric mean 
radii to values of effective radii using the relation given in Appendix A.3). 

The third source of error in the mass scattering efficiency relates to the contribution of other fine particle constituents 
to the total aerosol scattering. The values shown in the table were computed from the differences in the urban and 
continental aggregate aerosol models (Table 2.1.4.3). As pointed out in Section 4.22 (c), this source of error alone 
induces variations in the mass scattering efficiency which are significant. The implication of these findings is that a 
working value of this t^arameter must be measure d or inferred on a regional scale. 

For the conditions defined in Table 4.3.1 ground surface reflectances need to be known to an accuracy of about .001. 
The surface reflectance uncertainties tabulated in Table 4.2.1.1 indicate that this level of accuracy is only achievable 
over water surfaces in the visible wavelength region. The use of other surfaces is not excluded but would require a 
dedicated calibration program to reduce the typical oror levels associated with remotely sensed surface reflectance. 
The most pffomising surface other than water is vegetation in the chlorophyll absorption wavelength band centered at 
about .65 \m. Table 4.3.1 excludes anv calculations of errors near the critical surface reflectance since the results of 
Section 4. l demonstrate clearlv that the associat ed inversion errors would be prohibitive, 

For small surface reflectances the addition of soot particles to a relatively pure sulfate mixture has minimal effect on 
the computed apparent reflectance (c.f. Sections 3.4.1 and 4.2.3). Accordingly the required single scattering albedo 
accuracy in Table 4.3.1 is achievable by a few well placed ground measurements (for example). For large surface 
reflectances, howevCT, the apparent reflectance is strongly dependent on changes in soot concentration and the 
required accuracy is prohibitive (better than .01). This dependency can however be turned lo advantage by developing 
reflectance inversion techniques for extracting soot indicators over highly reflecting surfaces while emphasizing 
sulfate inversion over low reflecting surfaces. 

In the absence of a single scattering albedo estimate an inversion applied as per equation 4-1 yields a scattering 
optical depth which is overestimated but which, at least for large surface reflectances, will be an a^woximation to 
the total optical depth. Accordingly the nominal sulfate mass computed using equation 4-lc is actually a value 
related to the combined concentration of soot and sulfates. 

The noise and radiOfnetric resolution requirements of Table 4 3 1 are nnlv achieved hv the 10 hit NOAA an^ 
SEAWIFs sensors and the 6 bit METEOSAT and GOES isensnr s in terms of the snecificatinn.s nresent<^d in Table 
4.2.^.1. The rather coarse radiometric resolution of the GOES and METEOSAT sensors is in fact well suited to 
sulfate detection levels since the range of these sensors is adapted to the low light signals associated with 
atmospheric phenomena. These sensors are however constrained by the low saturadon levels which necessarily result 
from their smaU dynamic range. 

The accuracy required in the sulfate fractional number density is derived from Figure 4.2.5 and indicates the 
uncertainty permitted in the sulfate fractional number density (t^ defined in Appendix A.1) to achieve the 1.6 (ig/m^ 
figure of merit. The procedure for obtaining these entries involved a coarse approximation of converting the optical 
depths in Figure 4.2.5 (.2 and .4) to equivalent sulfate concentrations at 70 % and 90 % relative humidity (using 
Table 2.1.4.3). and then, assuming bnearity of apparent reflectance and sulfate concentration, extrapolating the Ap* 
obtained from Figure 4.2.5 (at 0.55iim) to 25 and 50 |ig/m3. The figure quoted in Table 4.3.1 is accordingly the 
allowable error in the sulfate fractional number density between an aerosol atinospheric slate of pure sulfates (ys = 1) 
and an aerosol aunosphere characterized by an equivalent optical load of pure dust particles (7i = 0). The computed 
values for 0ys/^p*kfi are shown in Table 4.3.2 below: 
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Table 4.3-2 (Sys / 5p*)E,p values (per percent in p*) 

rel. hum 70% 90% 

sulfate cone. 
[|ig/m3] 
25 Id M 

m M M 

An estimate of Ys for (at least for the fine particle fraction for which sulfates are the dominant component) is 
possible in the light of the discussion presented in Section 4.2.5(b) above. The implementation of such an approach 
requires two or more sateUite bands which span the range between .4 and 1.6 jim and a noise/radiomctric resolution 
superior to the figures presented in Table 4.2.4,1. 

The figure of merit for sulfate mass scattering efficiency is simply derived from equation 4-lc where one assumes an 
accurate scattering optical depth has been extracted and hence where the ratio Ae/e (expressed as a percent) is simply 
equivalent to the relative uncertainty in equivalent sulfate concentration. The mqjor contributor to Ae is the 
uncertainty asscx:iated with the relative contribution of sulfa tes to scattering bv the fine particle mode. Natural 
variations of e (c.f. Table 2.1.4.1 and Section 2.23) considerably exceed the required precision of Table 4.3.1. The 
most pragmatic method of reducing this uncaiainty is to correlate satellite derived values of aerosol scattering optical 
depth with ground based measurements of dry sulfate concentration over measuring sites where the other input 
optical parameters have been evaluated. The verification of this approach would best be accomplished by 
simultaneous nephclometer measurements. Altonatively a season of nq)helometer measurements at rural and urban 
sites could be employed to establish empirical values of the mass scattering efficiency as a function of a minimal set 
of meteorological parameters. 

The scale height entry indicates the level of accuracy required in this parameter if one attempts to extract a real value 
of sulfate surface concentration. As above this error was derived from equation 4-lc assuming no error in the other 
dependent parameters. 

4.4 Operational Considerations 

4.4.1 Geometrical Effects 

(a) Observation Geometry 

As Stated in Section 3.4.6 the principal effect of non nadir geometry in a zero reflectance atmosphere is to induce 
apparent reflectances which vary roughly as 1 / ji (where \l » cosB is the cosine of the observers zenith angle) for 
small surface albedo. In keeping with the ideas prcsentwi in Appendix A.3 one may consider the apparent reflectance 
as being roughly proportional to the extended mass scattering efficiency t(K, ji) s c(X, 1)^1. One then expects the 
enor rates to scale by a factor of the ordCT of \i for otherwise equivalent conditions (at least for small to intomediate 
values of the observer's zenith angle where non linear multiple scattering effects and non spherical geometry 
conditions are less impoitant). It is advantageous then to favour larger zenith angle measu rements provided dve errors 
Il jgmselves do not covarv with zenith angle. In principle this is not the case unless the input paramters to the 
inversion procedure are themselves determined from non nadir inversions of the satellite data. 

At large surface reflectances the non linearities induced by the competing effects of scattering and attenuation suggest 
that alternative methods for sulfate inversion should be investiga^. The use of apparent reflectances d ifferentials 
over increments of viewing ancle bears investigation (c.f. Figure 3.4.6.3). This approach which reljttes zenith angle 
dependent apparent reflectance differentials with sulfate concentration is well suited to NCAA and SEAWIFs data 
whose overlapping images at large viewing angles (inter image frequency of about 100 minutes) are jarticularily 
amenable to studies of angular variation. In addition, because equivalent pixels can be selected as tai^ets, the 
problems of surface reflectance calibration are considerably reduced. The excellent signal to noise characteristics of 
the SEAWIFs sensor enhance this multi-image capability. 
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(b) Solar Geometry 

As discussed in Section 3.4.6 the effects of variations in solar zenith angle are analogous to variations in the 
observer's zenith angle. Accordingly the same statement concerning larger zenith angles applies in the case of solar 
geometry as well. For sun synchronous orbits such as those associated with the NOAA, Landsat and SPOT satellite 
sensors the solar zenith angle for a charactoistic 1500 GMT midlatitude crossing varies between a minimum of 30 
degrees in summer to 70 degrees in winter. The meteorological satellites (GOES and Meteosat) which acquire images 
at half hour intervals arc subject to solar zenith angle variations from 90 degrees to a minimum depending on the 
latitude and solar declination (time of year). 

4.4.2 Wavelength Passband 

In terms of first order models such as the one presented in Appendix A,3 the sulfate error rates vary as 1 / c in 
atmospheric states well removed from the transition point between scattering and attcnufttion dominated atmospheres. 
As such the error rates increase with wavelength roughly as (K / Xref)" where a is the commonly employed 
Angstrom coefficient (see Shaw et al. , 1972 for example) and A^f is a reference wavelength (0.55 ^m in our case). 
The value of the Angstrom coefficient is typically of the order of 1.3 ( in comparison the Lcaitch pure sulfate log 
normal model of Table 2.1.1.1 yields a value near 2.2 for wavelengths between .4 and .860 \xm). 

This spectral tendency which reflects a decreasing sensitivity with increasing wavelength was observed for error rate 
calculations at different wavelengths (not presented in the text) and is qualitatively evident in apparent reflectance 
computations such as the satellite passband dependent curves presented in Figure 3.4.5.1. Clearly however the 
i ncreasing su l fate error rates cannot be the only ontimizatinn criteria since dec reasing renectance noise valii<-; ^n fl 
decreasing Uncen ai ntV in t he s urf a ce reflec tance ca n e ff ectively offset the en-nr rate degradation. This offsetting efifeci 
in terms of a decrease in reflectance noise is indeed the case for the visible bands of the TM sensor and for SEAWIFS 
(Table 2.4.5.1). In terms of surface reflectance considerations the decreasing reflectance of water pixels with 
increasing wavelength effectively induces a decreasing uncertainty (c.f. Table 4.2.1.1) which can be effectively taken 
advantage of in the near ER spectral region. Of the bands presented in Table 4.2.4.1 the SEAWff's band 4 is clearly 
one of the most relevant in regards to exploiting the above considerations. 
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CHAPTER 5 
GENERAL CONCLUSIONS AND RECOMMENDATIONS 



We list here the general conclusions and lecommendations extracted from the radiative transfer simulations and the 
computations of error rates. The reader is referred to the main body of the text for details. The principal reference 
source in the text is given in parenthesis. 

o Inversion algorithms most be referenced to low or high surface reflectances well removed from the 
Kaufinan critical reflectance (Section 4.2.7} 

o The measured satellite signal is reasonably independent of the carbonaceous soot component over low 
reflecting surfaces and strongly dependent over high reflecting surfaces. This infers an inversion scheme 
would be best optimized by concentrating or weighting the sulfate inversion over low reflectance 
surfaces and the soot inversion over highly reflecting surfaces (Section 4.2.S). 

o Heavy mass loadings of dust can be discriminated from sulfates if two or more wavelengths in the 
visible and near IR (.4 to 1.S |im) are av^Iable and provided the sensor noise specifications are 
adequate (Section 4.2.5). 

o In principle 10 bit radiometric resolution is the minimum requirement for sensors not optimized to the 
low intensity levels characteristic of atmospheric scattering phenomena. Of the sensors inv^tigated 
only the NCAA series and the SEAWIFs sensors meet these specifications. Coarser resolution 
sensors which are optimized METEOSAT and GOES are restricted by saturation over highly reflecting 
targets. 

o For a nadir look direction, a wavelength of 0.S5 pin and a range of sulfate concentrations between 

and SO |ig/m3 a cumulative precision in the order of S jigMiS is obtained given the following 
constraints (Section 4.3): 

- fractional sulfate scattering efficiency known ~ (2.8 ng-m-3 / SO4") (quadrature sum 

of 3 independent errors). This is equivalent to the conditions (1.6 iig/m^ error per 
component): 

(a) relative humidity known to - 3% @ 70% rJi. and ~ .3% (3> 90% r.h. 

(b) dry eOfective radius smaller than .02 |jm 

(c) sulfate fractional numbo' density >- .OS 

• surface reflectance known to - .002. Presently this is only possible for waier 

surfaces. The most promising candidate surface other than water is vegetation in the 
chlorophyll absorption band centered near .65 |im. 

- single scattering albedo known to better than - .05 for a low reflectance surface. The 
accuracy requirement for a highly reflecting surface is prohibitive 

- fractional number density of sul&tes known to - .OS in the context of a pure 
sulfaie/dust atmosphere 

o At non nadir observing angles the atmospheric scattering is more responsive to sulfate presence. The 
sulfate error rates are correspondingly decreased as are the computed errors, if the input parameter 
uncertainties do not themselves covary with the observation geometry (Section 4.4.1). The NOAA and 
panic uiarily SEAWIFs sensors are well suited to non nadir investigations because of their large swath 
width and because the overlap obtained permits multi-angle investigations of the same pixel. 



72 



At large surface reflectances and non nadir angles the satellite signal is non linear with sulfate 
concentration. In this case, methods of apparent reflectance differentiais over zenith angle bear 
investigation (Section 4.4.1(a) ). 

The preferred operating wavelength (or the key wavelength in a tnultispectral inversion scheme) is 

constrained primarily by error rates which vary roughly as X°- and the spectral form of the surface 
reflectance. For water this implies an optimum wavelength in the red/near IR region of the 
spectrum (Section 4.4.2). A band centered in the chlorophyll absorption minima at .65 jim could be 
wnployed in an inversion procedure over vegetation surfaces (Section 4.2.1). 

If the relative humidity effects on the mass scattering efficiency can be specified sufficiendy accurately 
then a high relative humidity atmosphere is the preferred state for measurements of sulfate 
concentradon (Section 4.3). 

The 5S model is sufficiendy accurate for nadir computations but must be improved for large off nadir 
angles (Section 4.2.6) 

High spectral resolution imaging systems such as the EOS sensors HIRIS and MODIS can be 
exploited to detect molecular absorbers by differential absorption techniques. The high spectral 
resolution capability is less effective for scattering based phenomena where the spectral information 
content is coarse (Section 4.2.5). 

The achievement of the inversion precision targeted above requires at the very least a coordinated 
summer and winter campaign of ground based measurements tied to a multi-temporal set of NOAA or 
preferably SEAWrs overflights. The components of the ground campaign would include: 

- sulfate mass concentration measurements 

- nephelometer measurements 

- selection and reflectance calibration of standard homogeneous test sites (ground based 
and airborne reflectance measurements) 

- carbonaceous soot measurements 

- multi-spectral sunphotometer measurements 

- fine piffticle size distribution measurements 

- measurements of relative humidity and other meteorological parameters (pressure, 
temperature, wind magnitude and direction) 



Since the lion's share of these measurements are already performed in existing networks (CAPMON 
network, AES sunphotometer network etc. ) a ccnsiderable part of the wrak would involve 
coordinating and collating network measurements with supplementary measurements and with the 
satellite data. 
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ERRATUM 



During the editing of the above report it was realized that a small but persistant error was imbedded in the 
calculations. This error does not influence any of the conclusions but would be understandably annoying for anyone 
trying to duplicate our calculations. The source of the error stems fi?om the manner in which the aerosol optical depth 
was entered into the 5S model. The value actually entered was the aerosol scattering optical depth while the input 
required by 5S is in fact the total aerosol optical 6epth. These two values generally differ only slightly at least in the 
visible portion of the wavelength spectrum. In any case, as explained below, the error is effectively set to zero by 
reinterpreting the number densities or concentrations used as inputs to the optical and radiative transfer calculations. 

The aerosol scattering optical depth (c.f. Appendix A.1) is expressed as: 

where a, s, c, and d refer resp«;tively to aerosol, sulfate, carbonaceous soot, and dust like, and where N ref^s to the 
particle number density in units of cm-3. If ta*^* is assumed to be a total aerosol optical depth then one can adjust 
the number density definitions such that it in fact does represent a total aerosol optical depth viz: 

"ta" = TaS'=a = Cls (asS'=a/as)Ns + ac(CTc^a/ac)Nc + Od(adSca/ad)Nd 

= CTs (0s«=«/cts)Ns + ac(ac«='/CTc)Nc + ad(CTdsca/ad)Nd 

= as (coo.sNs) + ac(0Jo,cNc) + CTd(a)o,dNd) 

where the Mo.i*s are the pure single scattering albedo's for each component. This last expression indicates that the 
nominal aerosol scattering optical depths entered can be viewed as total aerosol optical depths if the number densities 
are modified by a factor equal to the pure single scattering albedo for the component of interest Table 6.1 lists the 
precise values of these modificadons. In general the sulfate concentration changes by less than 3% while the soot 
concentration decreases by a factor of the order of two. This latter decrease is the only significant change as far as the 
simulation studies are concerned. In those few cases (chapter 3 to 5) where references are made specifically to the soot 
concentration, an appropriate modification based on the above argument is included in the text and referenced to this 
section. Otherwise no attempt was made to perform a general correction according to the above relation. If more 
precise calculations are required the reader need only apply the corrected values in Table 6.1. 
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Table 6.1 Corrected number densities: Ml (cm**<3) 



a) Values used initially as per Table 2.1.4.3 

<JVSt-liK9 Sfifll I sulfate 

C.V. 2.3229 74.1400 | 

U.V, 0.2682 775.4200 1 430 

U.S. 0.4646 14.8280 j 860 

C.S 0.0536 155.0850 j 1290 

I 1720 

I 2150 

I 2580 

I 3010 

I 3440 

I 3870 

I 4300 



b) Effective values if the scattering aerosol optical depth entered is treated as total aerosol optical 
depth 

I 





dust-like 

1.69172 


soot 1 


0.0% 


50.0% 
0.00 


70,0% ,. 

0.00 


80.0% 


90.0% 
0.00 


95.0% 
0.00 


C.V. 


32.3263 1 


0.00 


0.00 


U.V. 


0.19532 


338.0963 1 


414.9 


417.22 


420.40 


421 .98 


424.83 


427.07 


C.S. 


0,33834 


6.4653 1 


829.82 


834.45 


840.80 


843.95 


849.66 


854.14 


U.S. 


0.03906 


67.6197 1 


1244.74 


1251.67 


1261.20 


1265.93 


1274.48 


1281.21 








1659.65 


1668.89 


1681.60 


1687.90 


1699.31 


1708.28 








2074.56 


2086.12 


2102.00 


2109.88 


2124.14 


2135.35 








2489.47 


2503.34 


2522.40 


2531 .85 


2548.97 


2562.41 








2904.38 


2920.56 


2942.80 


2953.83 


2973.80 


2989.48 








3319.30 


3337.79 


3363.20 


3375.81 


3398.62 


3416.55 








3734.21 


3755.01 


3783.60 


3797.78 


3823.45 


3843.62 








4149.12 


4172.23 


4204.00 


4219.76 


4248.28 


4270.69 



-yvmVT'wy • ' VTW*! 



REFERENCES 



Ackerman, T. P., O. B. Toon, (1981), Absorption of Visible Radiation in Atmosphere Containing Mixtures of 
Absorbing and Nonabsorbing Particles, Applied Optics, 20, pp. 3661-3668. 

Ahem, F. J., Goodenough, D. G., Jain, S. C, Rao, V. R., G. Rochon, (1977), Landsat Atmospheric Correction at 
CCRS, 4th Can, Symp. on Rem. Sens., Qudbec City. May (1977). 

Ahmad, Z., Sutton J., Stowe, L., C. R. N. Rao. (1988), A Regression Method for Estimating Atmospheric Aerosol 
Properties from Current and Modified Versions of the NOAA/AVHRR Instrument, InL Radiation Symp.. LiUe, 
France, 18-26, August 1988. 

Baker, J. (chairman), (1987), System Concept for Wide-Field-of-View Observations of Ocean Phenomena Crom 
Space. Report of the Joint EOSAT/NASA SeaWifs Working Group, NASA. Washington D.C., EOSAT, Lanham 
Md. 

Bergstrom, R. W.. (1973). Comments of the Estimation of Aerosol Absorption CoefTicicnts in the Atmosphere, 
Contributions to Atmospheric Physics (translation of Beitrage zur Physik der Atmosphare), 46, pp. 198-202. 

Begni. G.. Henry, P., Lcvoy, M., et M. Dinguirard (1987), La qualitd des images SPOT, Compte rendu du 
Colloque SPOT-1, Utilisation des images, Bilan etrdsultats, Paris, France, novembre 1987, pp. 1505-1515. 

Carlson. T. N.. Wendling. P., (1977), Reflected Radiance Measured by N0AA3 VHRR as a Function of Optical 
Depth for Saharan Dust, J. App. Meteor.. 16, pp. 1368-1371. 

Chandrasekhar, S., (1960), Radiative Transfer, Dover Pub. Inc., New York. 

Colwell, R. N., (1983), Manual of Remote Sensing, American Society of Photogrammetry, Sheridan Press, Fails 
Church, Virginia. 

Dalu. G.. (1986), Satellite Remote Sensing of Atmospheric Water Vapour, InL J. Rem. Sens., 7, pp. 1089-1097. 

Dedieu, G., Deschamps, P. Y., and Y. H. Kerr, (1983), Solar Irradiance at the Surface from Meteosat Visible Data, 
Machine Data of Remotely Sensed data Symposium, pp. 1-4. 

Deiimendjian, D. (1969), Electromagnetic Scattering on Spherical Polydispersions, American Elsevier, New York, 

Deschamps, P. Y.. L Phulpin, (1980), Atmospheric Correction of Infrared Measurements of Sea Surface 
Temperature Using Channels at 3.7. 1 1, and 12 mm. Boundary-Layer Meteorology, 18, pp. 131-143. 

Deschamps, P. Y.. Homan. M.. D. Tanrd, (1983), Modeling of the Atmospheric Effects and its Applications to the 
Remote Sensint of Ocean Colour, Applied Optics, 22, p. 3751-3758. 

Deuze, J. L., Devaux, C. Herman, M., Santer, R., D. Tarn*. (1988), Saharan Aerosols over the South of France: 
Characterization Derived from Satellite Data and Ground Based Measurements, J. of Applied Meteor., 27, p. 680- 
686. 

Egan, W. G., (1986). Satellite Remote Sensing of the Optical Properties of the Sahara and the Climatological 
Significance, 20ih Inter. Symp. on Rem. Sens, of Envir.. Nairobi, Kenya, pp. 95-114. 



-'*»»^ npcn 



76 



Elterman, L., (1970), Vertical-Attenuation Model with Eight Surface Meteorological Ranges 2 to 13 kilometers, Air 
Fbrce Cambridge Research Laboratories, Bedford, Mass. 

Fraser, R, S., (1986), Satellite Measurements of Mass of Sahara Dust in the Atmosphere, App. Optics, 15, pp. 

2471-2479 

Fraser R. S. , Kaufman, Y. J., (1984), Satellite Measurements of Aerosol Mass and Transport, Atmospheric 
Environment. 18, pp. 2577-2584. 

Freemantle, J., Miller, J. R., K. S. Jammu, (1984), Atmo^heric Aerosol Variability over the Great Lakes fix)m 
NOAA-7 AVHRR Imagery. Proc. of the 9th Can. Symp. on Rem. Sens., Sl John's, T'M.D, 13-17 Aug. 1984, pp. 
263-270. 

Frouin, R., Gauthier. C, Calibration of NOAA- AVHRR, GOES-5 and GOES-6 VISSR/VAS solar channels, Rem. 
Sens, of Envir., 22, pp. 73-101. 

Griggs, M., (1987), Satellite Measurements of Tropospheric Aerosols, NASA CR-3459, NASA-15898, 58p. 

Gordon. H. R., (1978), Removal of Atmospheric EffecB from Satellite Imagery of the Oceans, Atmlied Optics, 17, 
pp. 1631-1635. 

Hanel, G., (1976), The Properties of Atmospheric Aerosol Particles as Functions of The Relative Humidity at 
Thermodynamic Equilibrium with the Surrounding Moist Air, Advances in Geophysics, Academic Press pp 73- 
188. 

Hay, J. E., (1971). IVeciptable Water Over Canada, AtmosphCTe, 9, pp. 128-143 

Hindman, E. E.. Durkee, P. A., W. Kuenning, (1985), Properties of Aerosol Particles Detected by Satellites in 
Coastal regions, J. Rech. Aonos., 19, pp. 315-322. 

Husar, R. B., D. E. Patterson, (1980), Regional Scale Air PoUution: Sources and Effects, Annals New York 
Academy of Sciences, pp. 399-417. 

Isaac, G. A.. Leaitch, W. R., J. W. Strapp, (1988), The Vertical Distribution of Aerosols and Acid Related 
Compounds in Air and aoudwatcr, submitted to Atmospheric Environment. 

Ivanchik, M. V., Klyvslinikov. S. I., Krovotynisev, V. A.. Serebrennikov, A. N.. (1985), Remote Sensing of Dust 
Formations in the Atmosphere, Sov. J. Remote Sensing. 4. pp. 375-380. 

Japar, S. M.. Brachaczek, W. W., Gorse R. A., Jr., Norbeck, J. M., W. R. Pierson. (1986). The Contribution of 
Elemental Carbon to the Optical Properties of Rural Atmospheric Aerosols, Atmospheric Environment. 20 dd 
1281-1289. ' *^' 

Jennings, J.G. J.B. Gillespie (1978), Mie Theory Sensitivity Studies - The Effects of Aerosols Complex Refractive 
Index and Size Distribution Variations on Extinction and Absorption Coefficients part 11: Analysis of the 
Computational Results, ASL report, VS Army Electronics Research and Development Command. 

Kaestner, M., Koepke, P., H. Quenzel, (1982), Monitoring of Saharan Dust over the Atlantic Using Meteosat-VIS 
Data, XXIV COSPAR, Ottawa, Canada 

Kattawar, G. W., G. N. Plass. (1967), Electromagnetic Scattering from Absorbing Spheres, Applied Optics. 6. dd. 
1377-1382. p . .FK 



77 



Kaufman, Y. J., J. H. Joachim, (1982), Determination of Surface Albedos and Aerosol Extinction Characteristics 
from SateUite Imagery, J. Geoph. Res.. 87, pp. 1287-1299. 

Kaufman, Y. J., R. S. Fraser, (1983), Light Extinction by Aerosols during Summer Air Pollution, Journal of 
Climate and Applied Meteorology, 22, pp. 1694-1706. 

Kaufman, Y. J., (1987), Satellite Sensing of Aerosol Absorption, Jour. Geophys. Res., 20, pp. 4307-4317. 

Kaufman, Y. J., C. Sendra, (1988), Algorithm for automatic atmospheric concctions to visible and near-ER satellite 
imagery. Int. J. Remote Sensing, 9, pp. 1357-I38I. 

Kerr, J. B., McElroy, C. T., Olafson, R. A., (1980), Measurements of Ozone with the Brewer Ozone 
Spectrophotometer, Proc. Inter. Ozone Symposium, Boulder, Col., pp. 74-79. 

Kerr. J. B., W. F. J. Evans, (1980), Monitoring of Sulphur Dioxide at Downsview with the Brewer Ozone 
Spectrophotometer. Internal Report APRB 1 1 1X35, Atmospheric Environmental Research, Environment Canada 

Kim, Y. J., Sievering, H., J. F. Boatman, (1988), Airborne Measurement of Atmospheric Particles in the Lower 
Troposphere Over the Cen&al United States, Jour. Geophy. Res., 93, pp. 12631-12644. 

Lauritson, L., Nelson, G. J., F. W. Porto, (1979), Data Extraction and Calibration of Tiros-N/NOAA Radiometers, 
NOAA Technical Memorandum NESS- 107, NOAA, US Dept. of Commerce. 

Legrand. M., Bertrand J. J., M. Desbois, (1985), Dust Clouds over West Africa: A Characterization by Satellite 
Data, Ann. Geophys., 3, pp. 377-786. 

Legrand, M., Desbois, M., VORN?, K., (1988), Satellite Detection of Saharan Dust: Optimum Imaging During 
Nightime, J. Climate, submitted for publication. 

Leitch, W.R., G.A. Isaac (1988), Tropospheric Aerosol Size Distribution in Eastern and Northern North America 
Between 1 000 and 500 Mb: Signiuures for Condensation and Stratospheric Intrusions, Atmospheric Environment 
Service internal report 

Leaderer, B. P., J. A. Stolwijk, (1980), Optical Properties of the Urban Aerosol and Their Relation to Chemical 
Composition, Annals New York Academy of Sciences, pp. 70-85. 

Lyons, W. A.. (1980), Evidence of Tansport of Hazy Air Masses From Satellite Imagery, Annals of New York 
Academy of Sciences, pp. 418-433. 

Macias. E.S., Blumenthal, DX., Anderson, J.A., B.K. Cantrell (1980), Size and CcHnposition of Visibility- 
Reducing Aerosols in Southwestern Plumes, Annals of New York Academy of Sciences, pp. 233-257. 

Malila, W. A., D. M. Anderson, (1986), Satellite Data Availability and Calibration documentation for Land Surface 
Climatology Studies. International Satellite Land-Surface Qimatology Project (ISLSCP), North American ISLSCP 
Retrospective Analysis Program (IRAP). ISLSCP Report No. 5, University Coiporation for Atmoq)heric Research, 
Boulder Colorado. 

Meckler, Y., Quenzel, H., Ohring, L., I. Maums, (1977), Relative Atmospheric Aerosol Content from ERTS 
Observations, J. Geophys. Res., 82, pp. 967-970. 

Middleion. W. E. K., (1952), Vision Through the Atmosphere, Univwsity of Toronto Press, Toronto. 

Morel, A., L. Prieur, (1977), Analysis of Variations in Ocean Color, Limnology and Oceanography, 22, pp. 709- 
722. 



78 



O'Neill. N. T., Miller, J. R., F. J. Ahem, (1978), Radiative Transfer Calculations for Remote Sensing 
Applications, 5th Can. Symp. on Rem. Sens., Victoria, August (1978). 

O'Neill, N. T., (1982), A Radiative Transfer in an Urban Atmosphere, Ph. D. thesis, York University, Toronto. 
Ont. 

O'Neill. N. T., J. R. Miller, (1988), On Calibration of Passive Optical Bathymetry through Depth Soundings: 
Analysis and Treatment of Errors Resulting from the Spatial Variation of Environmental F^rmaters, accepted for 
publication in the Int. Jour, of Rem. Sens. 

Otternian, J., Eraser, R. S., O. P. Bahethi. (1982), Characterization of Tropospheric Desert Aerosols at Solar 
Wavelengths by Multispectral Radiometry from Landsat, J. of Geophys. Res., 87, pp. 1270-1278. 

Penndorf, R., (1957), Tables of the Refractive Index for Standard Air and the Rayleigh Scattering Coefficient for the 
Spectral Region Between 0.2 and 20.0 microns and their Application to Atmospheric Optics, J. Optical Soc. of 
America, 47, pp. 176-182. 

Picrson, WJl., Brachacczek, W.W., Truer, T.J.. BuUer. J.W.. TJ. Komiski (1980), Ambient sulfate Measurements 
on Alleghenny Mountain and the Question of Atmospheric sulfate in the Northeastern United States, Annals New 
York Academy of Sciences, pp. 145-173. 

Richards, L. W.. R. W. Bergstrom Jr., (1986), The Optical Effects of Fine-Particle Carbon on Urban Atmospheres. 
Atmosphere Environment, 20, pp. 387-396. 

Roosen, R. G., R. J. Angione, C. H. Klemcke, (1973), Worldwide Variations in Atmospheric Transmission: 1. 
Baseline Results from Smithsonian Observations, Bulletin American Meteorological Society, 54, pp. 307-315. 

Rosen, H., A. D. Hansen, L. Gundel, T. Novakov, (1978), Identification of the Optically Absorbing Component in 
Urban Aerosols, Applied Optics, 17, 3859-3861. 

Rosen, H., T. Novakov, (1983), Optical Transmission Through Aerosol deposits on Diffusely Reflective FQters: a 
Nfcthod for Measuring the Absoiting Component of Aerosol Particles. Applied Optics, 22, pp. 1265-1267. 

Royer, A.. Bemier, R., Adler, U.. E. Pilon, (1985), Calcul de diffusion de la lumi&re par Ics aerosols: Descriptions 
dcs programmes "MDELOR" et r&ultats obtenus. 

Royer, A.. Charbonneau, L., Brochu, R., Murphy, M. M., P. M. TeUlet, (1987), Radiometric Comparison of the 
Landsat-5 TM and MSS sensors, InL J. Remote Sensing, 8. pp. 579-591. 

Royer, A,. Charbonneau, L., (1988a), Interannual Landsat-MSS Reflectance Variation in an Urbanized Temperate 
Zone, Remote Sensing of Environment, 24, pp. 423-446. 

Royer, A., Davis, A., N. T. O'Neill, (1988b). Analyse des effets atmosph^riques dans les images HRV de SPOT, 
Canadian Journal of Remote Sensing. 14, pp. 80-91. 

Ruby, M. G.. A. P. Waggoner, (1981). Intercomparison of Integrating Nephelometer Measurements. Environmental 
Science and Technology. 15, pp. 109-1 13. 

Shaw, G. E., J. A. Reagan, B. M. Herman, (1972), Investigation of Atmospheric ExUnction Using Direct Solar 
Radiation Measurements with a Multiple Wavelength Radiometer, Journal of Applied Meteorology, 12, pp. 374- 

Slater. P. N., (1980), Remote Sensing: Optics and Optical Systems, Addison- Wesley, Reading Mass. 



'(~-^y;'. - —IT' ' •?^^,'^" ' l^l,mw 



Sloane, C.S. (1983), Optical Properties of Aerosols - Comparison of Measurements with Model calculations. 
Atmospheric Environment, 17, pp. 409-416. 

Sloane. C.S. (1984), Optical Propenies of Aerosols of Mixed Composition, 18. pp. 871-878. 

Sloane, C.S. (1986), Effect of Composition on Aerosol Light Scattering Efficiencies, Atmospheric Environment, 
20, pp. 1025-1037. 

Shcttle, E.P., R.W. Fenn (1979), Models for the Aerosols of the Lower Atmosphere and the Effects of Humidity 
Variations on their Optical Properties, AFGL-TR-79-0214, Air Force Geographies Laboratory, Hanscom AFB, 
Massachusetts. 

Syed, M. Q., A. W. Hsvrison, (1980), Ground Based Obsovations of Stratospheric Nitrogen Dioxide, Can. J. Phys., 
58, pp. 788-802. 

Takayama, Y., T. Takashime, (1986), AmxjsoI Optical Thidoiess of Yellow Sand over the Yellow Sea Derived from 
NOAA Satellite Data, Atmos. Envir.. 20, pp. 631-638. 

Tang, I. N.. Wong, W. T., H. R. Munkelwitz, (1981), The Relative Importance of Atmospheric Sulfates and 
Nitrates in Visibility Reduction, Atmospheric Environment. 15, pp. 2463-2471. 

Tanrd, D.. Deschamps P. Y.. Devaux C, M. Herman, (1988), Satellite Measurements of the Saharan Aerosol 
Optical Thickness over Land Surface, IRS88. 

Tomasi, C, (1979), Weak Absorption by Atmospheric Water Vapour in the Visible and Near-Infra-Red Spectral 
Region, 11 Nuovo Cimcnto, 2C, pp. 511-526 

Toon, O.B., PoUackJ.B., BJsT. Khare (1976), The Optical Constants of Several Aunospheric Aerosol Species: 
Ammonium Sulfate, Aluminium Oxide, and Sodium Chloride, Journal of Geophysical Research, 81, No 33, pp. 

5733-5748. 

Tsonis. A., (1987), A Method for Inferring Aerosol Properties From Satellite Data over Forested Terrain, 
Atmospheric Environment, 21, pp. 2467-2471. 

Tsonis, A., WJ(. Leaitch (1986), Minimum Detectable Pollution Levels from Satellite Imagery, Geophys. Res. 
Lett., 13, pp. 56-59. 

Turner, R. E., M. M. Spencer, (1972). Atmospheric Model for the Correction of Spacecraft Data. Proc. of the 8th 
Int Symp. on Rem. Sens, of Environ., Ann Mtot, pp. 895-934. 

Vossier, TX.. E.S. Macias (1986), Contribution of Hne Particle sulfates to Light Scattering in Sl Louis Summer 
Aerosol, Environmental Science and Technology, 20, pp. 1235-1243. 

Young, A. T., (1981), On the Rayleigh-Scattering Optical Depth of the Atmosphere, J. App. Met, 20, pp. 328- 
330. 

WCP-55, 1983, World Climate Research Program Report of the Experts Meeting on Aerosols and Their Climatic 
Effects, WMO. Williamsburg. 28-30 March 1983, 107 p. 

Wickramasinghe, N. C, (1973), Light Scattering Functions for Small Particles and Applications in Astronomy, 
Wiley, New York. 

Wolfe, L. W., G. J. Zissis, (1978), The Infrared Handbook. Prepared by the Environmental Research Institute of 
Michigan for die Office of Naval Research, Deiartment of the Navy, Washington DC. 



I 
I 



APPENDIX A 



A.l Symbols and Deflnitions 

A. 2 Relative Humidity Modelling 

A. 3 Some Useful Approximations 



APPENDIX A.1: SYMBOLS AND NOMENCLATURE 

Optical Parameters for a Given Aerosol Component 

The size distribution of a heterogeneous collection of aerosol particles is commonly rcfexred to in terms of the 
number density per cm"^ pet tmit radius int«Tral viz: 

iJiW = NiHifr) 

where J n;(r)dr = N. 


is the total number per cm"^ and hence where the normalized size distribution has the property: 



J n.(r)dr = 1 



t 

The total extinction cross section, assuming independent scattering, is the sum of all Mie cross sections for aU 
particles in the distribution: 



Oj(X) = J <^mj, a)n j(r) dr 



where m j is the refractive index and a s 2ia / A. is the size parameter. 

The contributors to the extinction cross section are pure absorption and outscattering from the incident beam: 

where 





Oj""(»^ j <^*(mj,a)n.(r)dr 



are the particle size distribution averaged absorption and scattering cross sections respectively. Multiplying by the 
total number of particles per unit volume yields the volume extinction, scattering and abs<xption coefficients: 

K^Q.) = Niai(X) 

KfU) =Nia**(X) 
In terms of the Mie scattering nomenclature the conservative scattering phase function is given by: 



i,(mj.a,cosx) + i,(m., a,cosx) i 
pjnij, a,cosx) = 



2 „2, 



cr"(m., a ) 
where (^ = ^——^^ 
a(mj.a) 



is the scattering effeciency and il and i2 (or i^ and i|, respectively) are the standard Mie intensity functions (Bora and 

Wolf, 1965, Wickramasinghe, 1973). The conservative scauering phase function is constrained by its integral over 
solid angle viz: 



J p^(nip a, cos X ) — = 1 
A 4»t 



The average conservative phase function over the size distribution is given by the optical avaage: 



PoiCcos X ) = J P„(mj. a, cos x ) n ^(r) o»°^m., a) 



^ o^Q.) 

J p^(mj,a,cosx)%(r)o*^(nij.a) * 



^Q^) 

where an integration over solid angle clearly yields: 

f dQ 

I Poi(cosx)— = 1 
4ji 

The phase function definition which includes the effects of absorption is simply: 

p (ffij, a, cos X ) = (ojim., a ) p^(m., a, cos x ) 
where: 

©^(mj.a) = 

©(mj.a) 

is temed the single scauering albedo. The optically averaged representation of this phase function is: 

Pi(cosx) = J p(mi,a.cosx)nj(r)o(mj, a)-^ 
^1 (^) 

= Jp(rn.. a.cosx)n.(r)o(m,.a)-i- 
o 'q(X) 

or hence: 

Pj(cosx) = aJX) P„i(cosx) 



where the PSD averaged single scattering albedo is given by: 

The Lengendre polynonnial representaion of the Mie phase function is writtea: 

p (mj. a, cos X ) = jl^^A^v^^^ii'^'^X) 

where P^ (cos x ) are Lcgendre polynomials (Ambramowiiz and Stegun, 1972). Integrating over the size distribution 
yields the optically averaged representation: 



Pi (cos X ) = X ®^i W P^ (cos X ) 



where: 

f A d" 
Q) (X) = J to (nij, a)nj(r)<<m., a) 

' opC) 
= I 0) (m., a)n.(r)<<m., a) 

Aggregate Aerosol Model Parameterization 

For aggregate collections of indivichial aerosol components one can define an aggregate nonnalized size distribution: 

n(r) = XYini(r) viae X Y^ = 1 

y, = ^ «1 N = ZN; 
With this definition we clearly have: 

J n(r)dr = 1 
Q 

The actual number density is then simply: 

i 

= NXrinVr) 



The volume of the PSD normaiized to a single particle is: 



The fractional volume associated with all the particles of component i is: 



t^ sVj/V = N;V./V 

The ftactional volume and fractional number density arc accordingly interrelated by the expressions given below: 

A 



<1 



Yi - 



ZlTiVi 



Ci/Vj 



Volume extinction and scattering coefTicients for the aggregate aerosol distribution are given by: 

i<X) = Na(X,) 

where: 

<J(^) = S^tJjC^) and a**(X) = ^y.a^CK) 

i i 

The optically averaged conservative phase function takes on the form: 

X«»"(X)p„i(cosx) 



Po(cosz) = 



ICS 

o (X) 
X "T^PoiCcosx) 



In a similar feshion as above the phase function is expressed as an optical mean weighted by the component 
extinction coefficients. This yields: 

p (cosx) = (o^(X) p„(cosx) 
where the aggregate single scattering albedo is given by: 

CO (X) = ^^ = J^ 
" g^ K(X) 



Tf 



Finally, the aggregate aerosol Lengdre formulaiion can be expressed as: 

p(cosx) = X 03.WP,(cosx) 

By integrating over solid angle and exploiting the orthogonality properties of the Lengendre polynomials one finds: 






o(X.) K(X) 



Log Normal Size Distributions 

A common and convenient means of describing the particle size distribution is in terms of the log nomial size 
distribution viz: 

6 



•^^"8^ yiiJiogc 



where the log is referenced to the base 10. The parameter rj^ is commonly called the geometric mean radius while ais 
the geometric standard deviation. As written this relation has the following mathem^cal properties: 

f (N 

dlogr = N 



J dlogr 



J '°8f -Ti ^ = logr^ 

J dlogr N " 



f n 1 v2_d^ dlogr , 2 
I aogr - logr^ — -A- = logCT 

J "^ dlogr N 

The form of these expressions is clearly independent of the logarithmic base for the (defined) ciut^c parameters i^ •■ 
lO^oS^N and a- 10^°8O'. The volume (or mass) distribution is given by: 

dV 4 3 dN 

= — itr 



dlogr 3 dlogr 
After some manipulation this can be written as: 



2 2 
dV 4 3 N - (logr - logry) /21og o 



= — Jt Tj — ^ e 



dlogr 3 ^y^ 



logo 



where rj is an effective mean volume radius per particle in as much as the integration of the above ejcpression yields: 



A 



1 f dV ., 4 3 

V = — dlogr = — jir- 

Ul dlogr * 3 3 

Because the volume distribution retains the log-nonnal form, ry can be interpretated as a geometric mean radius 
(weighted by the volume distribution) viz: 



f , dV dlogr 

I logr- -^ = logry 

_^ dlogr V ^ 

As well it will be noted, from the form of the exponential term for dV/dlogr above, that the log normal volume 
distribution retains the same standard deviation (log a) as the log normal distribution for the number density. The two 
volume based mean radii are related to the geometric mean radius r^j by: 

3 2 3 2 

logr3 = logr]^j + -InlOlog o ot lnr3 = lnrj^+-lno 

2 2 

log Ty = log Tjj + 3 InlO log c or In Ty = In r,^ + 31n o 

These expressions were derived by explicitly deducing the log normal fcwm for dV/dlogr above and thence equating 
terms. It will be noted that the log and In expressions are self consistent (i.e. derivable one from the other). 

In a similar fashion one can obtain analogous surface expressions. The form of these expressions is identical to the 
volume expressions except that one replaces S for V. rs for ry. r2 for r3. and n T^:n r^ for 4/3n T^-Af3m^'^. The two 
radii parameters are given by: 

2 2 

logrj =logrj^ + 21nl01og a or In Tj = In rj^ + 2hi a 

2 2 

log tj = log rj^ + InlO log o or In r^s In rj^ + In a 

In a good deal of the fundamental aerosol research the particle dimensions are reported in terms of diameter rather than 
radius. This implies some trivial changes in the above expressions which for the sake of completeness are noted 
here: 

dN tf^ N -(tegD-k,gDj^)/21og^a) 

— e 



d log D d log r r~ , 

y/lK logo 

whae log o is the same as above and where: 

r 1 cN 

logD^=| logD— --— — dlogD = log2 + logr^, 
/ N d log D " 



Columnar Parameters 

The total number of particles in a one cm^ atmospheric column between the radii r and r + dr per unit radius 
increment is tenned the di£Ferential abundance. Per a single aerosol component the diffeiential abundance is related to 
the differential number density by: 

aj(r) = J ni(r,h)(fli 


= n.(i)J Ni(h)dh 


. nj(r)Aj 

where Aj is the total abundance for all radii. Clearly: 

J ij(r)dr « Aj 


The extinction optical depth is expressed as: 

1^ = J KjOl. X,)(fll 


= Ai<jj(X) 

In a similar fashion one can write for the scattering optical depth (similar expression for the absorption optical 
dq)th): 


For the aggregate aerosol the differential abundance is clearly: 

•ft 

a(r) = J n(r.h)* 



= I N(h)nO-)dh 

= An(r) where A = J N(h)dh 



is the total number abundance for the aggregate aerosol. The aggregate aerosol optical depths are written: 



a 



t = J i<h,X)dh 



s Ao(X) 

T^= J K"(h.X)dh 


Parameterization of Altitude Pronies 

The number density altitude profile is given by: 

n. (e, r ) 

h(z,r) * 

nj(0.r) 

If the size disuibution is independent of altitude one can write: 

Ni(z) ic(z,X) 



hi(2.r) = h;(z) 



Ni (0) Kj (0, X) 

The equivalent column length referenced to the surface wmcentration is defined by: 

as 

H = fh,(z)dz - — ^ 
J ' N.(0) 

where Hj is termed the scale height represents the height of an atmospheric column if all the panicles of component 
i were compressed lo the number density at the surface. If h^iz) varies exponentially as exp(- hk ) then clearly k-1 is 
the scale height. The optical depth am be written: 

Tj . Ni(0)H.Oi(X) 
and the ag^gate optical depth as: 

t = N(0)X^(0)a,a)Hi 



APPENDIX A.2: RELATIVE HUMIDITY MODELS 
Particle Growth Model 

As per Hanel (1976) one assumes that the total volume of a wet particle is the sum of dry particle volume and the 
volume of the pure water condensed on it viz: 

V = Vo + V^ 
If pQ, Mo and p^, Mw "»« *« densities and masses of the dry particle and the water respectively we may write: 

X = 1 + fii^ 

\ Pw M^ 

By employing the equivalent radii r(0) = (3Vq('4ji)1^ and r = (3V/4jc)1/3 this expression can be transformed into an 
explicit equation for the wet radius: 

r = r(0)tl + ] 

The wet particle mass m^ is a strong function of the water activity a^ given by: 

2ov^ 

where f is the relative humidity (o s f < 1). <yis the surface tension on the wet particle surface, v^ the specific 
volume of pure water. T the absolute temperature in degrees Kelvin, and Rw the specific gas constant of pure wate*. 
An evaluation of the exponential term for room temperature conditions yields (Shettle and Fenn . 1979): 

, , Q001056 , 

a^ = fcxp[ — --] 

r(a^ 

According to Shettle and Fenn variations in curv^ure effect due to typical departures of atmospheric temperature finom 
rtwm temperature are less than 2%. Given the ccMipled wcpressions for a^ and r, input values for f and r(0), and 

representative values of the water activity depaidcm ratio \x^{3i^)f\i^ and the density ratio p^p^ one can iterativety 

solve for the wet particle equivalent radius (r) , 
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Variation of Refractive Index With Relative Humidity 

As indicated by Shettle and Fenn (1979) the re&active index of the wet particle is simply the volume weighted mean 
of the dry and wet refractive indices: 

m = — 

V 

= i\ + (ni, - n^) — 



= nV + («b • "V)[t] 



3 
1 

r 



where r is obtained torn the iterative procedures described in the fwevious section. 
Variation of Log Normal Parameters 

Except for very small particle sizes ( r < .01 pm) the variation of water activity with radius is quite weak (a^ is 
approximately given by the relative humidity f) and the particle growth ratio r / r(0) is practically constant. If the 
governing size distribution is log normal and if no significant contributions to the size distribution arise &x)m the 
very small particles then one may continue to describe the wet particle distribution by a log normal expression with 
an appropriate shift in r^ towards larger radii. This follows from the definition of ij^j as given above in die discussion 
of the log normal parameterization, viz: 



logr,(0 = J logr(f)-f-il2^ 
•' dlogr N 



= J IlogKO)+logi^]-*Ldlogl 
i 1(0) dlogr N 

= logrf,(0) + log A 

I(U) 

where r(0 represents the particle radius at relative humidity f. Employing the above relation for the wet to dry particle 
radii ratio yields: 

TnCO - r^(0)[l + ^^l"" 

The simple proportionality of the geometric mean radius (r^) ^^^^ substituted into the defining equation for the 
geometric standard deviation (c^ yields: 

log2o(f) = log2o(0) 
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Cumulative Particle Size Distribution Effects 
The volume averaged refractive index can be written: 

n^(dV - dVg) + mjdVo 



d "Qti^ — 



dv ("I - "^^0 



V V V 

or hence by integrating over volume and employing the log normal descriptiais of the dry and wet particle size 
distributions we obtain: 

•cn> = m^ + (mg - m^) 

where r3 is the effective mean volume radius per panicle as defined above in the discussion of the log normal 
function and its associated parameters. As noted above the calculations are considerably simplified since one is only 
required to find a functional form for rj^a^) and o(a^) as a function of relative humidity or a^. 

The definition of a functional form for rj^Ca^) and o(a^) necessitates the development of representative mean model 

for the mass ratio m^ / p^^ as a function of water activity and an evaluation of the density ratio p^ / p^ (the dependent 

parameters of the r / r(0) ratio above). Figures A.2.1 (a) and A.2.1 (h) show the variation of the mass ratio as a 
function of (1 - a^) for two different types of aerosols (urban Model 1 and rural Model 6 respectively) as measured by 

Hand (1976). 

An interesting phenomena observed in the analysis of Hanel's data was a type of hysterisis effect wherein the 
variaticHi of the mass ratio versus water activity was different depending on whether the measurements were made in 
conditions of decreasing or increasing water activity. In practice this phenomenon manifests itself as an abrupt 
increase in particle size with increasing water activity (deliquescence point) foUowed by a more gradual increase as the 
log of (1 - with increasing relative humidity (see Tang ct al., 1981 for example). If however the relative humidity 
is decreased down to a point below the deliquescence point (79.5 % relative humidity for (NH4)2S04 particles) a 

state of supersaturation lingers such that a significant humid mass exists well below the deliquescence point. 

In the light of die above remarks on the lack of dependence between water activity and particle radius (and as 
evidenced by the resemblance of the two different water activity curves) we arbitrarily selected the urban Model 1 
curve tabulations of Hanel to represent the water activity dependence of the Leaiich water soluble aerosol (ci. Section 
2.1.1) along with a dry density of 1.8 g/cm^. The latter value represents a density typical of sulfate particles (see for 
example Sloane, 1984). 

In order to facilitate the incorporation of the watCT activity data into our particle prtwth model some simple empirical 
relations were deduced for the data of Figures A.l .2 (a). These relations served as well to produce mean values of 
mass rado at those values of water activity characterized by an uncertainty error bar. In the case of the region below 
the deliquescence point this approach produced dsi effective average between the states of supersaturation and pre 
dehquescoice. 
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These relations are given by; 

M^/Mo = 0.0217a^ + 0.3753^2 + o.08333a^3 0.00 <a^< 0.30 

= 0.00810 exp( 5.6091 a^ ) 0.30 <a^< 0.70 

= 0.1083 (l-aw)-^-^21 0.70 < a^ < 0.997 

Effective Relative Humidity Over Altitude 

In [he main body of the text we associate relative humidity with the ground level value commonly measured by 
meteorological stations. In this section an effective relative humidity over altitude is defined and shown not to be 
substantially different from the ground level value. 

The relative humidity is defined as: 

fsevCD/e^CT) 

where ev is the ambient vapour pressure aid e^ is the saturation vapour pressure. The saturation vapour pressure is 
given by the empirical relation (AES, 1987): 

e^ni = aexpilill^l 
' c+(T-273) 

where T is in degrees Kelvin and CjCT) is in milli barrs (a = 6.1 121; for T > 273, b = 17.368, c = 238.88 while for T 
< 273, b = 17.966, c = 247.15). The ambient vapour pressure can be expressed in terms of the ideal gas law as: 

ey(T) = .01R^p^T 

where Ry is the gas constant normalized to the molecular weight of water [.461 J / (Kg-°K) ] and p^ is the water 

vapoiff density [g / m^]. These parameters were computed for two standard midlatitude atmospheric models (summer 
^d winter versions as per McClatchey et al., 1971) in order to demonstrate the variation of relative humidity with 
altitude and with season (Table A.2.1). 

The pertinence of relative humidity for atmospheric scattering problems stems from its influence on aerosol particle 
growth. An effective relative humidity should be an appropriately weighted avoage over altitude which accounts for 
aerosol optical activity. Accordingly, we defined the weighted integral: 



W = 



Jf(h)K^)dh 


J K**(h)* 
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Table A. 2.1 Relative humidity variations with altitude for standard-r 

Altitude Pressure Temperature Sat. vapour Water vapour Rel. humidity 

pressure density 

(kml [mbarr] (deg. Kelvin] Imbarr] lg/m**Z] [%] 



Mldlotltude Summer atmosphere 






1013 


294 


24,87 


14.0 


76.3 


1 


902 


290 


19.38 


9.3 


64.2 


2 


802 


285 


14.03 


5.9 


55.3 


3 


710 


279 


9.35 


3.3 


45.4 


4 


628 


273 


6J1 


1.9 


39.1 


5 


556 


267 


5J1 


1.0 


31.5 


6 


487 


261 


2.44 


.61 


30.1 


7 


426 


255 


1.49 


.37 


29.2 


6 


372 


248 


.8! 


.21 


29.6 


9 


324 


242 


.645 


.12 


20.8 


10 


281 


235 


.234 


.064 


29.6 






flldlatUude Winter atmosphere 









1018 


272 


S.77 


3.5 


76.2 


1 


897 


269 


445 


2.5 


69.6 


2 


790 


265 


3.40 


1.8 


64.7 


3 


694 


262 


2.58 


1.2 


56.0 


4 


608 


256 


1.58 


.66 


49.2 


S 


531 


250 


.94 


.38 


46.4 


6 


463 


244 


m 


.21 


43.3 


7 


402 


238 


.31 


.085 


30.4 


8 


347 


232 


.17 


.035 


22.5 


9 


299 


226 


.087 


.016 


19.1 


10 


257 


220 


.044 


.0075 


17.4 



I 
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where K^^^ is the volume scattering coefficient for hydroscopic aerosols. For an exponentially decreasing aerosol 
number density of scale height H one can write: 



f = - 



|f(h)exp[-h/H]dh 



Jexp[-h/H] 



dh 



Applying this equation to the profiles of Table A.2.1 and assuming an aerosol scale height of 1 km yields f^ff / f(0) 

ratios of .95 and 1 .00 for the summer and winto' atmo^ho^s resi)ectively. A somewhat extreme scale height of 3 
km yields values of.85 and .93 respectively. 

These computations clearly indicate that the effective relative humidity may without serious error be replaced by the 
relative humidity commonly measured at ground level. A slightly more precise approach would be to employ a mom 
ratio value (i.e slightly less than unity). 
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AJ: SOME USEFUL APPROXIMATIONS 



The purpose of this section is to present some relevant approximations which help to simplify the physical 
interpretations presented in the main body of the text. These approximaUons were not used to create the tables and the 
figures in the main body but rather serve as a means of quantifying variations to the first order or better. 

(a) Mass scattering efTiciency for sulfates 

Figure A3.1 shows some calculations of the ratio of wet to dry geometric mean radius (computed from first 
principles without assuming a constant ratio r(0 / r(0), for a geometric standard deviation a ^ 1.24 and a density ratio 
of Po / P^ = 1 .8 g/cm3) as a function of relative humidity. These points can be fitted with a set of sectionaMy linear 
functions similar to those developed for the ratio of wet to dry mass in Appendix A.2. However for the level of 
approximation employed in this section a single monotonic relationship suffices. Accordingly an expression of the 
form: 

rN(0 = rN(0) [1 - '-^^ 

yields an accuracy of better than 6% relative to the points shown in Figure A.3.1 (relative humidities between and 
95%). ThU model is apphcable within the range of validity ascribed to the simple log normal results deduced in 
Appendix A,2 (average particle size » .Ol|im). 

Hansen and Travis (1974) defined an effective radius for scattering given by: 

J rjcr n(r)dr 



I 



J Jt r n(r) dr 



with an effective variance given by: 



^*ff - T 



1 



r 2 2 

J(r - T^) Kr n(r)dr 



■rf f 2 

I itr n(r)dr 



In the c^ of a log normal distribution the effective variance reduces to a particularily simple and familiar foim, 



viz: 



1 1 5,2 

Inr^ = lnrj^ + — Ino «• 

In order to obtain a simple and rapid approximation to ihe computation of scattering coefficients one can substitute 
the expression for the effective radius into the anomolous diffraction approximation (Deirmendjian, 1969): 

(^ = 2(1 + i.). i(sinp + ^^) 
pi p p 

2jtr ., . . ■ 



wtwre p = 2 {=-^ ) ( ra - 1 ) 



I 
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and where Q^^ (defined in Appendix A.l) is the unitless scattering effeciency. The range of validity for this 
approximation is nominally that the complex pan of the refractive index is zero or negligible, the real part of the 
refractive (rOj) is close to unity and that the size parameter a = (2 jc r / X) is » 1 . In practice it yields useful numbers 

beyond the nominal range of validity (in particular if a correction factor introduced by Delrmendjian (1969) is 
employed). For our purposes we employed the above expression without refinement Table A.3.1 is a comparison of 
scattering effeciency calculations using the anomolous diffraction approximation and exact Mie calculations applied 
to log normal size distributions. In both cases the relative humidity effects were computed using the algorithms 
defined in Appendix A.2. The table shows that the approximation becomes ccmsiderably less reliable as the 
wavelength increases into the iKar IR spectral region. 

The dry sulfate mass scattering effeciency is defined as (c.f. Section Z2.3 of the main text): 

Ak'"'(X) 



em = 



ASO4 



For a change in atmospheric state characterized by a simple increase in the total number of aerosols (c.f. continental 
standard or urban standard models of Table 2.1.4.3) the mass scattering effeciency for a three component atmosphere 
(sulfates, soots and dust like particles) becomes: 



sq, *"" ""' SO4 



e(X) = : + . + e(X) 



30; 



so. 



where the denominator of the first two terms is simply the mass of a single dry $04" particle and the Ys represent 
number density fractions (c.f. Appendix A.l). The last term can be written explicitly as: 



'^^ i«r3(0)TiD 

where T2 and rj arc area and volume radii defined for a log normal distribution in Appendix A.2, T| is the molecular 
weight fraction of sulfate to ammonium sulfate particles and D is the density of ammonium sulfate (1 ,8 g/cm^). 
Employing the expressions for die log normal radii deduced in Appendix A.1, the relative humidity approximation 
defined above and the expression for effective radius one finds: 



Q„.(i- f)"* 



A) - 



SO4 



so; 4 

3 



* =.TlDr^^(0) 



For fine particles defined by 1.5 < pS 4.5 the ratio Q / p is approximately constant ( ~ .8 using the anomolous 
diffraction approximation). Accordingly the sulfate scattering effeciency term above varies ^proximately as .8 p fAf) 

for small-particle size distributions. In this region then, we may characterize the scattering as being dominated by 
volume effects since the scattering cross section varies roughly as the cube of the radius. The mass scattering 
efTiciency expression becomes: 

2.4(m -l)(l-f)"'^ 



SQ 



'4 
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TabU A. 3. 1 

Comparison of oxact tcattarlng tfriclinu computetloni with valuai obtalnad uting 
tht AnomoTout DUfractlon Approximation (log normal ttza distribution) 

A 
Wevtisngth 0.4 0.iS 0.633 

[mlcrom.l 

m Q Q' m Q Q" m Q Q' 

1.530 2.066 1.908 1.530 

1.495 2.096 1.905 1.492 

1451 2.106 1.946 1.449 

1.430 2.246 2.005 1.427 

1.396 2.494 2.205 1.393 

1.371 2.086 2.562 1.360 



Re 


Bttve 


humidltg(«) rN 


raff 









.12 


.135 






50 


.129 


.145 






70 


.144 


.161 






80 


.154 


.173 






90 


.100 


,202 






95 


.245 


.234 



1.011 


1.120 


1.530 


.6033 


.906 


1.047 


1.125 


1.495 


.7103 


.904 


1.132 


1.155 


1.451 


.7946 


.920 


1.202 


1.196 


1 430 


BS48 


.964 


1.412 


1.345 


1.396 


1 036 


1.092 


1.773 


1.642 


1.371 


1.35 1 


1.353 



Notts 



m = rtal pari of Iha rafracllva Index (m(wat»r) s 1,339, U33,1.332 
for 0.4, 0,555, and 0,633 mlcrom. respectlvtiy) 

Q : txecl tcatltring effaclancies computed by applying ttw first principle 

particle growth equations , artilrsrily exact tlie calculations, and 
and subsequently Integrating Q(2"pi*r/1amda) over the entire 
size distribution 

Q' * scattering effeclency from the enomotous diffraction approximation 

A - log normal (sulfate) size distribution 

(rN(0) « .12 mlcrom., Sigma CO) = 1.24) 

rN : geometric mean radius, raff : effective radius Imlcrom.1 
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and hence is independent of particle sizing parameters. The other components of the total sulfate mass scattering 
efficiency (soot and dust) are sufficiently small at visible wavelengths that they do not significantly disturb this 
condition of approximate constancy. 

Actual values of the mass scattering effeciency (computed using the anomolous diffraction aj^noximation and the 
relative humidity algorithms of Appendix A.2) at a wavelength of 0.55 pm are graphed in Figure A.3.2 as a function 
of the dry effective radius. The dashed region for which approximate constancy of e appUes is indicated. The dry 
effective radius was allowed to vary by eithor changing Ac dry geometric mean radius ( .12 S rj^O) 5 .24 ) or the 
geometric mean standard deviation ( 1.2 < os 2.6) until the wet effective radius ( p^O ) exceeded the fine pjoiicle 
limits (indicated as dashed lines on the graph). 

(b) Single scattering albedo and the mixing ratio of soot to sulfates 

At a given altitude in the atmosphere the aerosol single scattering albedo is given by (c.f. Appendix A.l): 



^o 



Km 

where the subscripts C. SO4 and D refer respectively to carbonaceous soot, sulfate and dust like particles and the 
symbol "N" represents number density. If the contribution of the dust like particles is neglected along with 
^sorptivc component of the sulfates, the expression for o) reduces to: 



<w 


X + 1 


X • 


f 1 • 


Nc 


«fcw 



% = : , whee 



X s 



and where %^(k) is the single scattering albedo fw soot particles only. The relative humidity (f) dependence of the 
sulfate scattering cross section is shown explicitiy. This latter parameter along with the ratio of soot to sulfate 
scattering efficiencies is easily computed using Mie scattering computations. Accordingly there exists an approximate 
one to one relationship between Uie aerosol single scattering albedo and the soot to sulfate number density mixing 
ratio. 

This simplification can be exploited to extract the single scattering albedo given a value of the mixing ratio and 
relative humidity or conversely a value of the mixing ratio given an estimate of the single scaOering albedo and the 
relative humidity. Figure A.3.3 is a graph of the aerosol single scattering albedo as a function of the mixing ratio 
parameter "x" computed for the urban standard and continental standard aerosol models (see Chapter 2 in the main 
body of the test) at two different relative humidities and a wavelength of 0.633 \m. The dashed carve is a plot of the 
simplified expression presented directiy above. The difference between the two curves, which is not large, could be 
further improved by inotxlucing empiric^ correction factcvs. 

(c) Error rate approximations 

The error rate calculations performed in the main body of the text can be, at least in a first order sense, vastiy 
extended with respect to the variation of the input parameters. To this end we introduce the analytical/empirical 
function for apparent refiectance: 

p* - pTexpL-kjCSOj] + kjpeSO^ + (k,-k2)p^£S0j + p^ 
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where kj and k2 are constants to be detennined, p^ is the critical reflectance, p is the surface reflectance, e is the 
aerosol mass scattering (per unit mass of sulfates), p^, is the background apparent reflectance in the absence of 

sulfates and T is a transmission function which accounts for the effects of attenuation due to molecules and aerosol 
attenuators other than sulfates. This simplistic model is constructed around the assumptions (i) that the apparent 
reflectance is linear with surface albedo (ii) that the apparent reflectance is linear with sulfate concentration in the 
presence of a non reflecting surface (iii) that the addition of a reflecting surface adds an exponential attenuation term 
plus a term linear with sulfate concentration and (iv) that a critical reflectance exists where all sqiparent reflectance 
lines intersect independent of the sulfate concentration (c.f. Section 3.4.2 of the main text). The latter assumption is 
relaxed somewhat in the model since a unique intCTsection point exists only for small sulfate concentrations. A 
further simplification (which is not an intrinsic limitation of the model above) was to assume that the cridcal 
reflectance was independent of the scattering effeciency cr hence of the relative humidity. Finally, the expression 
includes the tacit assumption that changes in the apparent reflectance due to variations in the phase function are 
second order. 

Given die above analytical equation it is a simple matter to derive the three basic partial derivatives (error rates) 
introduced in Chapter 4. For simplicity we drop the standard partial derivative subscripts unless their inclusion is 
essential to an understanding of a particular equation. The first error rate is obtained by sought differentiation, viz; 



1/e 




•^iPc + l^(P-Pc) - kjPeTcxp(-kjeS04) 

This error rate shows most the principal properties of main text Figure 4.1.1 including the change of sip across the 
critical reflectance, and the reduction of the error slope with increasing optical activity ( increasing £). It does not 
respond to the abrupt non linear dependences cm sulfate concentration which are responsible for some of the more 
extreme excursions seen in the same Figure. Table A.3.2 below summarizes some of the comparisons made. 

Given the above expression and given that the denominator is not strongly dependent on e one can express the 
dependence of this error rate in terms of a reference enor rate according to the approximate relation: 



8S0j t^ /SSO^N 



For the surface reflectance error rale one can write (constant p*): 



. (Jel)^^o: - (^)„.^ 

5SO4 ^ 5p SO^ 

If the sulfate dif ferendal is due to a change in surface reflectance under conditions of constant apparent reflectance one 
can replace the differendals by their appropriate partial derivatives and after some rearrangement obtain: 

(f2^)...(!!°I)(^). .,„., 

6p P 5p* P 5p so; 

^^^so: = ^S0;k2 + Texp(-kj£S0;) 
Unless the sulfate concentration is very large the latter term is of the order of T. The surface reflectance error rate is 
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consequently strongly associated with the apparent reflectance enor rate with a change in sign, viz: 

^ Sp f^ ^ ^ gp* ''p 

and one may effectively describe its behavior in reference to the apparent reflectane error rate. This approximate 
correspondance can be verified by comparing the main text error rate curves in Figures 4. 1 . 1 and 4. 1 .3. 

The sulfate mass scattering cffeciency error rate reduces to the expression: 

5SO4 so" 

Se t 

This simple rel^on results because the analytical expression for p* defined above is constant if the product eS04= is 

constant The sign of this error rate is always negative b«;ause an increase in e implies a greater optical thickness 
which must be offset by a decrease in sulfatge concentration if the apparent jeflectance is to remain constant. 

If the change in e is brought about by a change in relative humidity then the actual value of the error rate (absolute 
value) will be less than that indicated above. This is because an increase in relative humidity implies a decrease in 
phase function at large scattering angles typical of satellite remote sensing (see the m^n text Figure 2.2.2.3 for 
example). Accordingly an increase in c due to an increase in relative humidity induces a decrease in the phase function 
and hence a reduced requirement for the offsetting sulfate reduction. Figure 4.12 in the text clearly ilustrates this 
effect 

In the case where the change in e is due to a change in relative humidity the mass scattering efficiency error rate is 
easily converted into a relative humidity error rate, viz: 

5S0j . 5S0j 
Sf ~ Sf 5e 

As noted above in the previous section, the mass scattering effecicncy can be expressed approximately in terms of its 
dry value as: 

e = mn-ff sothat; 

■*■ = and hence; 

Sf 1-f 

8S0r _ ^_n_ 

Sf ' " M-f 

where f is in the fiBcdonal representation. The exponent "n" varies with particular (dry) size region and generally 
increases with decreasing (dry) particle size. It can be computed from numerical values of e such as those presented in 
Figure A.3.1. 

The linear increase of mass scattering efficiency error rale magnitude with equivalent sulfate concentration is a 
phenomenom which can be observed (except for the abrupt non linear changes) in the main text Figure 4.1.2. If a 
given atmospheric state differs from a reference state only in that £ is different at the same relative humidity, the mass 
Kaocring efficiency eiror rate can accwdingly be related to the reference error rate by the expression: 



5S0: e^ /5S0! 



& e 6e 5 



Table A.3.2 Apparent reflectance and apparent reflectance enor rate (dS04- / r*) comparisons between 5S 

calculations and a simple analytical model (see text). Conditions are nadir geometry, solar zenith angle 
= 45°, 0.55 jjm wavelength and the urban variant model. Ihe selected model parameters were p = .2, p 
= .036, kj = 2.88E-4 and k2 = (the model was largely insensitve to the exact values of kj and k2 if 
k^ - k2 was kept ctMistant). The mass scattering efficiencies for 70%, 90% and 95% relative humidity 
wa« respectively 6.84, 13.0, and 23.4. 

(a) Apparent reflectance 
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(b) apparent reflectance error rate (ng/m^ 


per percentage increment in p*) 
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5S04= 


/p* 






[%1 


(5S) 


(analytical approximation) 



0.0 70 38 25.4 

0.0 SO 70 21 25.4 

0.0 50 90 10 13.4 

0.0 SO 95 $ 7.41 

0.5 70 -5 -23.7 

0.5 SO 70 -m -28.9 

0.5 SO 90 -37 -18.6 

0.5 SO 95 <-100 -153 
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APPENDIX B 

Comparison of Radiative Transfer Models Used to Determine Atmosplieric Optical 

Parameters from Space 
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Comparison of Radiative Transfer Models used to Detemiine Atmospheric Optical Parameters 

from Space 

Alain ROYER, Norman T. O'NEIUL, Anthony DAVIS * and Uurent HUBERT 

CARTEL. Universiai dc Shcrbnx>ke, Sheibrooke, Qudbec, Canada JIK 2R1 
* Department of Physics, McGill University, Montreal, Quebec, Canada, H3A 2T8 

ABSTRACT 

An array of models and techniques e^isi for the calculation of the aimospherk: backscatteied radiance and ground 
reflected radiance received at satellite altitudes. In remote sensing applications, where one deals in mep-byte units of 
data, it is essential that these models be computauonally fast while retaining a reasonable degree of accuracy. We 
have evaluated a number of such models (Lowtran 6, Turner, Discrete ordinates method, 5S) relative to an accurate 
multiple scattering, multi-layer (Dave) model in order to assess the performance of these models in an inversion 
scheme f<^ aerosol optical depth. 

The Turner and single scatter Lowtran 6 models generally produced large errors in apparent reflectance. Overall, the 
Turner model was not significantly better than the Lowtran 6 model except at near nadir geometries and non zero 
albedos. The 5S model which is orders of magnitude more rapid than the DOM model was significantly more 
accurate than the L6 and Turner models. The accuracy of the inversion procedures for the extraction of aerosol optical 
depth from satellite apparent reflectance was then analyzed for the two most precise models (DOM and 5S). Fat 
typical measurement conditions, the 5S inversion errors were found to be of the order of .1 in a turbid atmosphere 
case (aerosol optical depth approximately .5). The DOM poduced the most impressive results in terms of 
comparisons with the Dave computations. It's time of execution, however, is a serious constraint with respect to 
saieUiie remote sensing applications. 

INTRODUCTION 

Arbitrarily exact numerical methods for the solution of the radiative transfer equation in planetary atmospheres 
have evolved in parallel with the processing power of digital computers (see for example Hansen and Travis) . In 
recent years there has been a natural and beneficial tendency to rationalize the cumbersome numerical overhead 
associated with generalized solutions by extracting the essential physical characteristics governing a particular 
problem in order. to achieve solution accuracies commensurate with measurement precision. In terrestrial remote 
sensing applications where one deals in megabyte units of data the rationalization approach is a necessity bom of a 
need to achieve practical data processing throughput rates. This problem is further complicated by the existance of 
spatial inhomogenieties in the atmosphere which, short of a full scale 3 dimensional radiative transfer solution, 
require at least a mosaicing of individual solutions. Finally, if the intent is to invert remote sensing data for 
purposes of studying the atmosphere itself, the tractibility of the inversion methodology is clearly dependent on the 
speed and simplicity of the radiative transfer solution. 

The objective of this paper is to evaluate a number of commonly used radiative transfer models in terms of relative 
precision and speed and to investigate the implication of these constraints on inversion algorithms for the extfaction 
of aerosol optical parameters from passive remote sensing imagery. Aerosols arc here defined to be particles whose 
dimension is the order of a wavelength and whose optical variation is largely confined to the troposphere. The 
wavelength region of interest extends from the near UV (.4 ^un) to the near IR (1 .6 \m). 

The a^roximate radiative transfer models which simulate radiance signals received at satellite altitude r^resent 
various levels of sophistication ranging from single scatter models to numerical single layer solutions of the 
radiative transfer equation. The strategy employed in this study was to analyze these models relative to a multilayer 

multiple scauering numerical solution (Dave) . Our goal was thus to clarify the uncertainty associated with 
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modeling contributions to the total eiror budget and hence facilitate tradeoff decisions of speed vctsus accuracy in 
model selection. 

In the first section of the paper we present a review of past radiative transfer work in the context of remote sensing 
and inversion methodologies. The second section deals with the relevent radiative transfer concepts followed by a 
description and physical interpretation of the radiative transfer models which we selected for comparison. These 
models id^tified as L6 (Lowtran 6 single scattering model (Kneizis ct al.) .the Turner modd (Turner and Spencer) . 
5S ("Simulation of the Satellite Signal in the Solar Spectrum". Tanr< ct al. ) and DOM (discrete ordinatcs method, 
Liou', O'Neill ) arc eithor readily available as computer codes or are well described in the literature. In the third 
section the models are compared with the baseline DAVE results in terms of apparoit reflectance difl^erenc^ as a 
function of a variety of geomcirkal and environmental parameters. The siibse<iuent section reviews the impUcation (rf 
these differences in terms of inversion procedures for optical depth. We conclude by presenting a strategy for the 
exploitation of the models investigated. 

1. BACKGROUND 

For remote sensing applications the driving force behind the development of radiative transfer models was a 
requirement for atmospheric corrections to remove what was considered to be a contaminating influence on the signal 
reflected from the tencstrial surface (see rcf. Bariou ct al. for a comprehensive survey of remote sensing inspired 
radiative transfer models). In as much as the estimation of atmospherically scattered signals (path radiance) is critical 
whether one speaks of atmospheric correction or inversion, the modek developed in the former application have 
proven valuable in the latter. Note that, within the context of this present study, we arc excluding from 
consideration those scmi-cmpirical or ratio types of approaches which extract values of the path radiance from other 

radiometric measurements (see ref. Gordon et al. for example). 

The first step towards a pragmatic radiative transfer algorithm was the development of the Turner model (Turner 
and Spencer)' . This model which was intended to be an improvement over single scattering models was appealing 
in its analytical simplicity and rapidity of execution. An advancement over this type of model was realized by 
Gordon ' and Tanr6 et aL* who revived the concept first introduced by Sckara of effectively separating the 
molecular and non molecular radiance fields. Further improvements included a capability for modeling the effects of 
spatially inhomgeneous surface reflectance (Tanrd et al.) , more realistic 2 layer modeling of vertical 
inhomogeneities (Deschamps et al.) and a simple method to account for errors associated with the decoupling 
assumption of the Raykigh and aerosol radiance fields (ibid). 

The inversion of radiative transfer models to extract inherent optical parametns is a process susceptible to 
modeling errors and radiometric c^bration errors. Modeling etmt% can be further classified into two types: numerical 
modeling crrore which result from approximations applied to shotcut the computational work necessary to achieve 
an arbitrarily exact solution of the radiative transfer equation and physical modeling errors which result from 
limitations in the characterization of the optical parameters input to the radiative transfer equation. The first type is 
the subject of this paper. 

The physical modeling errors are coupled to the amount of apriori knowledge and degree of physical reality 
incorporated into a given input parameter to the radiative transfer inodel. In the context of an inversion procedure the 
parameters to be inverted (which are normally input parameters to the radiative transfer equation) are forced to 
converge to values which reflect the reality of their individual constraints and the nominal values assigned to the 
other input parameters. Depending on the type of inversion which is being performed the physical modeling errors 

can include: phase function variations related to changes in size distribution (Manellacci et al.) , variations in the 
real and complex parts of the refractive index (Fraser and Kaufman ; Hand ), the neglect of polarization effects 
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(Hansen and Travis ; Kattawar ct al. ' ^). departures front Mic scattering phase functions (Schuerman)' * . and the 
effects of atmosphOTc stratificaiion (Liou) as well as horizontal inhomogencity (Mishin)^ ^ . 

Other sources of physical modeling errors are related to the influence of thin cirrus clouds (Saunders and 
Kriebel ). the effects of spatial non homogenieties in surface reflectance (Tanr6 et al. ^ ' ; Mekler and Kaufman ^ ^) 
and departures from the Lambertian assumpUon (Tanr6 et al.)". Gaseous absorpUon and the way in which it is 
incorporated in the modeUng may also play a significant role. Real perturbations from nominal assumed values for 
absorpuon due to amongst others, water vapour, nitrogen dioxide, and ozone (visible spectral region) may dependine 
on the precision required and the wavelength range selected, bias the radiative transfer invexsion algorithm. 

Errors of ra^omctric calibration have increased in importance with the evolution of modeling and sensor 
technologies. The fust gcnerauon of broad band sensors (GOES-VISSR. NOAA-AVHRR. and NUMBUS-7-CZCS) 
produced useful results which however suffered from the absence of post launch calibrations (Fraser et al.^\ Frouin 
ct Gauthier ; Gordon ) and limited radiometric resolution (Griggs) ^ ^ . 

The aerosol optical parameters which are most readily extracted by inversion techniques include the aerosol 
scattenng optical depth used as an mdicator of total columnar aerosol concentration, the single scattering albedo as a 
measure of absorption (Kaufman) and opucal depth spectral coefficients as indicators of the average atmospheric 
size distribution (Aranuvachapun)*'. The principal influence on radiance variation and hence the most easily 
determined opucal parameter is the aerosol scattering opucal depth. Inversion algorithms for the extracUon of this 
parameter are further simphfied by it's strong linearity versus radiance over large ranges. This relaUonship is well 
known in modeling studies and has-been observed in empirical comparisons between satelUte data and ground based 
extinction measuremaits (Griggs ;Meklcretal. °; Norton el al.). 

2 . THE RELEVENT RADIATIVF. TRansffr TH E ORY AND model DF.srRTPTTni^ 
2.1 Statement and Red uction of the prohlerr^ 

dii?.,!!in^«t."h! ^"^^ °'''"!fT °^'^*^^« ""^f"' ^^^ ^ presented i" ordtt to facilitate and complement the 
discussions of the various models which wci« studied. More detailed and comprehensive reviews can be found in 

fundamental texts such as Chandrasekhar , Lenoble and Hansen and Travis^ ' 

TTie radiation field in and around the eartlis atmosphere is best described by the concept of radiance I(r.t;flA); the 
amount of EM flux (energy/surface /time) propagating into direction Q (unit vector) at wavelength X per unit of solid 
angle per unit of wavelength, measured at position r and time L I obeys the following kinetic equation in a medium 
of constant refractive index: 

1 dl 

-5j.a.Vl--a(I-S) ^ 

where C is the speed of light. V the gradient operator, and a the volume extinction coefficient (i.e. absorbtion plus 
scattering). The source function S is given by: 

S(Lt;flA) - ^ ^p(aa')i(i.tii'.x)da' (2) 

where the integral symbol represents an integration in solid angle over all directions (4tc steradians). 
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•mis describes coherent scalu^ring processes (i.e. excluding redistribiiion in wavelength). TTie single-scaaering 
albedo tUo is defined as the ratio of the scattering (a) to exiincUon (a) coefficients. TTie phase function p depends on 
a»fl'- cos(scattcring angle). We adopt the usual normalization for p, viz: 

jpiOmdO: ' An '■ ^^ 

To determine I completely, boundary conditions must be spccmed, these boundary condiUons take the form o. 
inimtng Z^ce^ ^1 interfaces ofSte medium with its environment In ti^a^of phnetary a^spf^^^J; 
i^S c^ditions are specified by an external souixx of (coUimated) radiauon bom above and. often, a pamaUy 
reflective surface below. 

TUc general three-dimensional Radiauve Transfer (RT> problem is extremely difficult to solve although progress is 
being made towards a general medKxl (Stevens)". Nevertheless many solutions have been obtained m pmicuUr 
cases (one of which is mentioned further in connection with the 5S model). TTic time-dependent problem has been 
applied to lidar signal analysis. 

In the rest of this study (and most RT research to date), stationarity and horizontal homogeneity are assumed. For 
geometry which is plane (wrallel the RT equatirai becomes; 

H^.-a(I-S) M 

GZ 

where tt is the direction cosine relaUve to the vertical axis (Qz) and z is the vertical position coordinate. I and S 
are now functions of z. n. <|> (azimuthal angle) and X (wavelength). The X dependency enters through the opucal 
parameters a. a, tUo and p which may as weU be functions of Z. In this context, one can define the dimensionless 
Space coordinate known as optical depth: 



t" 

1 



Ja(z')dz' ^ 



Hence the standard RT integro-dif ferential equation: d^ " ^ ' 47 J ^^^^' 

with boundary condiUons KOji.*) - E5 6(H+Hs) m. H<0 I(t>,<>) - ^ F.(t*). i^>^ C?) 

2n ±1 
where F±(t) - I |l(t»i,^)jidMd« ® 

, . 

designate up/downwcUing irradiances. Es is the extra-terrestrial irradiance incident at zenith angle Gs- cos Hs »> 
the plane ♦ = 0. f is the optical thickness of the atmosphere (z: 0-) and p is the reflectance of a Umbcruan 
surface kxated at t-xV For simplicity we have dropped the explicit dependence of I on X. 

Notice that (6) admits the integrating factor exp{-T/n). Accordingly the formal solution to the RT equation 
yields the (total) radiance: 
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2.2 Varifiiis annmarhe.s for -living the RT pr^^ f m 

sorrrllT °T^ "'"'" ^'^ "^ '^'^ '" "'^ ^""^"^ ^'"" '^"^ i;-" <^'ffi-'»y "- - evaluating the 
source function S. A number of approaches arc reviewed in L«nnhl..'^ ch ^ • • ""^ ""^ 

00. X and p. The phase funcuon .s often expressed as the Ugendre polynomial summation: ^" 

WC0S9) . 2ftPl(COS9) (,jj 

Rayleigh scaoolng by panicles much smaller to U,e wavclenph of incidem mdiadon is dcscHbcd by ,™ ie^s (B„ 
n"i£of2r° ■"— '^-P'-' f-uonscha^^rislc of aunospteic ..sols ,«,ui. a'^Mv^lj^ 

coSrsfa?:S'o"^s'yt^L;rotch':;s^"=r"75! "-^ ■^^ «- -""™ 'f '-gh 

expression being aken for Sc^anSim^H^"^ ^i' °"'«°^ *"''='' Wically involve «, approximaie 

milels iha, udlL an aZoxL^^'?"^ „h ' ?? "T^ 'TT " «''^'»' ®- TV second .^ to„ 

app^ch by ^cSSg d» ugtf^n':;s5:rdKX fuTcSS^r iT ■" "^^ '''' '"■""'^ •^■' 

2. 3 Applicai i nn in P>moir Sm-in. . r,^„..,ri, ^ | „ ^ R,Hin,„„ri, r-pp ,.fn||„„. 

a." ^ olTtf l!;T«lTlrn7,^ "> P-iic, Uje pad, radiate and hence d« „u. »iianc above ^ 

,-«n^c„i.deJS;s«rnS:^Sfri3X"hic;:^aSr^tr *^ "-"^ -^ ^"'™ 

sa^^L-siSl^rS^lijcet ra^rS^-cllTneJbr "'" '""' ^"^"'"" "■'^' "" "^ 

.. ill 

•^ "usEs (11) 

bonr^TonSs 7^':, *: :T'L7T"' "^^ ''^ '""^•^' « » f--^- or «» wavelengd,, d. 
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variables of RT theory were assigned the values ev-0»(10°)70° and (t>-0'(30°)180» 
space were computed per atmospheric state and per model ( see Table 3-1). 



In all, 8400 points in parameter 



v 



Figure 2.1 SchemaUc diagram of satellite mcasuremem geometry ; Ss is the sun's zenith anglc.ev is the 
observer's zenith angle and <b is the obscrvCT's azimuth angle relative to the solar plane. 

2.4 DescriPt inn of radiative transfer models 

In tciTCStrial remote sensing applications the total optical depth (t*) employed in a given model is composed of 
scattering contributions from Raylcigh or molecular particles (Tr), absorbing and scattering contributions from 
aerosols (Xa). and absorption contributions from particular molecular constituents (Ztgo). The models employed in 
the comparative study were either multi-layer (Dave and L6) or single layer (Tunicr. DOM, and 5S). The single Uyer 
models ignore vertical variations in the phase function and assume a unique representation which is a weighted mean 
of the Raylcigh and aerosol phase functions where the weights are proportional to the total Raylcigh optical depth 
and the aerosol (scattering) optical depth respectively (see Turner for example). The multi-layer models admii 
vertical variations in the phase function which is taken as a Rayleigh/aerosol mean weighted by the volume 
scattering coefficients characterizing each layer. The justification for employing the single layer approach is that, tt 
least at the boundaries of the atmo^ere. Uie radiance contributions arriving from all points along the line of sight 
are effectively represented by an atmosphericaUy avaaged phase function. Aerosol and Raylcigh optical depths which 
are the critical optical parameters in detcnnining atmospheric scattering contributitms must however be idenucal for 
all models. 

Dave Model 

Dave' describes an extensive dau base of atmospheric radiance calculations produced using the spherical harmonics 
approximation. We shall briefly review die main ideas, some of which are used in the development of the other 
models discussed below. First of all, I is decomposed into its "direct" and "diffuse" components. The solution for the 
diffuse component is obtained by summing the separate contributions of a radiauon field rcsulung from only a 
collimated (solar) source only (input conditions defined by (7) with p =0) and a radiation field whose only source of 
iUumination is an underiying Lambcrtian surface of reflectance p (boundary conditions (7) with Eg = 0). The source 
functicHi for die collimated source is given by: 
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S(ilD - ^ [ ^p{Q.a')I(iii')dQ' + p(fl^)Ese- 1. ] 



m 



which accounts for primary scattering of sunlight as an internal source. Invoking the linearity of (6). the two 
solutions are combined to yield I anywhere in the atmosphere from x = to t = t*. 

TTie method of solution begins with a (truncated) Fourier series in of both the radiance and phase function p{ n. 4,; 
H'. n Each coefficient of the radiance Fourier series is then solved for independently by expanding the corresponding 
phase function coefficients into an associated Legcndre polynomial series in ji and n' . The derived RT expression for 
each radiance Founer coefficient is then integrated over n to yield radiance coefficient moments. TTie result is a system 
of simultaneous Ist-ordcr ordinary differential equations in terms of the radiance cocOlciem moments. 

J^i^li^wIJi *r "'*^ "'.'^'l?? ^^.fr'^ '■""''^°" ^°' substitution into an azimuth independent version of 
equation (9). Finally, having solved for all the azimuth coefficients the Fourier series is computed to yield the 
radiance at a given value of «. TTiis is done for all levels in an atmosphere which may been subdivided into any 
number ot layers to account for optical stratification in the vertical direction. 

Models to he tg«;^f^ 

f^l°S^^^' **^*' *^ '""^'^ "''^ '*'* ^^^ '"'^'' '^ ^^'^ ^^^ '^«'"« 2- 1 summarizes the salient 

Tibkll: Modc]f( 
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TTiis model was proposed Originally by Oiandrasekhar". applied by Liou' and refined by Stamnes'* Likethe 

SfiriiT''"'" ^PP^^'"'^ ^" *«=«^ above the method employs a (truncated) Fourier expansion of I in <|>. 
One proceeds m an identical fashion as above but stops short of the integration of each azimuth independent RT 
?r.cZ ""^H ^ "^ "^^ ^"^"^ component of the azimuth independent RT equation is replaced by a 

Gaussian quadrature approximauon to eventually yield an algebraic eigenvalue problem. *> r 

The number of Gaussian pivot points used in the quadrature is generally referred to as the number of streams and 
defines the precision of a given DOM^a^lication. The solution of the eigenvalue problem (wh^expSSSn 
matrix symmetnes noted by Stamnes ) yields a general soiuUon to which the boundary value conditions (7) are 
applied. The resulting linear system of simultaneous equations yields values of the radiance azimuth coefficients at 
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discrete (Gaussian) values of |jl These discrete cocfncients are then substituted into the source functon and a second 
quadrature performed to yield the source function at any desired value of ji (Stamnes) . One then proceeds as 
outlined above to evaluate the azimuth coefficients and Fourier series. 

1 9 
Although the DOM can be formulated as a multi-layer model (Liou) the version used in this study was single 

layered. The number of streams employed in the computations was 32 (see Section 3). The effect of gaseous 
absorption was folded in by multiplying p* by exp[- (— + -) Ztgu]. 

(ii) Turner 

This model, formulated by Turner and Spencer , yields a closed-form expression for I*. The authors derive a 
simple analytical expression for the source function S by substituting in an angularily primitive expression fot the 
integrand radiance I in equation (2). The first term of this expression is the sum of a double delta-function aligned 
with the solar and anti-solar directions of the solar beam and weighted by F±(t) respectively (equation 8). The second 
term is an isotropic expression which attempts to account for the source function contributions of surface-reflected 
iiradiance. The relevance of this intermediate radiance expression is thai its integral over solid angle yields fairly 
reliable two flow estimations of the downwelling and upwelling irradiances at every level in the atmosi^ere. 

The Turner model was originally developed for conservative scatiering (tOo=0' Although the code could be easily 
extended to accomodate non conservative scattering, we chose to retain the model in its original form. Gaseous 
absorption is included following the same recipe employed for the DOM. 

(Ui) L6 

This is the 6* version of a series of codes produced by the AFGL (Kneizis el al.) . The primary purpose of this 

B 
model is to compute atmospheric transmittSKe [cxpt- fa(s)dsl between two arbitrary points A and B. This is done 

for a variety of atmoqjhcric components including ozone, N2. HNO3 and wata^ vapour at moderate spectral resolution 
(20 cm-l). from the thermal K (350 cm'^ ■ 28.6 ^m) to the near UV (40000 cm'^* .25 urn). At IR wavelengths the 
model utilizes an assumption of local ihemodynamic equilibrium to compuw the contribution of thermal emission. 

For our purposes the most relevant feature of L6 is its incorporation of a single scatter code to permit calculations 
of cxo-atmospheric radiance (or hence apparent reflectance). The single scatter contribution is obtained by 
substituting the second torn of equation (12) into equation (9). 

The L6 code as utilized for our calculations accounted for only the reflected component of direct solar radiaticxi. 
Accordingly the contribution of downwelling diffuse radiation to the reflected term of equation (9) plus the 
contribution of reflected and then scattered radiation to the second term of (9) is ignored. As well the code was 
modified to contct for a fk;tor of 2 differeiices between the L6 surface reflected radiance and the radiance expected from 

a Lambertian (constant reflectance) surface. This discrepency was also noted by Richter . 

Other features of L6 include it's incorporation of spherical refractive geometry and the treatment of vertical 
stratification in terms of a multi-layer atmosphere. For the range of geomcffies employed in this study, the exclusion 
or inclusion of spherical refractive geometry does not play a critical role in the evaluation of the single-scatter 
dominated radiance computed with L6. The multi-layer aspect complicates comparisons made with L6. In order to 
impose user-defined a and p values above the boundary-layer (0-2 km), it was necessary to modify the driver program. 
As well, care was required to ensure that the optical depths utilized in our comparisons were compatible with the L6 
optical parameters. 
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iv)5S 

I, 
This model (Tanrd el al.) departs from the previous approaches since ihe authors literally build, picce-by-piece. an 
analytical expression for the total apparent reflectance: 

P'( P.e, ;ev.4> ) - patmCe, :9v,«) + p-^^sHIM (,3^ 

The parameters of this expression are: 

- Intrinsic atmospheric reflectance (patm): This term is the sum of separately derived formulae for Rayleigh and 
aerosol apparent reflectances. TTic aerosol contribution is obtained using a Sobolev approximation combined 
with a single scatter model. Aside from geometrical dependencies the form derived is an explicit function of the 
aerosol parametCTS ta. the phase function p, and the asymmetry factor (g - < cos X > = Pi / 3. X == scattering 
angle). 

- Total atmospheric transmittance ( T(6) ); This factor is the sum of a direct transmission term plus a diffuse 
radiation transmission term. The latter which is obtained from a delta-Eddingion approximation is an explict 
function of t* and g. 

- Spherical albedo (a) is a function oft* alone and is obtained using a semi-<mpirical expression based on single- 
scattering and Monte-Carlo simulations. 

Like DOM and Turner, 5S is a single-layer model, gaseous absorption is also treated in an independent fashion 
by the application of a transmission factor to equation (14). The 5S code includes its own routine for computing 

gaseous transmittance . One noteworthy extension of 5S, discussed elsewhere (Davis and Roycr)*'. is the 

accomodation of a simple but effective model for simulating inhomogeneous ground conditions. 

3. COMPARISON OF MOnFT S 

The four models were compared with the Dave model over a variety of input geometrical parameters and four 
different atmospheres (Table 3. 1). These atmospheres range from a pure molecular or Rayleigh atmosphere (identified 
as aerosol type 1) to a fairly turbid atmosphere with a total (absorption plus scattering) aerosol optical depth at .5 
lun of about 0.5 (aerosol type 4). The aerosol type 4 atmosphere is characterized by weak aerosol and molecular 
absorption and a Deiimendjian Haze L (continental) size distribution. The aerosol type 3 atmosphere is identical 
to aerosol type 4 except that its aerosol abundance is only one fifth of the type 4 value. Aerosol type 5 is 
chaiacterizcd by an optical thickness which is less than that employed in aerosol type 3 and a Deirmcndjian Haze M 
(maritime) size distribution. The phase functions utilized in the study were computed for each aerosol type and each 
of the S input wavelengths using citho- Mie scattering functions (Royer ct al.) or Lcgendre coefficients (O'Neill) * . 

l£ and 5S are delivered with a selection of standard molecular and aerosol atmospheres. To establish a significant 
comparison basis, all models were run for identical atmospheres. This necessitated the development of a common 
input data base of c^tical parameters (Table 3. 1) which had to be adapted to fit the input requirements of each model. 
Moreover. 5S, and to some extent, L6 feature a variety of user-friendly interpolation and integration schemes as well 
as techniques for calculating geometry, aerosol mixtures, ground surface reflectances from standard reflectance datasets 
and conversions from visibility to aerosol optical thickness. These features were bypassed or eliminated or adapted 
for our study. 

Data bases of total apparent reflectance ( p* ) were computed and subtracted from the DAVE data base to yield an 
ensemble of difference files for each of the four models as a function of four different aerosol types. Representative 
plots of these differences are displayed in Figures 3.1 to 3.5 respectively as a function of the geometrical parameters 



126 / SPIB Vol. 928 Modeling of the Atmosphere 119881 



TYPE OF ATMOSPHERE 



MM lalUlsU MBMr wtiT TSTMr nataat; IM (.«>1; 

VwUul pronta: I( lw*to '« pmNrt, Itapwilan. MM* u4 « 



■ter Ttpw KroB • I* M fea. 



TYPE OF AEROSOL 



K*(ncllvt >■<« 



I Tkta Hu« L 

4 Tklck Ban t 

5 Hiu M 



1 J - I •.•! 
IJ - I •.•! 
IJ - I •.•! 



a . I, b ■ If-lIM, T ••■S 
a ■ 1, k ■ 1(.11M, T -••3 
a > 1. k ■ l.»44], T ■•■< 



TtM Aanwl • fcaiHT 

If.lIS 
4.(T3 



Vtrtksl proni*: uiM U Imli far a(r«*l« t dauNT 









RELATED PARAMETERS 






YYiTtlcn*"' '"■'^- 


.US. 


Aaroaal SIb|I> Sc«IKrta| Alb*4* 


Mil 


LH. 








Aarowl 


Typai 












3.4 

J 




.I»T11 
.134*1 


.fll«4 
.15*1* 

Aaraaal JUjmm 


.tl311 
.17114 

alrr Factan 


.*34«l 
.1*1** 


.•41*1 
.)344i 


1.4 
5 




.TIT»1 
.TS2*( 


.71173 
.737*4 

Aaraaal ScattertH 


.71*14 
.71(47 

OpUtal DapUu 


.7*417 
.71711 


.*Sli5 
.7***3 


1 

i 




.ilflt 
.44«t* 
.•4I1« 


.t>II* 
.41S44 
.•41*1 


.•f«4* 
.4SI«I 
.*44«* 


.*I4S* 

.41111 
.•4*11 


.•4I*T 
.14331 

. 14111 








Aaraaal AkaarpUea Optical DapUi 






3 




.•1«1J 

.•1113 

.••its 


.**(T1 
.•435» 
.**711 


.•*7S3 

.•37*4 
.•**4t 


.*«Sil 

.•1*41 
.••13* 


.**2tl 

.*14t7 
.••I*! 



1 (tajlal|k) 3,4,3 .1132* 



l„3.4,S 



l,J,4,f 



1,3,4,5 



.**IST 



,••••• 



.••*•* 



KayWik ScaUarlac Optical Paplki 
.*f44( .MII3 

Otaw OpUcal Daptka 
.*iT71 .*It4t 

Water Vapaw OpOcal Daptka 

.•*••* .••••• 

Carkaa MasMa Optfcal Paptk i 
.*••*• .*••*« 

Tatal Optical Itaptki 



%TM 



.***•• 



.*•»> 



.***•• 



.••U9 



.••••• 



.*••*• 



.•14** 



11477 


.1111* 


.•((71 


.till* 


.(1(11 


3144* 


.111*1 


.1*4(3 


.113*3 


.*«77l 


7ll»» 


.(1111 


.SS«3I 


.4747* 


.17311 


1(171 


.17141 


.1177* 


.•741* 


.•*t44 



IMPUT PARAMCTEB RANGES FOR RADIATIVE TRANSFER CALCULATIONS 



Vlaalat aaila *, 1*, », 1*. 4*. 1*, *•, 7* 

Ailaalh aaila *, 3«. **, **. U«, 15*. 1(* 

Solar taallk aa|la 34, 45, (*, 7*, 75 

Wa*el*Btlk 
Croaad alk«4o 



.455*. .5»*, .(55*, .1455, 1.(1 

*. .2, A, .(, J 



{ngmt* 1.1) 

(ntara 1.1) 

(ngara 3J) 

(ngara 3.4) 

(ncara 3.S) 



Tabic 3.1 SttBaar; oC lk« apUca-pkjilcal alBsapfccrlc paranlart aa* lapirt paraaalarc 
r«r Ika radlaUn traaafar cMpariaan. All aM>M ir* ■■ ^afraaa. 



SPIE Vol. 928 Modeling of the Atmosphere (1988} / 127 



(Ov. <{> and 65), the wavelength and the surface renectanoe. We describe below the general behavior of the reHectancc 
differences for each model. 

With certain exceptions the errors tended to increase with decreasing wavelength or effectively increasing total 
optical depth (Figure 3.4), increasing slant range of the observer's line of sight or solar beam (Figures 3 13 3) and 
increasing aerosol opucal depth (i.e. as a funcuon of model type). Except for the L6 model (and the Tumcr'model to 
a less severe exient) the errors were reasonably independent of surface albedo ( Figure 3.5 ). The errors in the figures 
and the rms errors reported below can be interpreted in a relaUve sense by noting that the apparent renectance in a 
zero ground reflectance atmosphere is of the order of .1 for short wavelengths (.5pn; sec Figure 4.1a for example) 
and .05 for near IR wavelengths (.85 jim). The perspective which is most pertinent to the present study is how the 
Ap* errors translate into errors in aerosol optical depth obtained by inversion (Section 4 below). 

1. Lfi_m2d£l • Not unexpectedly, the single scattering compulations of the L6 model produced large apparent 
reflectance errors relative to the Dave model. Of note in the zero albedo figures, where backscattcring is important 
and the optical line of sight is large relative to nadir ((t>=0 in the zero albedo case of Figure 3.1a for example) is the 
predominance of single scattering over multiple scattering. This effect persists for aerosol type 3 and aerosol type 5 
unul the optical depth is sufficienUy large (aerosol type 4) that multiple scattering masks the low order scaaering 
influence. The large negauve errors observed for the non zero albedo cases ( e.g. - -.25 in the aerosol 4 case of 
Figure 3.1a ) arc related to the fact that the single scatter model accounts only for the direcUy reflected term of the 
total apparent radiance and as well underesUmates the direct term by ignoring the downwelling diffuse irradiancc at 
ground level. 

2. Turner Moric I - Our computaUons showed that the Turner model was not substanUally better than die single 
scaaer model except in near nadir configurations and conditions of non-zero surface albedo. Apparent reflectance 
errors were largest m the neighbourhood of the region corresponding to solar beam backscatter (Figures 3.1a and 3.2 
: <t>=0. aerosol type 4). This effect which is particularily evident in the more turbid atmospheres (i e the 
atmospheres characterized by more anisotropic phase funcUons) is an unfortunate artifact of the delta function 
representauon taken for the source funcuon radiance (see Section 2 above). Rms averages performed on ihe figures 
presented in this section indicated diat for the Rayleigh case A p* 5 .03 at X = .555 ^m and < .006 at X = .g455 
fim. For aerosol type 4 the anologous parameters were respectively A p* < .09 and .08. 

3. 3S Model - This model produced the most satisfactory results given the simplicity of its formulation The 
greatest errors occurred at large values of the observers or solar beam slant range and at large optical depth This 
level of dispanty is not unknown to the authors of 5S and can be improved upon by means of a relatively simple 
modification to the code which accounts for the coupling of the Rayleigh and aerosol radiance fields (Deschamps et 
al.) . Rms errors of apparent reflectance in the Rayleigh atmosphere case amounted to A p* ^ .003 at both .555 
and .8455 jun wavelengths. For aerosol type 4 the maximum errors were respectively .02 and .03 for the .555 and 
.8455 (im wavelengths. For n«ar nadir configurations these errors decreased by about one ordCT of magnitude. 

4. t >QM Model - Although not a computationally simple model in the sense of 5S or Turner, the DOM model is 
significantly faster than the DAVE code given that computations are performed for a single layer only It was thus 
of some mi^est to investigate the level of error obtained for an arbitrarily accurate single layer model relative to the 
mulu layer DAVE model. e, j «.i. vnuuic 

To test the convergence of the DOM model computations for the standard number of sti-eams employed (32) 
compansons were made with the Dave calculations at 16, 24. and 48 streams. Above 24 streams no significant 
improvement was noted in the relative precision. For the majority of nadir cases tested (with varying solar zenith 
angle and albedo m the aerosol type 4 atmosphere) the small reflectance errors in going from 32 to 48 streams 
ctenged by less than 50%. This suggests than these errors are primarily due to mulii-layer versus single layer 
differences between the Dave and DOM models. r / j g c wye. 

For the Rayleigh atmosphere the difference between DOM and DAVE results was < .0002. Wc note however 
that the neglect of polarization in both models can lead to extreme errors < 23% relative to a complete solution in 
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aerosol type 4 (AERO 4) atmospheres as a function of the cosine of the obterver's zenith angle for a 

wavelength of .555 urn and a solar zenith angle of 45 degrees. * = represents a radiance ray in the 
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Apparent total renecunce differences 
(relative to the Dave model) computed for 
the Rayleigh and aerosol type 4 (AERO 
4) atmospheres as a function of the 
observer's azimuth angle for a 
wavelength of .555 urn. The observer's 
zenith angle is 50 degrees, the solar zenith 
angle 45 degrees, and the surface albedo 
zoo. 
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Apparent total reflectance differences 
relative lo the Dave model) computed for 
the Rayleigh and aerosol type 4 (AERO 
4) atmospheres as a function of the 
cosine of the solar zenith angle for a 
wavelength of .555 urn. The geometry 
was nadir viewing and the surface albedo 
was taken to be zero. 



terms of Stokes vectors for a Rayleigh aunosphere (Katuwar ct al )' ^ F. 

computed for the figures in dtis section were C^5 J^. ^i L ' """'"'' '^ ' ^* "^ "^™«"" 

polarization in this latter case are con^i^S; fL <^L (1^' ^ ''' "' '''' ^ "^^^^y- -^ ^^^ of 

4. nWFR<:fn| ^ANATY , S r'T 
Our principal interests in evaluating the accn^rv ^f ,»,- ^ , 
aversion algorithm for the extractio^^'of^^Si'^^^^^^ ^ ^'^ine their relative performance in an 

denve aerosol opucal depths from curvcTv^usT^aSr^n. n ''''°" '"^thodoiogy emploj^d ^WSmoW 
sateUue data). In practice, this approadfSS Tdo^ZT 'f ^'^^ <'•«• ^ actual reflJcta^S^ve?Sm 
Pamn^representsa source of Sn?yo^ri^Tf Iht^^^^ °" ^°""' ""^"'=«- ^« -timaUonl/^S 

Figure 4.1 shows some rqirescntative nms of all a m«^.i- 
gcometnc and optical parameter.. I„ c^Z^i^lLT^l z^'^'tS^m^'C "' '^'''' '''''' ''' ^^'^ 
to the same phcnomenom observed by Kaufman * ' a. T ,^,^^° ^^"'^ ""<= ^ "«« a reverse in slope due 
albedos where the apparent reflectance is indcLdent ort^f" T"" ""'"'" " ''"^"'^'^ intermediate to the two 
pomt an inversion scheme for extracung opS deoU, ?' '"S °^"'^ '''^- '^^' « ^^ ^^^ '^ectance 
singulanty is however opUmal for extractinVrsfngSatSglS^^^^^ ''"^'^"^^ *'" "°' ^^-e^ge.^U 

models^ TTie speed of the 5S model pcr^itSe p?(LrcZJ» ^"f f P'" ''^'="'°" ^™" o" '^e 5S anTSSj 
to the DAVE model for the completeW^fSJput nS^^^^^ '^^ ^^, °i °P^'^ ^^"^ i"v"sion errors reladve 

mversion errors deduced from a number if LTSSf tS^s. '" '*"' ^■'- "^"^ '^^""^ *"« compared w^IJ 
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3.4 Apparent total reflectance differences 
(rebdve to the Dave model) computed for 
the Rayleigh and aerosol type 4 (AERO 4) 
atmospheres as a function of wavelength 
for a 4S degree solar zenith angle. The 
geometry was nadir viewing and the 
Hirface albedo was taken lo be zero. 




Figure 3.5 Apparent total reflectance dilTercnces 
(relative to the Dave model) computed 
for the Rayleigh and aerosol type 4 
(AERO 4) atmospheres as a function of 
surface albedo for a 45 degree solar 
zenith angle. The geometry was nadir 
viewiitg and the wavelength .SSS \ur. 



Figure 4.2 shows a rcprescniaUve set of the computed inversion errors for the 5S model as a function of various 
input parameters. The variation versus wavelength tends to decrease with incrcasuig wavelength although this was 
not always a consistant feature with the 5S inversions. TTw small peak at .655 jim reflects an equivalent (absorbing) 
moleculi optical depth disparity of about .1 between the DAVE and 5S codes. This disparity which was maximal at 
655 irni effectively demonstrates the level of error which can result from differences in molecular absorbing modcb 
if the gaseous contribution is interpreted as an aerosol optical depth variation . Note however that ihe d^t« 
between 5S and DAVE absorption models are fairiy extreme ; the corresponding differences between 5S and L6 arc 
considerably less. 

Variations in the correlation of inversion error magnitudes ( At^ ) with the apparent reflecumce errors { Ap* ) of 
Section 3 were most strongly influenced by the behavior of the slope dx^ /dp* in the neighbourhood of the 
Kaufman critical reflectance ( p^ ). The error magnitudes seen on the plot veisus albedo show an expected rapid rise 
on both sides of the critical surface albedo defined by Kaufman. Values of dXA / dp* at surface reficctance values 
removed from die critical reflectances were typically of the order of 10 or larger. From simple geometnc 
considerations (c.f . Figure 2 of ref. Kaufman)* ' an increase in surface albedo to a value p results in an increase of 
eiror by a factor of approximately [ - Pc / (P - Pc) 1- 
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(b) Same parameters as (a) for a surlace albedo of 0.4. 



i.v^" r^ '^-LTS^U^oS^raifS^irnS^^ r^*^^ -^ - -^ viewing geotnetry .^e 
angle tncrca^s the inverston error can i„crea« toTuesTLTr^e atTcT .T''^^" ^^ '^^ --^^-^ 
cxTDTs were, because of effecu related to the slotv. d^TZ ^ '^ "* ^^ ^-^^^ ^>- ^ DOM 

n^ladve behavior in the errors computed ft^e Lfu:^!' "nn'^r J"^' ^ '^^« ^ '^^ 5S results. However 
degrees solar zenith angle) were fairly well r^reZL^ Tt . 7 ^"^^^^ selected (nadir viewing and 45 
expected that the DOM inversion ^1t\2T^,!l'^'' ^^\T^' """^ ^^-^- For off nadir angles it .s 
the ratio of dp- enors found in Sec^"' "" ^^ '"°"' "'" ^ ^"^^ ^^ « f«^ i" rough prop^,S)n ^^ 



5. Si;MMARYANr>rnfvrrj^icf,^^|c^ 



evaluated in terms of leir iinpaci ^ an ftn^eJsl^ri^"'^/'™°^''"«- '^'^ ^o-^P^ted dispariues were Z 
measures of total apparent reflccS;«r P'"^""" ^°^ "'^^^^''^ 0P''=^' depth frori, ex<«,Jo^hS 
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pyure4.2 (A): Competed aerosol optical depth error of the 5S model relative to the Dave computations for the 
aerosol type 4 atmosphere (corresponding to xa = ^ as a funcuon of the cosine of the observers 
zenith angle. The wavelength is .455 |im. the solar zenith angle 45 degrees and the surface albedo 
zero. The values of the observes azimuth angle are noted above the figure. 

(B): Same as for (A) but for varying azimuth angle (observer's zenith angle noted above the figure). 

(Q: Aerosol optical depth error as a function of wavelength for a nadir viewing geometry. The solar 
zenith angle is 45 degrees and the surface albedo zero. 

(D): Varying surface albedo for the three aerosol types. The wavelength is .455 nm and the geometry 
is nadir viewing. Pc defines the critical ground reflectance (see text). 



The Turner and single scatter (Lowtran 6) models generally produced large errors in apparent reflectance which led 
to unacceptably large inversion errors. Overall the Turner model was not significanUy better than Lowrnn 6 model 
except at near nadir geometries and non zero albedos. This weakness of the Turner model implies that hybnd vernons 
which renormalize the Turner model to more accurate nadir results do not significanay improve the model away from 
the point of normalization. 

Although the DOM produced the most impressive results in tcnns of comparison with the Dave model and inversion 
errors ia time of execution (in particular for phase funcuons which are more anisotropic than those employed m this 
study) is a serious constrainL However the incorporaUon of this model in an interpolaUon algonthm which requires 
only a limited number of solutions per image (TeUlet et al.) ^ is a compromise approach which makes numerical 
solutions more competitive with analytical models such as 5S. Furthermore the single layer feature of the model we 
employed represents a significant savings in execution time relative to mulu-Iayer modcU. The excellent agre«ncnt 
which was obtained relative to the Dave model bears further invesUgation in terms of developing a bias correction 
algorithm which can be used to adjust single layer computations to mulu-layer results as a funcuon of aerosol 
altitude profile parameters. 
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in the neighbourhood of the critical fKaufmLwXJoSl. ^ inversion (toial apparent reflectance) is , at least 
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APPENDIX C 



Particle Growth Model Considerations 
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Appendix C: Particle Size Considerations 

A drop growth program was obtained from Richard Leaitch of the Cloud 
Physics branch Atmospheric Environment Service. It is designed to calcu- 
late the scattering phase function of an internally /externally mixed polydis- 
persion of aerosol particles comprising of up to 3 different types of particles. 
It was suggested by Leaitch to conduct a comparison of this model with 
the externally mixed mie model of Royer, described in the main text and 
refered to here as the CARTEL model. 

Changes made to AES program 

• Modes changed so that size distributions can be overlapping. 

• Modified gamma size distribution added (for mode 1). 

• Size ranges increased for comparison with CARTEL. 

• Cosmetic changes to use 3 element vectors instead of 3 variables for 
modes. 

Summary of program 

The following is a brief pseudo-code description of the drop growth 
scatting program. 

Start 

Read in the parameters 

• Thickness of layer 

• Wavelength 

• updraft acceleration (0) f 

• number of time steps (large) 

• Last time 



I 
I 

I 

• Starting temperature and pressxure (STP used) 

• Do dry and wet particle calculation ? (yes) I 

• Type of particle (chose up to 3 of the following) 
NaCl, iNH,)iSO,, NH^HSO^, HjSO^, NaNO^, H^O, Ca(7V0s)a 

• soluble mass fraction 

• Density of insoluble part 

• Index of refraction of insoluble part 

• Time increment 

• Bin sizes for particle size ranges 0.01-0.09, 0.1-0.9, 1.0-39.9,40-99 100- 
200 

• Total number of aerosols each type 

• Effective radius of aerosol each type 

• Effective variance of radius each type 

• Ratio of primary peak to secondary peak for bimodal distribution or 
if modified gamma size distribution used gamma. 

• Lower boundary of distribution (one for each mode) 

• Upper boundary of distribution (one for each mode) 

• Size distribution type : gamma, bimodal, log-normal, power-law, 
modified gamma 

• Updraft velocity(O) f 

• Thermal accomomdation coefficient (set to 1) f 

• Base condensation coefficient (set to 1) f 

• Starting relative humidity 

• Number of carbon groups in the fatty acid surfactant (19) f 



• Reference radius, (turns on/off surfactant effect) f 

• Molar fraction of sites not occupied by organic material during initial 
growth (set to 0.99) f 

fThcse parameters were not relevent to this study, values chosen to have 
no effect on calculations. 
Define particles 

Constants are stored for each type of particle (constants changed to 
agree with CARTEL input parameters see tables 2.1.1.1 and 2.1.3.1). 

Generate size distribution 

Use input size distribution, start, stop radii and effective radius and 
variance to generate size distribution. 

Beginning of calculation loop 

If dry pMticles skip growth loop first time through and goto phase 
function. 

Beginning of growth loop 

Grow particles based on relative humidity, and salt concentration. Two 
methods based on salt concentration (low-high) used. References given in 
code are Mason, 1971 and Pruppacher and Klett, 1978. 

Calculate phase function 

Mie calculations (see for example Wickramasinghe, 1973). 

Printout results 

Print size distribution, phase function and scattering parameters. 

Return to beginning of calculation loop until time > end time. 
End of program 

Commients and discussion 

The AES model uses effective mean radius and effective variance. The 
standard deviation and the geometric radius given by the CARTEL model 
are related to these quantities in the following manner. 

In Appendix A.l the log normal size distribution is given as : 



dN N -<lotr-l.,r.)» 



dlogr ■v/5irlog<r 
The same equation is expressed by Hansen and Travis (1974) as: 



Ar(r) = 



%/2^( 



»»i 



7r<7-r 



where the total number of particles has been normalized to 1. 
Observing that: 

dN = n{r)dr |3) 

dlogr = dr 

r 

and 

log r = log e In r 

we can compare the two equation and deduce that 

<Tg = In a (6) 

Tn = r, (7) 



The effective radius is defined by Hansen and Travis (1974) as 

/~r»n(r)rfr 



r, = 



J^r^n(r)dr 
This can be expressed as: 



s 



rf 



(8) 



__ , C4/37rr»n(r)rfr 

= 3^ (10) 



(11) 



This gives us the expression: 

Inr. = 31iir3-21iira (12) 

= 31iir„ + 9/21n'ff-21nr„-2In'«r (13) 

= liir„4-5/21nV (14) 

= \iir, + 5/2erl (15) 

Hansen and Travis (1974) give a relationship between the effective radius 
r,, effective variance v* and geometric radius r^ : 



r„ = 



' (l+r.)»/» 
Taking the log of both sides we obtain : 



(16) 



Inrp = Inr, -5/21n(H-re) (17) 

Inr. = lnr, + 5/21n(l+v,) (18) 

comparing with equation 15 we must have : 

<r^ =ln(H-v,) = ln'o- (19) 

To summarize the geometric radius r, and the standard deviation a are 
related to the effective radius r, and the effective variance w, as follows : 

V, = e'-^'-l (20) 

r, = r,(l+v,)«/' (21) 

These relationships were used to convert the input values of table 2.1.1.1 
for use in the AES program. 

It was diiHcult to determine what the exact comparison should be for 
the internally mixed and externally mixed aerosol models. While the ex- 
ternally mixed models have a volume ratio the AES programme needs a 
mass fraction (scduble to total). A calculation was performed to determine 
the mass fractions from the volume fraction given in table 2.1.3.1. 

Mf = y . ,Y' . (22) 



where 

Mf = mass fraction, soluble to total (23) 

V, = volume of soluble aerosol (24) 

dt = density of soluble aerosol (25) 

'inioi = volume of insoluble aerosol (26) 

*Hntoi = density of insoluble aerosol (27) 

Using 2.00 gm/cc as the density of soot and 1.80 gm/cc as the density of 
the soluble aerosol the mass fractions for the continental and urban internal 
mixtures were calculated and are shown in table C.l 



Continent al 

^- H- (^°) 0~ 50 70 80 90 95 

Vf 0.29/0.1 0.33072/0.0098 0.3533/0.0091 0.4495/0.0078 0.5715/0.006 6466/0 005 

Ml 0-963 0.966 0£72 0,981 0988 0.992 

Urban 



V> 0.61/0.22 0.6811/0.1771 0.7760/0.1264 0.8155/0.1041 0.8821/0.0665 0.9320/0 0384 
Ml Q-714 0-776 0847 0.876 0.923 0.956 



Table C.l: Mass Fractions 



Program validation 

A comparison was made using the AES program with Deirmendjian 
Haze L, Haze H, Haze M models (Deirmendjian, 1969) and SRA + water 
soluble model (Masuda and Takashima, 1988). These phase functions and 
values for the scattering coefficients were reproduced accurately. 

The AES program must go through several time loops to grow the 
aerosols. The number of time steps needed was studied so as not to run 
the program beyond the time needed for convergence to a solution. Table 
C.2 shows the variation of the value for the phase function at degrees as 
a function of time for a water soluble aerosol at 95 % relative humidity at 
550 nm. Setting the last time step to be 0.04 sec was felt to be adequate 
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for this study. 

Time step (s) Phase function (0 deg.) % change 

0.00 8.13689 

0.01 . 10.03803 18.9 

0.02 10.41801 3.6 

0.03 10.50699 0.8 

0.04 10.53174 0.2 

Table C.2 : Time for convergence 

Figure C.l shows the variation of the phase function with relative hu- 
midity for {NHt)3S04 at 550nm. 

As part of the validation procedure the single particle phase functions 
where calculated and compared with the CARTEL model resuHi (table 
C.3). This comparison showed the differences in the particle growth proce- 
dures in the two programs. Below the deliquescence point of {NH4)2S04, 
79% the results are very close, less than 2 % difference. Above the deliques- 
cence point where the particles will grow the results are quite different due 
to the change in the particle size distribution caused by droplet growth. 

Relative Humidity % AES Program CARTEL Model 

95 15.3770 10.61567 

«0 10.5317 7.95237 

M 6.9779 6.37056 

to 5.7313 5.81405 

m 4.9148 5.00754 

4.5576 4.57088 

Table C.3: Phase function at degrees 

A more relavent comparison might be at for angles in the backscat- 
tered direction. Table C.4 shows a comparison for the phase function at a 



scattering angle of 136 degrees. 



Relative Humidity % 


AES 


Program 


CARTEL Model 


95 




0.06774 


0.07681 


90 




0.07170 


0.08429 


80 




0.08546 


0.09536 


70 




0.10633 


0.10503 


50 




0.13755 


0.13276 







0.15814 


0.15889 



Table C.4: Phase function at 136 degrees 

Once again the diflFerences are greater for relative humidities greater 
than 80%, here the % difference is between 10-15 % . The scattering cross 
section is an important parameter and is compared in table C.5. There are 
significant differences for relative humidities > 80 %. 

Relative Humidity % AES Program CARTEL Model 



95 


5.5483E-9 


2.7542E-9 


90 


2.4019E-9 


1.5306E-9 


80 


9.6486E-10 


9.7266E-10 


70 


6.8763E-10 


8.0620E-10 


50 


5.5293E-10 


5.9871E.10 





5.0436E-10 


5.0170E-10 



Table C.5: Scattering cross section cm^ 
The single scattering albedo was also compared and proved to be less 
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sensitive to the difFcrcnces in growth model (table C.6). 

Relative Humidity % AES Program CARTEL Model 



95 


0.9987 


0.9932 


90 


0.9965 


0.9880 


m 


0.9884 


0.9813 


n 


0.9775 


0.9777 


m 


0.9666 


0.9703 





0.9603 


0.9649 



Table C.6: Single scattering albedo 

Particle growth rate was evaluated by using the change in size distribu- 
tion maximum as a function of relative humidity (Table C.7). Except for 
the smooth transition around the deliquescence point the results are the 
same as those found by Tang et al 1981 for a simileir study. 

Relative Humidity % Size dist. r^y maximum 10~* m v/tq 
95 



0.251 


2.22 


0.197 


1.74 


0.152 


1.35 


0.133 


1.17 


0.119 


1.05 


0.113 


1.00 



w 
m 



Table C.7: Size distribution number density maximum radius 

Discussion 

The initial reason for this study was to compare the effects of treating 
the aerosols as an internal or external mixture. The results suggest that 
the different growth models may be more significant for this particular 
aerosol model than internal /external treatments. Figure C.2 is variation 
with relative humidity for continental model at 550nm using table 2.1.1.1 
Leaitch model and figure C.3 is urban model. The dust component has not 
been included since the AES model did not produce the same results for 
a dust aerosol alone. The differences in the growth model may confuse a 
direct comparison of phase functions. Figure C.4 is the size distribution 
as a function of relative humidity. Table C.8 contains the values for the 
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forward and backseat tering directions. The volume scattering cross section 
is tabulated in table C.9. 

Relative Phase Function 136 degrees Phase Function degrees 
Humidity Cont. Urban Cont. Urban 



95 


0.06729 


0.0661 15.07792 14.77224 


90 


0.06995 


0.06578 10.45494 9.99415 


80 


0.08281 


0.09815 6.89182 6.41811 


70 


0.10933 


0.14186 5.63640 5.27573 


50 


0.14705 


0.17993 4.83736 4.58726 





0.16958 


0.18893 4.49158 4.35756 




Table C.8 Internally Mixed Results (No dust) 


Relative 


AES model (no dust) 


CARTEL model 


Humidity 


Cont, 
5.4182E-9 


Urban 


Cont. Urban 


95 


5.1946E.9 


2.7423E-9 2.16817E-9 


90 


2.3751E-9 


2.2028E-9 


1.5608E-9 1.26902E-9 


80 


9.540E-10 


8.742E-10 


1.0219E-9 8.5897E-10 


70 


6.684E-10 


6.234E-10 


8.612E.10 7.3665E-10 


50 


5.442E-10 


4.914E-10 


6.608E-10 5.8417E-10 





4.954E-10 


4.465E-10 


5.672E-10 5.1289E-10 



Table C.9 Scattering cross section cm^ 

References 

Deirmendjian, D., 1969. Electromagnetic Scattering on Spherical Polydis- 
persions. American Elsevier Pub. Co, Inc., New York, 290pp. 
Hansen, J. E. and L. D. Travis, 1974. Light scattering in planetary atmo- 
spheres. Space Set. Rev., 16:527-610. 

Mason, B. J., 1971. The Physics of Clouds. Clarendon Press, Oxford, 
671pp. 

Masuda, K. and T. Takashima, 1988. Dependence of the radiation just 
above and below the ocean surface on atmospheric and oceanic parameters. 
Applied Optics., 27:4891-4898. 

Pruppacher, H. R. and J. D. Klett, 1978. Microphysics of Clouds and 
Precipitation. D. Reidel Pub. Co., Boston, 714 pp. 
Tang, I., W. Wong and H. Munkelwitz, 1981. The relative importance 



of atmospheric suKates and nitrates in visibility reduction. Atmospheric 
Environment, 15:2463-2471. 

Wickramasinghe, N. C, 1973. Light Scattering Function* for Small Parti- 
chs, with Applications in Astronomy. Wiley, New York, 506pp. 



I 
I 



EFFiET or REUnVE HUMHXTY ON SCAnERNG 




9% 

"hi'i'" 

90% 
188 " 



.:r:::;r:^-j^s^a 



MNUOIB 

MOC REPORT FIG 
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Figure C.2 Internal Mixture, continental, 550 nm. 
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Figure C.3 Internal mixture, urban, 550nm. 
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Figure C.A Particle growth. 
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